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TIE 741-MED-3042 


SUMMARY REPORT 
ON 

A BASIC MODEL OF CIRCULATORY, FLUID, 
AND ELECTROLYTE REGULATION IN 
THE HUMAN SYSTEM BASED UPON 
THE MODEL OF GUYTON 

DEVELOPED BY 
RONALD J. WHITE, Ph.D. 



INTRODUCTION 


The basic model of Guyton and co-workers is a model of circulatory, 
fluid, and electrolyte regulation. The model is functional in nature and is 
based almost entirely on experimental data and cumulative present knowledge 
of the many facets of the circulatory, fluid, and electrolyte regulatory systems 
of the human body. 

The attached Study Report (attachment 1) provides a detailed description 
of Guytons' model and modifications developed by Dr. Ronald White. Also 
included in the study report are descriptions of several typical experiments 
which the model can simulate to illustrate the model's general utility. Chapter 
IV of the study report includes a discussion of the problems associated with the 
interfacing of the model to other models such as respiratory and thermal regu- 
lation models which is of prime importance since these stimuli are not present 
in the current model. 

Attachment 2 (TIR 74I-MED-3017) provides a user's guide for the opera- 
tion of the model on the Xerox Sigma 3 computer. Two programs are described 
in the user's guide. Model A is the basic Guyton model and Model B is Dr. 

Ronald White's version of the Guyton model. 

Attachment 3 (TIR 741-MED-3026) presents a verification plan and pro- 
cedure for performing experiments with the model. 

MODEL DESCRIPTION 

The model consists of 16 distinct subroutines concerned with physiological 
function, and contains almost 100 independent parameters as well as more than 
350 mathematical relations, (see figure 1). Each function has only been modeled 
in a crude way with little attention being given to fine details. The systems 
analysis thus developed is successful in predicting the outcome of many varied 
stress experiments. This is only possible because of the extreme stability and 
built-in compensations of the body's actual circulatory system. 
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The model may be viewed as a controlled system plus controlling system 
with the controlling system having three major components: local control, hor- 
monal control, and autonomic control. These controls act to drive the controlled 
system to the appropriate level in response to stress. There are no thermal 
regulatory components present in either the controlled or controlling system. 
Respiratory elements remain to be added with the exception of the effect of 

pulmonary interstitial fluid on aortic oxygen saturation. Future plans include 

+ + 

the addition of hydrogen ion considerations. Only the major cations, Na and K , 
are treated presently. The model may be classified as an intermediate to long- 
term model with simulations of the order of days or weeks being the primary 
concern, although short-term simulations, as in the exercise experiment, may 
be conducted. 



ATTACHMENT I 


A BASIC MODEL OF CIRCULATORY, FLUID, AND 
ELECTROLYTE REGULATION IN THE HUMAN SYSTEM 


-Study Report- 


by 


Ronald J. White, Ph.D. 
Mathematics Department 
University of Southwestern Louisiana 
Lafayette, Louisiana 70501 



i. Introduction 


In recent years, it has become increasingly evident that systems analysis and control 
theory offer a convenient path to the goal of understanding the functioning and interrela- 
tion between various parts of complex physiological systems. For a system as large and 
complex as the human body, the systems analysis may be broken down into several large 
interacting subsystems and each subsystem may be treated somewhat independently of 
the others. These large subsystems, each with many subsystems of its own, may then be 
combined to produce a model of the overall functioning of the human body. 

This study report considers only a model of circulatory, fluid, and electrolyte regulation 
developed recently by Guyton, Coleman, and Granger (1). Other models of fluid and 
electrolyte regulation are briefly summarized elsewhere (2). The model of Guyton and co- 
workers is functional in nature and is based almost entirely on experimental data and cumulative 
present knowledge of the many facets of the circulatory, fluid, and electrolyte regulatory 
systems . 

The model itself consists of 16 distinct subroutines concerned with physiological function, 
and contains almost 100 independent parameters as well as more than 350 mathematical 
relations. Each function has only been modelled in a crude way with little attention being 
given to fine details. The systems analysis thus developed is successful in predicting the 
outcome of many varied stress experiments. Apparently, this is only possible because 
of the extreme stability and many built-in compensations of the actual circulatory system. 

A flowchart showing the interconnection of the basic subroutines is given in Figure 1. 

The subroutines PUTIN and PUTOUT are input-output routines and are not discussed here. 

(The use of this model is discussed elsewhere (3), (4).) The following chapter considers 
each of the physiological routines in detail, while Chapter HI presents results obtained 
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from typical experiments. Chapter IV contains a brief summary of the model characteristics 
as a whole and discusses the problem of interfacing 'the Guyton model with appropriate 
respiratory and thermal models. 

II. The Model 

The model of Guyton and co-workers will be examined subroutine by subroutine. For 
each subroutine, a color-coded flow chart is given and a line-by-line description is.pre- 
sented.* Red indicates a variable inputfrom another subroutine, blue indicates a variable 
output to another subroutine, green indicates a variable never calculated (independent), 
and black indicates a variable used only in the subroutine being considered. A complete 
glossary of terms is presented in Appendix A. Units used are: volume in liters, mass in 
grams, time in minutes, chemical units in mifliequivalents, pressure in millimeters of 
mercury, and control factors as ratio to normal. 


*This description was provided by Dr. A. C. Guyton. 



SUBROUTINE MEMO 
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Program Listing : See Program 1 . 

Flowchart: See Figure 2. 

Line 15. Addition of blood volume (VP+VRC) and subtraction of volumes of blood in all 
portions of the systemic circulation (VVS, VAS, VLA, VPA, VRA) to yield 
the net difference between blood volume and volume calculated in all the 
capacitative reservoirs of the systemic circulation; output of this line represents 
the correlation factor (VBD) that is added to the flow of the systemic circulation 
into the small veins, thus bringing the volume of blood in the circulation up to 
the appropriate level. This allows updating of blood volume when volumes 
pass through the capillary walls, when volume is gained by the process of drinking, 
or lost through the kidneys, and so forth. 

Line 16. Integration of rate of blood flow into the veins (DVS) plus correction factor 
(VBD) gives volume of blood in the veins of the systemic circulation (WS). 

Line 17. Integration of rate of change of volume in the pulmonary arteries (DPA) plus 
correction factory (VBD) gives the instantaneous volume in the pulmonary 
arteries (VPA). 

Line 18. Integration of rate of change of volume in systemic arteries (DAS) plus correction 
factor (VBD) gives actual volume in systemic arteries (VAS). 

Line 19. Integration of rate of change of volume in left atrium and pulmonary veins (DLA) 
plus correction factor (VBD) gives instantaneous volume in left atrium and pul- 
monary veins (VLA). 

Line 20. Integration of rate of change of volume in right atrium (DRA) plus correction 
factor (VBD) gives volume of blood in right atrium (VRA). 

Line 21 . Volume in systemic arteries (VAS) minus constant gives excess volume in systemic 
arteries (VAE) that causes stretch of the arterial walls. 

Line 22. Excess volume in systemic arteries (VAE) divided by compliance of the systemic 
arteries gives arterial pressure (PA). 

Line 23. Factor of 100 divided by arterial pressure (PA) gives arterial pressure multiplier 
factor for alteration of peripheral resistance caused by stretching of arteries 
resulting from arterial pressure (PAM). 

Line 24. Effect of autonomic stimulation (AUH) on loading effect of systemic arterial 
pressure (PA) to give effective arterial pressure on left ventricle (PA2). 
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Line 25. Function curve (See Figure 3.) showing effect of systemic arterial pressure (PA2) 
in loading left ventricle and determining its pumping effectiveness (LVM). 

Line 26. Volume of blood in right atrium (VRA) minus constant gives excess volume of 
blood in right atrium (VRE) that causes stretching of right atrium. 

Line 27. Excess volume of blood in right atrium (VRE) divided by compliance of right 
atrium gives right atrial pressure (PRA). 

Line 28. Curve (See Figure 4.) relating right atrial pressure (PRA) to output of right 
atrium under normal operating conditions of right atrium (QRN). 

Line 29. Volume in the pulmonary arteries (VPA) minus a constant factor gives the excess 
volume in the pulmonary arteries that causes stretch of the arteries (VPE). 

Line 30. Excess volume in the pulmonary arteries (VPE) divided by the capacitance of 
the pulmonary arteries gives the pulmonary arterial pressure (PPA). 

Line 31-33. Curve fitting process to empirically calculate resistance in pulmonary arteries 
to the midpoint of the pulmonary capillaries (RPA) from the pulmonary arterial 
pressure (PPA). 

Line 34. Calculation of the effect of autonomic stimulation (AUH) on the degree of 

loading of the right ventricle (PP2) caused by pulmonary arterial pressure (PPA). 

Line 35. Curve (See Figure 15.) relating effective pulmonary arterial pressure (PP2), and 
pumping effectiveness of the right ventricle (RVM). 

Line 36. Volume of blood in pulmonary veins and left atrium (VLA) minus constant factor 
gives excess volume (VLE) causing stretch of left atrium and pulmonary veins. 

Line 37. Excess volume in left atrium and pulmonary veins (VLE) divided by capacitance 
of left atrium and pulmonary veins gives pressure in the left atrium (PLA). 

Line 38. Curve (See Figure 6.) giving normal output o,f left ventricle (QLN) for each 
given value of pulmonary left atrial pressure (PLA). 

Line 39. Curve fitting process based primarily on waterfall effect to calculate resistance 
of pulmonary veins (RPV) whose change depends primarily on level of left 
atrial pressure (PLA). 

Line 40. Calculation of total pulmonary resistance (RPT) by adding pulmonary arterial 
resistance to midpoint of capillaries (RPA) and pulmonary venous resistance 
from midpoint of capillaries to left atrium (RPV). 
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Line 41. Pulmonary arterial pressure (PPA) minus left atrial pressure (PLA) gives the 
pressure gradient through the lungs (PGL). 

Line 42. Pressure gradient through the lungs (PGL) divided by resistance of the pulmonary 
circuit (RPT) gives rate of blood flow into the pulmonary veins and left atrium 
(QPO). 

Lines 43-44. Calculation of vasoconstrictor factor caused by angiotensin (ANU) (but does not 
fall below 0.8. ). 

V 

Line 45. Volume of blood in the veins (WS) minus vblume'dfSIood at zero venous 

pressure (VVR) minus vasoconstrictor effect of angiotensin (ANU) gives excess 
venous volume before correction factor for stress relaxation (VVE). 

Lines 46-47. Excess volume of blood in the circulation (WE) minus stress relaxation 
factor 0/V7) gives excess volume of blood in the systemic veins after stress 
relaxation factor correction (W8)(not allowed to fall below 0.0001). 

Line 48. Excess systemic venous volume (W8) divided by capacitance the veins (CV) 
gives pressure in the veFns (PVS). 

Lines 49-50. Right atrial pressure (PRA) but never negative (PR1). 

t 

Line 51 . Calculation of resistance between veins and right atrium (RVG) as determined 
by the level of small vein venous pressure (PVS). 

Line 52. Pressure gradient from the veins to the right atrium (PVS-PR1) divided by the 

large vein resistance (RVG) gives rate of blood flow into the right atrium (QVO). 

Lines 53-55. Curve fitting process to give effect of changing capillary pressure (PC) and 
autonomic stimulation (AVE) on venous resistance (RVS), showing principally 
a partial waterfall effect, the constancy of the pressure of which is determined 
by the constant CN7. 

Arterial pressure (PA) minus pressure in the small veins (PVS) gives pressure 
gradient of the systemic circulation (PGS). 

Addition of the resistance from the aorta to the midpoint of the capillaries 
to the resistance from the midpoint of the capillaries to veins to give the 
total resistance of the non-muscle, non-renal portion of the systemic circulation 
(RSN). 

Pressure gradient in the systemic circulation (PGS) divided by the resistance 
in the non-muscle, non-renal circulation (RSN) gives blood flow in the non- 
muscle, non-renal circulation (BFN). 


Line 56. 
Line 57. 

Line 58. 
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Line 59. 
Line 60. 

Line 61. 
Line 62. 

Line 63. 

Line 64. 

Line 65. 
Line 66. 

Line 67. 

Line 68. 

Line 69. 


Calculations of the resistance through the muscle circulation of the body (R5M). 

Pressure gradient in the systemic circulation (PGS) divided by resistance in the 
muscle circulation (RSM) gives blood flow in the muscle circulation (BFM). 

Addition of blood flow in the non-muscle, non-renal circulation (BFN) plus 
muscle blood flow (BFMj plus renal blood flow (RBF) plus A-V fistula flow (FIS) 
gives blood flow from aorta through the systemic circulation (QAO). 

Calculation of actual output of left ventricle (QLO) based on the following 
factors: Output of left ventricle under normal conditions (QLN), effect of 
arterial pressure loading factor on left ventricle (LVN), basic strength of left 
ventricle (HSL), degree of autonomic stimulation of left ventricle (AUH), degree 
of deterioration of left ventricle caused by low coronary blood flow (HMD), and 
degree of hypertrophy of the left ventricle (HPL), 

Calculation of actual output of right ventricle (QRO) similar to line 62, except 
that left ventricular pumping (QLO) plays a part. 

Stability check on calculation of rate of blood flow info the pulmonary veins 
and left atrium (QPO). 

Stability check on calculation of rate of blood' flow into right atrium (QVO). 

Blood flow through the systemic circulation (QAO) minus blood flow out of the 
small veins into the atria (QVO) gives rate of change of volume in the systemic 
veins (DVS). 

Output of right ventricle into pulmonary arteries (QRO) minus rate of blood flow 
from the pulmonary arteries into the pulmonary veins and left atrium (QPO) gives 
rate of change of volume in the pulmonary arteries (DPS). 

Actual rate of output of left ventricle (QLO) minus rate of blood flow from 
systemic arteries through systemic circulation (QAO) gives rate of change of 
blood volume in systemic arteries (DAS). 

Rate of blood flow into the pulmonary veins and left atrium (QPO) minus rate 
of blood flow out of the pulmonary veins and left atrium (QLO) gives the rate 
of change of volume in the left atrium and pulmonary veins (DLA). 

Rate of blood flow into the right atrium (QVO) minus rate of blood flow out of 
right atrium (QRO) gives rate of change of blood volume in right atrium (DRA). 


Line 70. 
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SUBROUTINE AUTQ 

Program Listing : See Program 2. 

Flowchart : See Figure 7. 

Line 9. Calculation of the bias on the setting of the autonomic drives in the central 
nervous system (EXE) due to the degree of exercise (EXC) and muscle Pq 2 
(P20). The factor EXE is normally zero, but increases with degree of 
exercise or decrease in muscle P20. The effect of muscle P20 is presumably 
mediated through such factors as release of lactic acid and pH, CO 2 , and 
P20 changes in blood carried from muscles to chemoreceptive areas. 

Lines 10-12. Calculation of the effect of tissue PQ 2 (POQ), limited between values of 
4 and 8, for determining the factor that drives autonomic responses, assuming 
that the tissue PQ 2 biases the setting of the effect of pressure on the central 
nervous system autonomic feedbacks. The Pq„ effect acts through direct 
effect of P 02 on fbe vasomotor center, througTi associated effects of CO 2 
that go along with P 02 changes, through possible cardiac receptors and other 
peripheral receptors that may be related to tissue blood flow and tissue 

Line 13. Calculation of (PA1 ) both the effect of exercise and arterial pressure (PA) and 
tissue PO 2 (POQ) that cause biasing of the factor for control of autonomic 
outputs. 

Lines 14-16. Effect of drive factor on the vasomotor center (PA1) of the autonomic system 
caused by pressure effects operating indirectly through the chemoreceptors 
(AUC). The function is expressed algebraically with two break points in the 
curve, at 80 and at 40. The autonomic output is expressed in terms of positive 
sympathetic drive and negative parasympathetic drive. 

Lines 17-19. Similar function as in line 14-16, but this time resprsenting pressure effect 
operating through baroreceptors to stimulate the autonomic system. Output 
(AUB) represents positive sympathetic drive and negative parasympathetic drive. 

Line 20. Adjustment of sensitivity of baroreceptor drive with output A1B. 

Lines 21-23. Similar function as in line 14-16 or line 17-19, but this time for central 
nervous ischemic response, with output AUN. 

Lines 24-25. These lines plus line 27 of SUBROUTINE MISC1 allow for adaptation of 
the baroreceptor system. The time constant for adaptation is determined by 
AUK. AU6 always reapproaches the value 1 with time because of adaptation. 
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Line 26. Summation of autonomic stimulation from cHemorec'eptors (AUC), baroreceptors 
(AU6), and CMS ischemic response (AUM)to give the final equilibrium 
summated effect that will be approached (DAU). 

Lines 27-29. Time delay circuit for full realization of autonomic drive. The output 
AUJ approaches the final equilibrium (DAU) with time constant determined 
by Z8. 

Lines 30-32. Calculation of the overall activity of autonomic system (AU) which represents 
the tendency to increase overall functional activity of the heart and to increase 
vascular constriction throughout the body. 

Line 33. Allows pre-set value (STA) to be substituted for AU. 

Line 34 Calculates departure of overall activity (AU) from normal (AUO). 

Line 35. Calculation of autonomic drive for peripheral circulation (AUP) from sensitivity 

control (AUQ) and autonomic level (AUO). 

Line 36. Same as line 35, except for heart (AUH). 

Line 37. Same as line 35, except for heart rate (AUR). 

Line 38. Sets sensitivity for control of systemic venous vascular volume (WR). 

Line 39. Determines sensitivity of autonomic drive to control arteriolar resistance in the 
muscle and non-muscle portions of the circulation, and also to control the 
degree of stimulation of the afferent arterioles of the kidneys (AUM). 
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SUBROUTINE HORMON 


Program Listing: See Program 3. 

Flowchart: See Figure 8. 

Lines 12-13. Determination of effect of the ratio of extracellular fluid potassium (CKE) 

to sodium concentration (CNA) by means of curve fitting on control of aldosterone 
secretion (AMR). 

Function curve (See Figure 9.) to determine' the effect of arterial pressure (PA) on 
aldosterone secretion (AMP). 

Calculation of total control effect on aldosterone secretion (AMI) by multiplying 
effect of potassium of sodium ratio (AMR), pressure (AMP), and stimulatory effect 
of angiotensin (ANM). 

Decay effect which specifies rate of buildup of aldosterone in the interstitial 
fluids. This level approaches the level set by the aldosterone control (AMI) 
with time constant AMT. The output is AMC, the concentration of aldosterone 
expressed as the ratio of the concentration to the normal value. 

Calculation of the degree of effect of aldosterone (AM) from the aldosterone 
concentration by empirical means. 

Lines 23-24. Subtraction of concentration of sodium in extracellular fluids (CNA) 

from a constant to give sodium concentration factor for control of angiotension 
secretion (CNE), a factor never less than I. 

Line 25. Determination of effect of glomerular filtration rate (@RN)., extracellular sodium 
concentration (CNA), and degree of normality of kidneys (REK) on renin 
output and subsequent formation of angiotensin (ANR). 

Lines 26-30. Curve fitting technique which allows for effect of renal blood flow and sodium 
level on angiotensin formation (ANP). 

Line 31 . Mathematical technique to allow for stability control of total factor for angiotensin 
control (AN1). 

Line 32. Decay effect which allows for a buildup of angiotensin in the circulation. 

The time constant for the delay is ANT and the concentration of angiotensin is ANC. 

Lines 33-34. Calculation of the degree of effect of angiotensin (ANM) from angiotensin 
concentration (ANC). ANM has a lower limit of 0.7. 


Line 14. 
Line 15. 

Line 16. 
Line 17. 
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- SUBROUTINE BLOOD 
Program Listing : See Program 4. 

Flowchart: See Figure 10. 

Line 9. Calculation of blood volume (VB) by adding plasma volume (VP) and red cell 
volume (VRC). 

Calculation of hematocrit (HM) by dividing red cell' volume (VRC) by blood 
volume (VB) and multiplying by 100. 

Calculation of the actual viscosity of blood caused by hematocrit (VIE). 

HMK and HKM are constants. 

Calculation of total relative viscosity of blood (ViB), assuming viscosity of 
wafer to equal one, by adding viscosity caused by red cells (VIE) to a constant 
factor representing viscosity of plasma. 

Calculation of viscosity multiplier (VIM) by multiplying relative viscosity times 
a constant. This factor is the viscosity multiplier factor that determines 
relative changes in vascular resistance with changes in viscosity from normal. 

Calculates rate of red cell destruction (RC2) by multiplying volume of red 
cells in circulation (VRC) times a constant (RKC). 

Lines 18-20. Calculation of the effect of non-muscle tissue P02 (POT) as a drive in 

causing formation of red blood cells. The drive is considered to be zero when 
POT equals the constant factor POl and to increase as the tissue PO 2 f Q lt s below 
this value. The drive factor (P02) has a minimum value. POY determines the 
sensitivity of the circuit and RC1 is the rate of red cell production. 

Line 21 . Calculation of net rate of change of red cell volume (RCD) from red cell 
production rate (RC1) and red cell destruction rate (RC2). 

Line 22. Calculation of volume of red cells in circulation (VRC) by integration. 


Line 10. 
Line 1 1 . 
Line 12. 

Line 13. 

Line 17. 



SUBROUTINE MUSCLE. 


. 1 ? 


Program Listing ; See Program 5. 

Flowchart: See Figure 11. 

Line 8. Calculation of aortic arterial oxygen saturation (OSA) by subtracting the fraction 
of desaturation of arterial blood (calculated by multiplying quantity of free 
fluid in lungs (VPF) times a constant) from maximum saturation (ALO) of one. 

Line 9. Calculation of volume of oxygen per liter of aortic arterial blood (OVA) by 
multiplying arterial oxygen saturation (OSA) times hematocrit (HM) times' 
a constant. 

Line 10. Calculation of actual venous oxygen saturation (OVS) by means of a delay 

mechanism which allows the venous saturation to rise to Its equilibrium value. 

The product of muscle blood flow (BFM) and volume of oxygen per liter (OVA) 
is the rate of delivery of oxygen to the muscle cells. Subtracting the rate of oxygen 
utilization by the tissues (RMO) gives the rate of oxygen delivery to the veins. 
Dividing this by the muscle blood flow (BFM), hematocrit (HM), and a constant 
yields the venous oxygen saturation after equilibrium ’has been established. 

The rest is the delay mechanism with the constant Z 6 controlling the time constant. 

Line 11. Calculation of venous oxygen pressure (PVO) by multiplying venous oxygen 
saturation (OVS) times a constant. 

Line 12. Calculation of rate of oxygen delivery to muscle cells from capillaries (RMO). 

The assumption is made that oxygen in the muscle capillaries is equal to oxygen 
in the venous blood (PVO). The pressure difference between oxygen in the 
muscle capillaries (PVO) and oxygen in the muscle cells (PMO) is multiplied 
by a constant (PM5) and divided by a resistance factor determined by the number 
of capillaries that are open at any given time. 

Line 13. Calculation of total quantity of oxygen stored in cells (QOM) by integration 
of rate of change of oxygen in muscle cells (RMO-MMO). The exponential 
factor represents a damped integration. Note that QOM represents oxygen 
stored in all of its energy forms, including dissolved oxygen, oxygen bound 
with myoglobin, and oxygen equivalents of energy compounds such as ATP and 
creati ne phosphate . 

Lines 14-16. Calculation of pressure of oxygen in muscle cells (PMO or PM1) by a curve 
fitting process from the quantify of oxygen in the cells (QOM). Note that PM1 
cannot fall below 0.001 mm Hg. 
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Lines 17-18. Calculation of muscle cell oxygen pressure effective in depressing rate of 
metabolism in cell (P20). 

Line 19. Calculation of the effect of the degree of autonomic stimulation (AUP) on 

the rate of metabolism expressed as an autonomic multiplier effect on metabolism 
(AOM). This value is normally unity. The sensitivity of the effect of autonomic 
stimulation on metabolism is determined by the constant 02A. 

Line 20. Calculation of rate of utilization of oxygen by the cells (MMO) by multiplying 

autonomic effect (AOM)/ exercise effect (EXC), basal level of oxygen utilization 
(OMM)/ and the effect of decrease in muscle cell Poo* This last effect is a 
curve fitting process involving P20 that assumes that the oxygen level must fail 
nearly to zero before very significant decrease in the rate of metabolism occurs. 

Line 21 . Calculation of difference between capillary ? 02 ' assuming this equals venous 
P02 (PVO), and normal capillary P02°^ 40. 

Lines 22-23. Calculation of sensitivity control for oxygen feedback loop (POE). The 
constant POM determines the degree of sensitivity. 

Line 24. Calculation of the autoregulation multiplier for the muscle vascular circuit 
(AMM) by a time delay mechanism. The constant A4K is the time constant 
for this delay. 
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SUBROUTINE AUTQRG 

Program Listing: See Program 6. 

Flowchart-: See Figure 12. 

Line 11. Calculation of the actual venous oxygen saturation (OSV) based on a delayed 
approach to its equilibrium value. The product of the volume of oxygen in 
each liter of arterial blood (OVA) and blood flow to the non-muscle tissues 
(BFN) gives the rate of transport of oxygen by arteries to the non-muscle tissues. 

By subtracting the rate of oxygen utilization of the tissues (DOB), the rate of 
oxygen delivery to the veins of the non-muscle tissue is obtained. This difference 
divided by the blood flow to the tissues (BFN), hematocrit, and a constant 
yields the equilibrium venous oxygen saturation. The time constant is Z7. 

Line 12. Calculation of venous oxygen PQ 2 (POV) from venous oxygen saturation (OSV). 

Lines 13-14. Calculation of resistance of diffusions of oxygen from capillaries to cells 

(RDO) assuming that far greater numbers of capillaries open up and the resistance 
decreases as the tissue P 02 (POT) falls below normal. 

Calculation of rate of delivery of oxygen from capillaries to tissue cells (DOB) by 
multiplying pressure difference between pressure of oxygen in tissue capillaries 
(assumed equal to POV) and pressure of oxygen in the tissue cells (POT) times 
a constant and dividing by the resistance for diffusion of oxygen (RDO). 

Calculation of the rate of oxygen utilization by cells (M02) by multiplying basal 
level of oxygen utilization (02M) by the autonomic stimulatory effect (AOM) and the 
tissue PO 2 e ff® c t on oxygen utilization. 

Calculation of actual total quantity of oxygen accumulated in the cells (Q02) 
by integration of the rate of accumulation of oxygen in the tissue cells. This rate 
is determined by subtracting the rate of utilization of oxygen in the cells (M02) 
from rate of delivery of oxygen to cells (DOB). 

Calculation of tissue cell Pq^ (POT) from quantity of oxygen accumulated in 
the cells (Q02). 

Lines 19-20. Calculation of .effective tissue Pq 2 for oxygen utilization (PIO). 

Line 21. Calculation of pressure difference that acts as control factor for autoregulation 
of non-muscle blood flow (POD) by subtracting reference value (POR) from 
capillary Pq ? in non-muscle tissues (assumed to equal POV). 


Line 15. 


Line 16. 


Line 17. 


Line 18. 
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Lines 22-23. Calculation of rate of change of rapid autoregulation vasoconstrictor 

effect (POB) from pressure difference POD. The sensitivity is set by POK. 

Line 24. Calculation of rapid autoregulation multiplier factor (ART) by time delay 
mechanism with time constant AIK. 

Line 25. Multiplication of the three autoregulation factors, short-time (ART), intermediate- 
time (AR2), and long-time (AR3), to give the total autoregulation factor (ARM) 
which multiplies the basic resistance for blood flow through the non-renal sector of the 
circulation. 

Lines 29-30. Calculation of rate of change of intermediate autoregulation vasoconstrictor 
effect (POA) from pressure difference POD. The sensitivity is set by PON. 

Line 31. Calculation of intermediate autoregulation multiplier factor (AR2) by time delay 
mechanism with time constant A2K. 

Line 35. Branching step depending on whether POD is positive or negative. 

Lines 36-37. Calculation of rate of change of long-time autoregulation vasoconstrictor 
effect (POC) if POD is positive. The sensitivity is set by POZ. 

Line 38. Calculation of rate of change of long-time autoregulation vasoconstrictor effect 
(POC) if POD is negative. The sensitivity is set by POZ. 

Line 39. Minimum value set for POC. 

Line 40. Calculation of long-time autoregulation multiplier factor (AR3) by time delay 
mechanism with time constant A3K. 
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SUBROUTINE ADH 


Program Listing: See Program 7. 

Flowchart: See Figure 13. 

Line 7. Calculation of the effect of extracellular fluid osmolarity on antidiuretic 
hormone secretion (CNB) by subtracting a reference value (CNR) from the 
concentration of sodium in the extracellular fluids (CNA) (taken to be a measure 
of the osmolarity of extracellular fluids). 

Line 8. Calculation of the partial effect of right atrial pressure (PRA) in controlling 
antidiuretic hormone secretion (AHZ). 

Line 9. Calculation of the degree of adaptation of the right atrial pressure mechanism 
for affecting antidiuretic hormone secretion (AHY) by time delay mechanism. 

Lines 10-11. Calculation of the effect of autonomic stimulation on the rate of antidiuretic 
hormone secretion (AH8) from the autonomic multiplier AUP. 

Line 12. Prevents sodium factor CNB (line 7) from being negative. 

Lines 13-14. Calculation of the equilibrium control value of antidiuretic hormone secretion 
(AH) by summation of the factors that cause antidiuretic hormone secretion. 

These factors are the sodium factor (CNB), the right atrial pressure factor 
(AHZ and AHY), and the autonomic factor (AH8). 

Line 15. Calculation of rate of buildup of antidiuretic hormone concentration (AHC) 

in the body fluids by time delay mechanism with time constant AHK. Normally 
AHC = 1. 

Lines 16-17. Calculation of antidiuretic hormone multiplier (AHM) by a curve fitting 
process from the antidiuretic hormone concentration (AHC). 
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SUBROUTINE MISC1 


Program Listing: See Program 8. 

Flowchart: See Figure 14. 

Line 10, Calculation of the rate of progression of stress relaxation (VV6) by subtracting 
from excess systemic venous volume (SR*WE) the reference volume (,301*SR) 
and the actual degree of stress relaxation (VV7). The factor SR is the adjustable 
intensity of stress relaxation. 

Line 11 , Calculates the actual stress relaxation volume (W7) by integration with time 
constant SRK. The amount of stress relaxation in the systemic veins is set to 
be somewhat higher than the normal stress relaxation of the venous system to 
make up for the fact that similar stress relaxation factors are not calculated 
for other parts of the circulation. 

Lines 17-19. Calculation of rate of water intake (TVD) by multiplying the tissue perfusion 
factor for thirst stimulation (STH) by the thirst center drive. The thirst center 
drive is calculated from the antidiuretic hormone multiplier (AHM) by curve 
fitting under the assumption that the same factors that drive antidiuretic hormone 
secretion play a similar role in causing thirst. 

Line 20, Calculation of total body water (VTW) by adding extracellular fluid volume (VEC) 
to intracellular fluid volume (VIC). 

Line 27, Calculation of effective adaptation of baroreceptor system (AU4) when coupled 
to lines 24-25 of SUBROUTINE AUTO. 



SUBROUTINE HEART 
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Program Listing: See Program 9. 

Flowchart: See Figure 15. 

Line 10. Calculation of rate of deterioration of the heart (DHM) by a curve fitting process 
that assumes the deterioration increases progressively as tissue PO 2 (POT) falls 
below 6 mm Hg. 

Lines 11-12. Calculation of deterioration multiplier factor (HMD)-^hich multiplies 
the strength of the two ventricles by integrating the rate DHM. 

Line 16. Calculation of mean circulatory pressure (PMC) by adding excess blood volume 
in systemic arteries (VAE), excess venous vascular volume (WE), excess 
volume in right atrium (VRE), excess volume in pulmonary arteries (VPE), and 
excess volume in left atrium (VLE), and dividing by a constant. 

Line 17. Calculation of the mean systemic pressure (PMS) by adding excess blood volume 
in systemic arteries (VAE), excess volume in venous system (WE), and excess 
volume in right atrium (VRE), and dividing by a constant. 

Line 18. Calculation of the mean pulmonary pressure (PMP) by adding excess volume 
in pulmonary arteries (VPE), and excess volume in left atrium (VLE), and 
dividing by a constant. 

Line 24, Calculation of heart rate (HR) by multiplying the summation of a basic heart 
rate factor (a constant), the reflex effect due to right atrial pressure (PRA), 
and the autonomic drive effect (AUR) by the effect of cardiac deterioration 
expressed in terms of the degree of normality of the heart (HMD). 

Line 25, Calcualtion of total peripheral resistance (RTP) by dividing the difference 

between aortic pressure (PA) and right atrial pressure (PRA) by the blood flow 
in the systemic arterial system (QAO). 

Line 26. Calculation of the stroke volume (SVO) by dividing the cardiac output (QLO) by the 
heart rate (HR). 
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SUBROUTINE CAPMBD 

Program Listing ; See Program 10. 

Flowchart: See Figure 16. 

Line 9, Calculation of total tissue pressure {PTT) from total volume of fluid in the interstitial 
compartment (VTS). 

Line 10. Calculation of free fluid in the interstitial Spaces (VIF)&y subtracting volume 
of gel fluid (VG) from total interstitial fluid volume (VTf>). 

Line 11. Calculation of solid tissue pressure (PTS) from volume of free interstitial fluid 
(VIF) by graphical interpolation (See Figure 17.). 

Line 12. Calculation of pressure of free interstitial fluid (PIF) by subtracting solid tissue 
pressure (PTS) from total tissue pressure (PTT). 

Line 13. Calculation of concentration of protein in free interstitial fluid (CPI) by dividing 
quantity of protein (IFP) by .volume of free fluid (VIF'j. 

Line 14. Calculation of colloid osmotic pressure of free interstitial fluid (PTC)by 

multiplying concentration of protein in free fluid (CPI) by a constant factor. 

Line 15. Calculation of concentration of proteins in plasma (CPP) by dividing quantity 
of protein in the plasma (PRP) by the plasma volume (VP). 

Line 16. Calculation of colloid osmotic pressure of plasma (PPC) by multiplying the 
concentration of plasma protein (CPP) by a constant. 

Line 17. Calculation of pressure gradient from midpoint of the capillaries to the veins 

(PVG) by multiplying resistance of the veins (RVS) times the blood flow through 
the non-renal, non-muscular portions of the circulation (BFN) and times a constant 
to account for blood flow through the other portions of the circulation. 

Line 18. Calculation of capillary pressure (PC) by adding the pressure gradient from the 
capillaries to the veins (PVG) to the pressure that is in the veins (PVS). 

Line 19. Calculation of the net pressure difference across the capillary membrane to cause 
movement of fluid molecules through the capillary pores (PCD) by adding capillary 
pressure (PC) and tissue colloid osmotic pressure (PTC) and subtracting plasma 
colloid osmotic pressure (PPC) and interstitial fluid pressure (PIF). 

Line 20. Calculation of rate of fluid movement through the capillary membrane (VTC) by 

multiplying pressure gradient across the capillary membrane (PCD) by the capillary 
filtration coefficient (CFC). 
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Line 21 . Calculofion of the driving pressure for moving fluid'into the lymphatics (PLD) 
by adding free interstitial fluid pressure (PIF), subtracting total tissue pressure 
(PTT), and adding a constant factor to account for lymphatic pumping. The 
total tissue pressure is considered to oppose lymph flow because of compression 
of the lymphatics while the interstitial fluid pressure is considered to promote 
lymph flow. 

Lines 22-23. Calculation of rate of lymph flow (VTL) from driving pressure for lymphatic 
flow (PLD). 

Line 24. Calculation of rate of change of fluid in interstitial fluid compartment (VT>D) 
by adding rate of movement of fluid into interstitial spaces from capillaries 
(VTC), subtracting rate of loss of fluid from the interstitial fluid compartment 
by way of lymph flow (VTL), and subtracting rate of movement of fluid from 
interstitial fluid compartment into cells (VID). 

Line 25. Calculation of total fluid in interstitial compartment (VTS) by integration of 
rate of change of fluid In interstitial compartment (VTD). 

Line 26. Calculation of rate of change of plasma volume (VPD) by adding fluid intake 
by drinking (TVD) and fluid return to circulation by way of lymphatics (VTL), 
and subtracting rate of movement of fluid through the capillaries (VTC), urinary 
output (VUD), and rate of fluid loss from the plasma through the pulmonary 
capillary membranes into the pulmonary spaces (DFP). 
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SUBROUTINE PULMON 

Program Listing: See Program 11. 

Flowchart: See Figure 18. 

Line 7. Calculation of plasma volume (VP) by integration of rate of change of plasma 
volume (VPD) (calculated on line 26 of SUBROUTINE CAMBD). 

Line 9. Calculation of pulmonary capillary pressure (PCP) from effects of pulmonary 
arterial pressure (PPA) and left atrial pressure (PLA). 

Line 10. Calculation of pressure of the free fluid in the pulmonary interstitial fluid 
spaces (PP1) from volume of free fluid in interstitial spaces of lungs (VPF). 

Line 11 . Calculation of protein in free fluid of lungs (CPN) by division of quantity of 
protein in pulmonary interstitial spaces (PPR) by volume of free fluid in lungs 
(VPF). 

Line 12. Calculation of Colloid osmotic £tessu fd 'bf, pro tel n in- free -fluitfi-effhe lungs (POS) 
from the concentration of protein in free fluid of lungs (CPN). 

Line 13. Calculation of pulmonary lymph flow (PLF) from driving pressure for pulmonary 
lymph flow. This driving pressure is determined gy adding pulmonary free 
interstitial fluid pressure (PPI) to a constant. 

Line 14. Calculation of rate of removal of protein from interstitial spaces of lungs by 
way of lymph (PPO) from pulmonary lymph flow (PLF) times concentration of 
protein in free fluid of lungs (CPN). 

Line 15. Calculation of rate of movement of protein- through pulmonary capillary membranes, 
into interstitial spaces (PPN) from the protein difference across the pulmonary 
capillary membrane. This latter quantity is determined by subtracting the 
pulmonary interstitial fluid concentration of protein (CPN) from the plasma 
concentration of protein (CPP). 

Line 16. Calculation of net rate of change of protein quantity in pulmonary interstitial 
fluids (PPD) by subtracting protein removal from interstitial spaces by lymph 
flow (PPO) from protein movement info interstitial spaces through capillary 
membranes (PPN). 

Line 17. This line tests to see if the total protein in the pulmonary fluids (PPR) as determined 
by integration of PPD calculated by line 16 is less than 0.025. If so, PPD is 
calculated directly from PPR. 
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Line 18. 

Line 19. 

Line 20. 
Line 21 . 

Line 22. 


Calculation of rate of flufd movement through pulmonary capillary membrane (PFI) 
by multiplication of pulmonary capillary filtration coefficient (CPF) times net 
pressure difference acorss pulmonary capillary membrances. The net pressure 
difference is obtained by adding pulmonary capillary pressure (PCP), and pulmonary 
interstitial fluid colloid osmotic pressure (POS), and subtracting pulmonary 
interstitial fluid pressure (PPI) and plasma colloid osmotic pressure (PPC). 

Calculation of net rate of change of fluid in pulmonary interstitial spaces (DFP) 
by subtracting rate of fluid movement from pulmonary interstitial spaces into 
pulmonary lymph and thence into the plasma (PLF) from fluid movement from 
pulmonary capillaries into interstitial spaces (PFI). 

This line tests to see if the volume of free fluid in the interstitial spaces of 
lungs (VPF) is less than 0.001. If so, DFP is calculated directly from VPF. 

Calculation of actual volume of free fluid in interstitial spaces of lungs (VPF) by 
integration of DFP. (Note that there is no calculation in the pulmonary fluid 
system for interstitial gel.) 

Calculation of total protein in the pulmonary fluids (PPR) by integration of PPD. 
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SUBROUTINE MISC2 


Program Listing: See Program 12. 

Flowchart: See Figure 19. 

Line 9. Calculation of hypertrophy of the left ventricle (HP L) by means of a time delay 

mechanism from the drive factor based on the arterial pressure (PA) and the strength 
of the left ventricle (HSL). 

Line 10. Calculation of hypertrophy of the right ventricfe (HPR) by means of a time delay 
mechanism from the drive factor based on the pulmonary arterial pressure (PPA) 
and the right heart strength (HSR). 

Lines 16-18. Calculation of the effect of tissue perfusion (expressed in terms of tissue 
oxygenation (POT)) on the mechanism for salt and water intake (STH). 



SUBROUTINE PROTEN 
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Progra m LisKng: See Program 13. 

Flowchart: See Figure 20. 

Line 11 . Calculation of rate of return of protein from interstitial spaces to the plasma 
(DPL) by multiplying rate of lymph flow (VTL) times concentration of protein 
in free interstitial fluid (CPI). 

Line 12. This line tests to see if the .capillary pressure (PC) is negative. If it is, PC is 
set to zero. 

Line 13. Calculation of the rate of protein movement through the capillary membrane 

(DPC) by multiplying the permeability of the capillaries to protein (considering 
that this permeability increases with the cube of capillary pressure (PC) and 
that its degree is set by a constant (CPD)) times the concentration difference 
between protein in plasma (CPP) and protein In the interstitial fluid (CPI). 

Line 14. Calculation of part of the rate of change of protein in the^free fluid'of the 

interstitial spaces (DPI) by subtracting the rate of return of protein to the plasma 
by way of the lymph (DPL) from the rate of movement of protein into interstitial 
spaces through the capillary membrane (DPC). (Note that the rate of movement 
of protein into the interstitial gel is not subtracted until line 25.) 

Lines 15-16. Calculation of undamped rate at which the liver produces plasma proteins 

(DLZ) from the difference between a reference factor (CPR) and the concentration 
of plasma proteins (CPP). 

Line 17. Calculation of rate at which the liver produces plasma protefa (DLP) by damping 
DLZ. 

Line 18. Calculation of the quantity of plasma protein (PRP) by integration of the rate 
of change of plasma protein as determined by adding the rate of formation 
of plasma protein by the liver (DLP) and the rate of return of proteins to the 
plasma by the lymphatics (DPL) and subtracting the rate of destruction or loss 
of plasma proteins by the body (DPO), the rate of loss of proteins through the 
capillary membrane (DPC), and the rate of loss of plasma protein through the 
pulmonary capillaries (PPD). 

Line 22. Calculation of the activity factor for protein in the interstitial fluid (PGX) 
by summing the effect of concentration of the protein in the gel (CPG) and 
the effect of concentration of hyaluronic acid in gel (CHY) to exacerbate the 
colloid osmotic pressure effect of protein in the gel. 
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Line 23. 

Line 24. 
Line 25. 


Calculation of the rate of protein movement into gel (GPD) by multiplying the 
activity difference between the free fluid and gel times a constant. This 
activity difference is itself calculated by multiplying the gel volume (VG) times 
the protein difference between interstitial fluid (CPI) and gel protein activity 
(PGXj. 

Calculation of the quantity of protein in the gel (GPR) by integration of the 
rate of movement of protein into gel (GPD). 

Calculation of the quantity of protein in the free rnterstifial fluid (1FP) by 
integration of the rate of increase of protein in the gel. This latter quantity 
is obtained by subtracting the rate of movement of protein from free interstitial 
fluid into interstitial gel fluid (GPD) from DPI calculated on line 14. 
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SUBROUTINE KIDNEY 


Program Listing: See Program 14. 
Flowchart See Figure 21. 


Lines 8-10. Calculation of degree of auto regulatory feedback at macular densa (GF3) from 
glomerularfiitrafion rate (GFN). This in turn, controls afferent arteriolar 
resistance. The factor GF4 controls the feedback gain of the autoregulatory 
loop. 


Line 1 1 . 


Line 12. 

Line 13. 
Line 14. 
Line 15. 
Line 16. 
Line 17. 
Line 18. 

Line 19. 


Calculation of the afferent arteriolar resistance to the midpoint of the glomeruli 
(AAR) by multiplying the autonomic effect by the viscosity of the blood (VIM) 
and by the degree of autoregulatory feedback at the macular densa (GF3). 

The autonomic effect is calculated from the autonomic multiplier (AUM) and 
the factor ARF which increases or decreases the effect of the autonomies on 
the kidneys. A value of ARF of zero will set the sensitivity to zero. 

Calculation of the renal resistance (RR) by addition of the afferent arteriolar 
resistance (AAR) and the efferent (postglomerular) resistance. The efferent 
resistance is calculated by multiplying a constant times the viscosity of the 
blood (VIM). 

Calculation of renal arterial pressure (PAR) by subtracting the Goldblatt 
parameter (GBL) from the arterial pressure (PA). 

Calculation of the blood flow through the kidneys (RFN) (assuming the kidneys 
are normal) by dividing the renal arterial pressure (PAR) by the renal resistance (RR). 

Calculation of renal blood flow (RBF) by multiplying the blood flow through the 
normal kidney (RFN) by the degree of normality of the kidneys (REK). 

Calculation of the pressure drop in the afferent arterioles (APD) by multiplying 
the normal renal blood flow (RFN) by the afferent arteriolar resistance (AAR). 

Calculation of glomerular pressure (GLP) by subtracting the pressure drop in 
the afferent arterioles (APD) from the renal arterial pressure (PAR). 

Calculation of glomerular filtration pressure (PPL) by subtracting plasma colloid 
osmotic pressure (PPC) and a constant value representing Bowman's capsule 
pressure from the glomerular pressure (GLP). 

This saves the value of GFN as GF1 . 
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Line 20. 

Line 21 , 

Line 22. 
Line 23. 

Lines 24- 
Line 30. 

Line 31 . 
Line 32. 

Line 33. 


Calculation of glomerular filtration (if the kidneys are normal) (GFN)'by 
multiplying glomerular filtration pressure (PFL) times a constant representing 
the glomerular filtration coefficient. The factors GF2 and Z represent damping 
effects. 

This is a test to see if the normal glomerular filtration has changed by more 
than 0.002. If it has, the calculation goes back to line 8 until stabilization 
is obtained. 

Calculation of actual glomerular filtration rate (GFR) by multiplying normal 
filtration rate (GPN) by degree of normality of the kidneys (REK). 

Calculation of total tubular reabsorption (TRR) by adding the amount of glomerular 
filtrate that is reabsorbed irrespective of control by aldosterone and antidiuretic' 
hormone (approximately 0.8 of GFR) and the maximum amount of fluid capable 
of being reabsorbed by the tubules each minute under the control of aldosterone 
and antidiuretic hormone, and by subtracting the amount of fluid not reabsorbed 
but could have been reabsorbed under the control of aldosterone and antidiuretic 
hormone. 

■25. Calculation of the rate of urinary output (VUD) by subtracting total tubular 
reabsorption (TRR) from glomerular filtration rate (GFR). 

Calculation of (undamped) rate of sodium loss in urine (NOZ) assuming a normal 
concentration of sodium in the urine of 100 meq/liter and assuming that there 
are three factors that affect this output; the volume of urine formed each minute 
(VUD), the aldosterone multiplier effect (AM). and the "third factor" effect 
related to the change in concentration of sodium in the extracellular fluid (CNE). 

Calculation of the rate of sodium loss in the urine (NOD) by damping NOZ. 

Calculation of the net rate of change of sodium in the extracellular fluid (NED) 
from intake of salt, expressed as basic intake of sodium (NID) times appetite 
factor (STH), anS sodium loss (NOD). 

Calculation of quantity of sodium in extracellular fluid (NAE) by integration of 
net rate of change of sodium in extracellular fluids (NED). 
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SUBROUTINE IONS 


Program Listing; See Program 15. 


Flowchart; See Figure 22. 


Line 7. 


Line 8. 


Line 9. 


Line 10. 


Line 1 1 . 


Line 12. 


Line 13. 


Line 14. 


Line 15. 


Line 16. 


Calculation of extracellular fluid volume (VEC) by addition of plasma volume 
(VP), volume of fluid in the interstitial spaces of the systemic circulatory bed 
(VTS), and volume of fluid in the interstitial spaces of the lungs (VPF). 

Calculation of concentration of potassium in extracellular fluid (CKE) by division 
of the quantity of potassium in the extracellular fluid (KE) by volume of 
extracellular fluid (VEC). 

Calculation of the rate of renal excretion of potassium (KOD) by multiplying 
the degree of normality of the kidneys (REK) by the sum of the non-aldosterone 
controlled portion of potassium excretion and the aldosterone (AM) controlled 
portion of potassium excretion. 

Calculation of total expected quantity of potassium in the intracellular fluid 
under equilibrium conditions (KIR) by addition of a constant value representing 
potassium in cells that is not dependent upon extracellular potassium concentration 
and the quantity of potassium inside the intracellular fluid that is dependent 
upon extracellular potassium concentration (CKE). 

Calculation of potassium gradient that causes potassium movement into the 
cells (KIE) by subtracting the actual level of potassium in the cells (KIR). 

Calculation of rate of movement of potassium through cell membranes (KCD) 
by multiplying difference between expected and actual potassium levels 
(KIE) times a constant for potassium diffusion. 

Calculation of quantity of potassium in the intracellular fluid (Kl) by integration 
of the rate of movement of potassium into the intracellular fluid (KCD) 

Calculation of net rate of change of potassium in the interstitial fluid (KED) 
by subtracting the rate of loss of potassium in the urine (KOD) and rate of 
movement of potassium into the cells (KCD) from the rate of potassium intake (KID). 

Calculation of total quantity of potassium in extracellular fluid (KE) by 
integration of the net rate of change of potassium in extracellular fluid (KED). 

Calculation of concentrations of potassium in intracellular fluids (CKI) by 
division of quantity of potassium in intracellular fluids (Kl) by volume of 
intracellular fluid (VIC). 
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Line 17. 
Line 18. 

Line 19. 
Line 20. 


Calculation of concentration of sodium in extracellular fluid (CNA) by division 
of quantity of sodium in extracellular fluid (NAE) by volume of extracellular 
fluid (VEC). 

Calculation of concentration gradient between intracellular and extracellular 
fluids by subtracting the concentration of sodium in the extracellular fluids 
(CNA) as an indicator of the osmolarity of the extracellular fluid from the 
concentration of potassium in the intracellular fluids (CKI) as an indicator 
of the osmolarity inside the cells. 

Calculation of the rate of movement of water into cells from the extracellular 
fluid space (VID) from the osmolarity factor difference (CCD). 

Calculation of volume of water in cells (VIC) by integration of the rate of 
movement of water into the cells (VID). 



29 


SUBROUTINE GELFLD 

Program Listing: See Program 16. 

Flowchart: See Figure 23. 

Line 7. Calculation of concentration of hyaluronic acid'In gel of interstitial spaces (CHY) 
by dividing quantity of hyaluronic acid (HYL) by volume of gel (VG). 

Line 8. Calculation of elastic suction of the hyaluronic acid in the tissues caused 

by elastic recoil of the gel (PRM) from the concentration of hyaluronic acid 
in gel (CHY). 

Line 9. Calculation of colloid osmotic pressure of the ge! reticulum caused by Donnen 
equilibrium (PGR) of hyaluronic acid (CHY). 

Line 10. Calculation of the concentration of protein in gel (CPG) by division of quantity 
of protein in gel (GPR) by volume of gel (VG). 

Line 1 1. Calculation of colloid osmotic pressure of the firotein in the gel (PGP) by 
multiplying the activity of the protein in the gel (PGX) by a constant. 

Line 12. Calculation of total colloid osmotic pressure of the fluid inside the gel (PGC) 

by adding that caused by the reticulum itself (PGR) to that caused by the protein 
in the gel (PGP). 

Line 13. Calculation of the volume of free fluid in the interstitial spaces (VIF) by 

subtracting volume of gel fluid (VG) from total interstitial fluid volume (VTS). 

Line 14. Calculation of solid tissue pressure (PTS) by graphical means (See Figure 17.) 
from volume of free interstitial fluid (VIF). 

Line 15. Calculation of pressure of free interstitial fluid (PIF) by subtracting solid 
tissue pressure (PTS) from total tissue pressure (PTT). 

Line 16. Calculation of concentration of protein in free interstitial fluid (CPI) by 
dividing quantity of protein (IFP) by volume of free fluid (VIF). 

Line 17. Calculation of colloid osmotic pressure of free interstitial fluid (PTC) by 

multiplying concentration of protein in the' interstitial fluid (CPI) times a constant. 

Line 18. Calculation of net mechanical forces attempting to cause movement into or 

out of gel (PGH) by summing the elastic recoil suction of gel (PRM), solid tissue 
pressure (PTS), and interstitial fluid pressure (PIF). 
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Line 19. Calculation of the rate of movement of fluid between gel and free interstitial 

fluid (VGD) by multiplying the resistance factor (V2D) by net pressure difference 
at the gel surface. This net pressure difference is obtained by subtracting the 
colloid osmotic pressure of the free fluid of the interstitial spaces (PTC) and 
the mechanical suction of the gel (PGH) from the sum of the total colloid 
osmotic pressure of the gel fluid (PGC) and the pressure of the interstitial fluid 
(PIF). 

Lines 20-21 . Calculation of the gel volume (VG) by integration of the net movement of 
fluid through the gel surface (VGD). 

Line 22. This is a test to see if the net movement of fluid through the gel surface (VGD) 
exceeds 0.0(12. If so the program returns to line 7 until VGD<0.012. 
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III. Typical Model Experiments 

In this chapter, the model is used to simulate a few typical experimental situations 
in order to illustrate the models general utility. 

Experiment 1: Hypertension in a salt loaded, renal deficient patient (Table 1). 

Variables monitored- extracellular fluid volume (VEC), blood-volume (VB), sympathetic 

stimulation (AU), cardiac output (QLO), total peripheral -resistance 
(RTP), aortic pressure (PA), heart rate (HR), angiotensin concentra- 
tion (normal=l)(ANC), urinary output (VUD). 

Changes made-After 2 hours, the renal mass was reduced to 0.3 normal v(REK=0.3). 

After 4 days, the salt intake was increased to five times 
normal (N!D=0.5). Total experimental time was 8 days. 

Observations-The initial decrease in renal mass had only a slight effect on variables 

monitored with the exception of a slight decrease in cardiac 
output and simultaneous increase in total peripheral resistance. 

The arterial pressure elevated a small amount. Increase of salt 
load caused more dramatic effects. The extracellular volume 
and blood volume rose, the cardiac output increased considerably 
and then stabilized, while the total peripheral resistance fell. 

The rise in cardiac output increased the arterial pressure. After 
120 hours, the cardiac output stabilized, while the peripheral 
resistance rose. The arterial pressure continued to increase, which 
demonstrates that the increase in total peripheral resistance, 
not cardiac output, was responsible for the long-term hypertension. 
Note that urinary output increased during salt loading because 
of the effect of high salt intake on thirst. 

Experiment 2: Nephrosis due to protein loss by plasma (Table 2). 

Variables monitored-urinary output (VUD), interstitial fluid gel volume (VG), 

total interstitial fluid volume (VTS), plasma volume (VP), total 
plasma protein (PRP), interstitial fluid pressure (PIF), aortic 
pressure (PA), cardiac output (QLO). 

Changes made-After 1 hour, the- rate of loss of -plasma-protein -was 'increased seven- 
fold (DPO=0.05>. After 108 hours the rate of loss of plasma protein 
was put back at three times normal (DPO^). 021). Total 
experimental time was 5 ]/2 days. 
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Observations-The initial decrease in plasma protein initiated slight decreases in both 

arterial pressure and cardiac output and marked decrease in 
in urinary output. The fluid thus retained caused swelling of 
the interstitial gel. The volume of free interstitial fluid (VTS-VG) 
remained relatively stable until the interstitial fluid pressure 
rose into the positive range. Then, marked edema occured with 
sharp drop of cardiac output. When the rate of renal loss of 
protein was increased to the point where the liver could 
increase the plasma protein level, the edema was relieved with 
high diuresis and increased cardiac output. 

Experiment 3: Severe muscle exercise (Table 3). » 

Variables monitored-urinary output (VUD), muscle venous oxygen pressure (PVO), 

muscle cell oxygen pressure (PMO), aortic pressure (PA)/ sympathetic 
stimulation (AUP), cardiac output (QLO), muscle blood flow (BFM), 
rate of oxygen utilization by muscle cells (MMO). 

Changes made-After 30 seconds, the exercise parameter was changed to 60 times its 

normal value (EXO=60), corresponding to a whole body metabolism 
increase of approximately 15 times. At the same time, the time 
constant for local vascular response to metabolic activity was 
reduced by 1/40 (A4K=0.025), the damping 'factor Z was increased 
5 fold (Z=5.), and the factors Z5, Z6, andZ8 were modified 
(Z5=l., Z6-=10., Z8=3.). The value of 13 was also set at zero 
to prevent long integration steps. After 2 minutes, the value of 
EXC was put back to normal (EXC= 1). After 5 minutes 13 was set 
back to normal (I3=£0). 

Observations-At the onset of exercise, cardiac output and muscle blood flbw increased 

considerably within seconds. Urinary output fell to obligatory 
level while arterial pressure rose moderately. Muscle cell and 
venous PO 2 fell rapidly. Muscle metabolic activity showed an 
instantaneous increase, but then decreased considerably because 
of the development of a metabolic deficit in the muscles. When 
exercise was stopped, muscle metabolic activity fell to below 
normal, but cardiac output, muscle blood flow, and arterial 
pressure remained elevated for a while as the person was repaying 
his oxygen debt. 

Experiment 4: Atrioventricular fistula (Table 4) 

Variables monitored-extracellular fluid volume (VEC), blood volume (VB), autonomic 

stimulation (AU), cardiac output (QLO), total peripheral 
resistance (RTP), aortic pressure (PA), heart rate (HR), angiotensin 
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concentration (relative to 1 at normal )(ANC)/ urinary output (VUD). 

Changes made-After 2 hours a fistula was created which would double-cardiac output 

(Fl$=0.05). After 4 days the fistula was closed (FIS=0.0). Total 
experimental time was 9 days. 

Observations-Opening the fistua! caused an immediate dramatic change in cardiac output, 

total peripheral resistance, and heart rate. Urinary output decreased 
to obligatory levels. As the body adapted, extracellular fluid 
volume and blood volume increased to compensate for the fistula 
with the result that after a few days arterial pressure, heart rate, 
and urinary output were near normal levels. Cardiac output 
doubled and peripheral resistance halved. When the fistula was 
closed, dramatic effects again occurred with rapid decrease in 
cardiac output, rapid increase in peripheral resistance, moderate 
increase in arterial pressure, and moderate decrease in heart rate. 
Marked diuresis reduced extracellular fluid volume and blood 
volume to normal or slightly' below. After several days, the patient 
was nearly normal. 

IV. Model Characteristics and Interfacing Problems 


The experiments contained in Chapter Hi, and others not reported here, demonstrate that 
the model of Guyton is capable of responding in a correct overall fashion to a variety of 
stress conditions despite the fact that the model itself is based on the gross function of the 
many different parts of circulation. Most of the subsections of the model are, in fact, 
developed at a crude level with minute details completely absent. The overall correctness 
of the model predictions is a result of the facts that the interactions between the basic 
regulatory mechanisms of the body possesses an inherent stability and that this stability is 
more important than the details of any one mechanism. 

The model presented and discussed here may be viewed as a controlled system plus 
controlling system with the controlling system having three major components: local control, 
hormonal control, and autonomic control. These controls act to drive the controlled 
system to the appropriate level in response to stress. There are no thermal regulatory components 
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present in either the controlled or controlling system. Respiratory elements are absent 
as well with the exception of the effect of pulmonary Interstitial fluid on aortic .oxygen 
saturation. Hydrogen ion levels are not considered. Only the major cations, Na + and K + , 
are treated. The model may be classified as an intermediate to long-term model with 
simulations of the order of days or weeks being the primary concern, although short-term 
simulations, as in the exercise experiment of Chapter III, are possible. 

It should be possible to modify the present model to accomodate stimuli to which the 
model does not presently respond. Probably the most important of these stimuli are concerned 
with temperature regulation in the body and the regulation of respiration. Related to the 
latter is the problem of hydrogen ion regulation. Large mathematical models of the thermal 
regulatory (5) and respiratory regulatory systems (6) are presently available and have Been 
subjected to considerable study. One basic question that immediately arises concerns the 
possibility of interfacing or combining these models with the circulatory model discussed 
here in order to include simulated responses by a patient to a wider variety of stimuli. 

Usually, it is not possible to directly interface different models. This is so because 
different models generally utilize distinct approaches to the study of their respective systems 
with the different types of both controlled and controlling equations resulting in partial or 
complete incompatability of models. Some models are developed as short-term models only, 
and their use in conjunction with an intermediate or long-term model would make little 
sense. Often the major controlling feature of one model is completely absent from a model 
of a different system, in spite of the fact that the second system plays some role in the 
regulation of this component. At times, two different controlling systems drive the same 
component. Here, a decision as to the proper way of combining these driving forces must 
be made and this decision may be difficult to arrive at. 
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There are several distinct approaches which may be utilized to interface subsystem 
models and so form an overall composite model. To begin with, the individual models 
could simply be run simultaneously under a single monitor system with no thought of making 
the models interact. Thus, each model would react independently and output would 
be selected only from the model of interest at any particular time. This solution is not very 
satisfactory, both from the point of view of resource availability and the point of view of 
physiological realism. Such a system would not really represent an overall regulatory 
model and would answer few, if any, questions that the individual model systems alone 
could not answer. A second approach to forming a composite model would be to identify 
all elements in common from each of the controlling elements of the respective subsystem 
models and to implement a new overall controlling system which receives dynamic data from 
each subsystem and then regulates each controlled system according to overall current 
information. This approach is much more attractive than the first approach mentioned, but 
suffers at least one serious disadvantage, other than the basic one of how to write the 
controller equations. This disadvantage stems from the fact that the subsystem models 
themselves are designed to be realistic for different time periods. Thus, the respiratory 
model may be reliable for experiment times of 20-30 minutes while the circulatory model 
of Guyton may be used for experimental times of days or weeks. This limitation would make 
a composite model of this type almost prohibitively expensive to run for any reasonable 
length of time because of the simple fact that the time limiting step size must be controlled 
by the model with the shortest response. 

A third path to the formation of a composite model appears to offer the most attractive 
alternative. This approach is similar, in some ways, to the second alternative mentioned 
above. Instead of implementing the detailed subsystem models, this third alternative would 
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utilize the Guyton model discussed in this report as a basic master model with the other 
subsystem models included in their gross function only. This new master model would be 

carefully planned so as to be compatible with the detailed subsystem models. The detailed 

/ 

subsystem models would be utilized only in the event that detailed response of a particular 
subsystem is of interest. Otherwise, only crude responses of the gross system would be 
calculated. Thus, the overall model would be capable of producing long-term regulatory 
features while the detailed subsystem models would be capable of examining short-term 
transient effects. This alternative presents a considerable challenge to modellers In the 
form of compatibility requirements, but it must be remembered that the overall system being 
modelled is capable of functioning as one unit. It is suggested that this alternative be 
explored in depth in future researches. 
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Table 1. Salt Loaded Renal Deficient Kidney 


HOUR 

VEC 

VB 

AU 

QLO 

RTP 

PA 

HR 

ANC 

VUD 

0 

15.0 

5.00 

0.989 

5.12 

19.4 

99.7 

71.7 

0.996 

1.03 

2 

15.0 

REK 

5.00 
= 0.3 

0.991 

5.11 

19.4 

99.4 

71.8 

1.002 

0.98 

3 

15.0 

5.03 

0.938 

4.87 

21.9 

103.6 

69.4 

0.286 

0.55 

6 

15.0 

5.06 

0.815 

4.76 

21.9 

105.1 

64.7 

0.209 

1.02 

12 

15.0 

5.06 

0.840 

4.74 

22.0 

104.0 

65.6 

0.224 

0.92 

18 

15.0 

5.04 

0.871 

4.69 

22.2 ’ 

104.0 

66.8 

0.238 

0.83 

24 

15.0 

5.04 

0.891 

4.67 

22.3 

104.1 

67.5 

0.244 

0.78 

30 

15.0 

5.04 

0.899 

4.67 

22.4 

104.5 

67.8 

0.245 

0.77 

36 

15.0 

5.04 

0.901 

4.68 

22.4 

104.9 

67.9 

0.244 

0.78 

42 

15.0 

5.04 

0.902 

4.69 

22.5 

105.3 

68.0 

0.241 

0.79 

48 

15.0 

5.04 

0.904 

4.71 

22.5 

105.6 

68.0 

0.240 

0.80 

54 

15.0 

5.04 

0.905 

4.72 

22.5 

106.1 

68.1 

0.239 

0.81 

60 

15.0 

5.04 

0.908 

4.71 

22.6 

106.3 

68.2 

0.238 

0.82 

72 

15.0 

5.05 

0.914 

4.71 

22.7 

106.9 

68.5 

0.237 

0.82 

96 

15.0 

N1D 

5.05 
= 0.5 

0.926 

4.72 

22.8 

107.9 

68.9 

0.236 

0.82 

97 

15.1 

5.07 

0.910 

4.96 

22.7 

110.5 

68.5 

0.221 

0.90. 

102 

15.6 

5.23 

0.646 

’ 5.44 

21.2 

117.0 

58.8 

0.012 

2.38 

108 

15.7 

5.31 

0.579 

5.73 

20.6 

119.3 

56.7 

0.010 

3.08 

114 

15.7 

5.34 

0.571 

5.84 

20.6 

121.1 

56.5 

0.010 

3.25 

120 

15.6 

5.35 

0.580 

5.88 

, 20.8 

122.5 

. 56.9 

0.010 

3.23 

126 

15.5 

5.35 

0.595 

5.87 

21.0 

123.8 

57.5 

0.010 

3.16 

138 

15.5 

5.35 

0.627 

5.86 

21.4 

126.3 

58.7 

0.010 

3.01 

150 

15.4 

5.35 

0.656 

5.87 

21.8 

128.6 

59.8 

0.010 

2.93 

162 

15.4 

5.35 

0.682 

5.87 

22.2 

130.8 

' 60.8 

0.010 

2.86 

174 

15.5 

5.35 

0.706 

5.87 

22.5 

132.7 

61.7 

0.010 

2.81 

192 

15.5 

5.35 

0.739 

5.87 

22.9 

135.4 

63.0 

0.010 

2.73 
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iOUR 


1 

2 

4 

6 

9 

13 

17 

21 

25 

31 

37 

43 

49 

55 

61 

67 

73 

79 

85 

91 

94 

95 
97 

103 

108 

109 

115 

121 

127 

132 


Table 2. Nephrosis Due to Plasma Protein Loss 


VUD VG VTS VP PRP PIF PA QLO 


1.03 

11.5 

12.1 

2.96 

DPO = 

■ 0.05 



0.94 

11.5 

12.1 

2.95 

0.36 

11.6 

12.2 

2.92 

0.20 

11.7 

12.3 

2.88 

0.20 

11.9 

12.5 

2.83 

0.20 

12.1 

12.7 

2.76 

0.20 

12.3 

12.9 

2.73 

0.20 

12.4 

13.1 

2.74 

0.20 

12.6 

13.2 

2.74 

0.20 

12.7 

13.4 

2.79 

0.98 

12.7 

13.4 

2.84 

0.84 

12.8 

13.5 

2.84 

0.75 

12.9 

13.6 

2.83 

0.69 

13.0 

13.7 

2.83 

0.69 

13.1 

13.9 

2.82 

0.56 

13.2 

14.0 

2.82 

0.41 

13.4 

14.2 

2.80 

0.20 

13.6 

14.5 

2.77 

0.20 

13.8 

14.8 

2.73 

0.20 

14.0 

15.0 

2.71 

0.20 

14.0 

15.2 

2.71 

0.20 

14.1 

15.2 

2.71 

0.20 

14.2 

15.3 

2.71 

0.20 

14.3 

15.6 

2.70 

0.20 

14.5 

15.9 

2.68 

DPO 

=0.021 



0.20 

14.5 

15.9 

2.68 

0.91 

14.6 

16.0 

2.75 

4.47 

14.2 

15.4 

2.92 

2.29 

13.9 

14.9 

2.90 

1.71 

13.8 

14.8 

2.88 


207.7 

-5.98 

99.7 

5.09 

205.4 

-5.93 

98.8 

5.04 

200.8 

-5.78 

96.4 

4.98 

196.7 

-5.60 

94.8 

4.87 

191.2 

-5.32 

92.1 

4.75 

184.8 

-4.89 

89.0 

4.66 

179.4 

-4.54 

90.6 

4.59 

174.8 

-4.24 

90.6 

4.56 

170.9 

-3.99 

90.7 

4.58 

165.9 

-3.74 

94.0 

4.63 

161.6 

-3.66 

94.3 

4.66 

157.2 

-3.43 

93.2 

4.66 

-152.7 

-3.15 

92.6 

4.66 

147.9 

-2.83 

91.8 

4.63 

142.9 

-2.50 

91.1 

4.61 

137.8 

-2.14 

90.0 

4.58 

132.3 

-1.76 

88.7 

4.53 

126.5 

-1.28 

86.2 

4.45 

120.4 

-0.78 

84.6 

4.38 

114.2 

-0.31 

84.2 

4.29 

111.2 

-0.08 

83.8 

4.28 

110.1 

0.00 

83.7 

4.27 

108.1 

0.15 

83.4 

4.25 

102.0 

0.55 

82.3 

4.19 

96.9 

0.91 

81.7 

4.12 

97.5 

0.96 

82.1 

4.18 

101.1 

1.08 

85.8 

4.42 

105.3 

-0.01 

89.9 

4.70 

108.6 

-0.44 

85.1 

4.67 

110.6 

-0.51 

85.4 

4.66 
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SECS 

VUD 

Toble 3. 
PVO 

Simulation of Severe Muscle Exercise 
PMO PA AUP QLO 

BFM 

MMO 

30 

1.07 

39.9 

8.00 

99.9 

0.98 

5.12 

0.99 

59.9 


EXC 

= 60. 







31 

1.00 

37.1 

5.25 

100.6 

1.14 

5.29 

0.93 

352.8 

32 

0.76 

29.1 

3.96 

103.6 

1.38 

5.60' 

1.09 

331.6 

33 

0.32 

25.7 

3.13 

108.5 

1.69 

6.26 

1.97 

307.5 

34 

0.20 

26.8 

2.70 

112.1 

1.94 

7.13 

3.13 

291.4 

45 

0.20 

28.1 

1.35 

137.5 

3.49 

14.0 

11.1 

210.2 

60 

0.20 

27.9 

1.08 

140.8 

4.10 

17.7 

14.6 

186.0 

75 

0.20 

27.9 

1.00 

142.2 

4.31 

18.9 

15.8 

178.1 

90 

0.20 

27.9 

0.97 

142.4 

4.40 

19.4 

16.4 

175.1 

105 

0.20 

27.9 

0.96 

142.4 

4.45 

19.7 

16.6 

174.0 

120 

0.20 

27.9 

0.95 

142.4 

4.49 

19.8 

16.7 

173.5 


EXC 

= 1.0 







121 

0.20 

27.9 

0.98 

142.1 

4.41 

19.7 

16.4 

29.5 

130 

0.20 

28.0 

1.32 

125.6 

3.37 

14.7 

11.8 

34.0 

140 

0.20 

28.4 

1.78 

118.8 

2.67 

11,3 

7.89 

39.7 

150 

0.20 

29.1 

2.29 

115.6 

2.19 

9.. 06 

5.39 

45.0 

165 

0.20 

30.3 

3.08 

111,1 

1.76 

7.31 

3.40 

51.3 

180 

0.20 

31.6 

3.85 

108.0 

1.52 

6.45 

2.42 

55.6 

195 

0.20 

33.0 

4.54 

106.0 

1.38 

5.99 

1.90 

58.3 

210 

0.20 

34.2 

5.16 

104.6 

1.29 

5.71 

1.60 

59.8 

225 

0.20 

35.2 

5.68 

103.5 

1.23 

5.52 

1.41 

60.6 

240 

0.24 

36.0 

6.12 

102.6 

1.18 

5.40 

1.29 

60.9 

255 

0.32 

36.7 

6.48 

101.9 

1.15 

5.31 

1.20 

60.9 

270 

0.38 

37.3 

6.78 

101.2 

1.13 

. 5.25 

1.14 

60.9 

285 

0.42 

37.7 

7.01 

100.8 

1.11 

5.20 

1.10 

60.8 

300 

0.46 

38.0 

7.20 

100.4 

1.09 

5.17 

1.06 

60.8 

360 

0.47 

38.4 

7.58 

100.4 

1.09 

5.17 

1.04 

60.8 

420 

0.49 

38.6 

7.72 

100.6 

1.09 

5.18 

1.02 

60.8 

480 

0.50 

38.8 

7.79 

100.7 

1.09 

5.18 

1.01 

60.8 

540 

0.52 

38.9 

7.82 

100.8 

1.08 

5.18 

1.00 

60.8 
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Table 4. Simulation of Atrioventricular Fistula 


iOUR 

VEC 

VB 

AU 

QLO 

RTP 

PA 

HR 

ANC 

VUD 

2 

15.0 

FIS 

5.00 
= 0.05 

0.99 

5.11 

19.4 

99.4 

71.8 

1.00 

0.98 

3 

15.1 

5.00 

1.67 

8.22 

9.83 

81.9 

100.4 

1.69 

0.20 

6 

15.5 

4.94 

1.97 

8.74 

10.2 

. 88.4 

112.2 

3.12 

0.20 

12 

16.2 

5.14 

1.51 

9.16 

9.82 

92.0 

95.0 

2.72 

0.20 

18 

16.6 

5.28 

1.22 

9.48 

9.53' 

92.7 

-84.0 

2.74 

0.20 

24 

16.7 

5.44 

1.01 

9.72 

9.45 

95.0 

77.2 

2.55 

0.26 

30 

16.9 

5.55 

0.96 

9.87 

9.33 

95.8 

76.6 

2.35 

0.71 

36 

16.9 

5.64 

0.93 

9.99 

9.23 

96.4 

76.5 

2.14 

1.01 . 

42 

17.0 

5.70 

0.92 

10.10 

9.15 

96.8 

76.6 

1.94 

1.17 

48 

16.9 

5.75 

0.91 

10.18 

9.09 

97.0 

76.9 

1.77 

1.24 

54 

16.9 

5.78 

0.92 

10.25 

9.04 

97.2 

77.3 

1.62 

1.25 

60 

16.9 

5.80 

0.92 

10.31 

9.00 

97.3 

.77.7 

1.50' 

1.24 

66 

16.8 

5.82 

0.93 

10.35 

8.97 

97.4 

78.0 

1.39 

1.22 

72 

16.8 

5.82 

0.93 

10.39 

8.95 

97.6 

78.3 

1.31 

1.19 

78 

16.8 

5.83 

0.94 

10.42 

8.94 

97.6 

78.6 

1.25 

1.16 

84 

16.7 

5.83 

0.94 

10.44 

8.93 

97.7 

78.8 

1.19 

1.13 

90 

16.7 

5.83 

0.95 

10.46 

8.92 

97.8 

79.0 

1.15 

1.10 

96 

16.7 

FIS 

5.83 
= 0.0 

0.96 

10.47 

8.92 

97.8 

79.1 

1.12 

1.08 

97 

16.3 

5.75 

0.44 

7.70 

14.4 

111.6 

54.9 

0.04 

13.12 

99 

15.1 

5.54 

0.48 

6.96 

15.7 

108.5 

54.2 

0.01 

11.70 

102 

14.1 

5.28 

0.69 

5.92 

18.1 

106.4 

60.8 

0.65 

6.35 

108 

13.9 

4.93 

1.10 

4.87 

19.5 

94.4 

75.7 

1.54 

0.20 

114 

14.2 

4.77 

1.44 

4.58 

19.9 

89.7 

88.7 

2.07 

0.20 

120 

14.5 

4.73 

1.42 

4.43 

20.5 

90.3 

87.9 

2.47 

0.20 

132 

14.9 

4.75 

1.24 

4.46 

20.4 

90.9 

81.0 

2.90 

0.20 

144 

15.1 

4.92 

0.81 

4.78 

20.3 

98.0 

64.3 

2.35 

1.74 

156 

15.0 

4.95 

0.86 

4.88 

20.0 

97.5 

66.4 

1.91 

1.33 

168 

15.0 

4.96 

0.89 

4.93 

19.8 

97.7 

67.7 

1.61 

1.19 

180 

15.0 

4.98 

0.91 

4.99 

19.6 

98.0 

68.6 

1.41 

1.12 

192 

15.0 

4.98 

0.93 

5.02 

19.6 

98.2 

69.4 

1.27 

1.08 

216 

15.0 

4.99 

0.96 

5.05 

19.4 

98.5 

70.3 

1.11 

1.03 
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Figure 1. Flow Chart for Circulatory, Fluid, and Electrolyte Regulation Model. 
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Figure 2. Flow Chart for Subroutine HEMO. 
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Figure 3. Relationship Between Effective Arterial Pressure (RA2) 
and Left Ventricular Pumping Effectiveness (LVM), 
Function 1 . 
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Figure 4. Relationship Between Right Atrial Pressure (PRA) 

and Normal Output of Right Atrium (QRN), Function 2. 
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Figure 5. Relationship Between Effective Pulmonary Arterial 
Pressure (PP2) and Pumping Effectiveness of Right 
Ventricle (RVM), Function 3. 
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Figure 6. Relationship Between Left Atrial Pressure (PLA) 
and Normal Output of Left Ventricle (QLN), 
Function 4. 
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Figure 7. Flow Chart for Subroutine AUTO. 



NASA 

S- 73- 31085 





56 


Figure 8. Flow Chart for Subroutine HORMON. 
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Figure 9. Relationship Between Arterial Pressure (PA) 

and Aldosterone Secretion (AMP), Function 7. 
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Figure 10. Flow Chart' for Subroutine BLOOD. 
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Figure 11. Flow Cjjarf for Subroutine MUSCLE. 
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Figure 12. Flow Chart for Subroutine AUTORG. 
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Figure 13. Flow Chart for Subroutine ADH. 
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Figure 14. Flow Chari - for Subroutine MISC1 . 
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Figure 15. Flow Chari- for Subroutine HEART. 
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Figure 16. Flow Chart for Subroutine CAPMBD. 
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Figure 17. Relationship Between Volume of Free Interstitial Fluid 
(V I F) and Solid Tissue Pressure (PTS), Function 6. 
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Figure 18 . Flow Chart for Subroutine PULMON. 
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Figure 19. Flow Chart for Subroutine MISC2. 



r?4«o 










80 


Figure 20. Flow Chart for Subroutine PROTEN. 
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Figure 21. Flow Chart for Subroutine KIDNEY. 
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Figure 22. Flow Chart for Subroutine IONS. 
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Figure 23. Flow Chari for Subroutine GELFLD. 
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Program 1. 


SUBROUTINE HEMO. 



IV (VER lib ) 
subroutine 


Su'URC* LISTING} HEMC SUBROUTINE 05/11/73 PAGE 
HC-MO f AMMjrANMj ANy,ANY* ANZ^ARMjAUHIaUMIaUY* AVE/ Bf-M, BFN» 

r k.-j _ r M'a . r M7 _ * w n*c . fM a . no a rv» a M/p*:Tc.uMn.l4DI . 


C, 

HPRjH 

qroTs 

RVS/U 
VRC * V 


^S^fiW'5 v jOAS#DI.Aj|3PA*DRA*OVSiFis^HeiDibfLf 
>ySL,HS^l2 *LVM,PA ,PAM*P42*PC # PGL# POs,PLA, 
T P f. 1 / P P 2 , » R A , P RT j P V S /QA 0 * 0 L NTO'LTl 7 U? 0 rQ R F 7 g R 


DIMENSION FLJH1(14)^FUN2(14)/PUAI3(U)«PUN4(14) 
REAL 1 2 > L V M 


u I l»U WLn l Wr\ * ^IIN«nu,3 U I, u 1 

hemodynamics 

VBD*VP*YRC-WS-VAS-VLA»VPA"VRA 


VPA B VPA+DPA#I2+VBD*. 155 
VAS = VAS+DAS^I2+V&D^» 261 

VLA«VLA*DU I ?I 2 + yBD^,l 28 . 

VKA»VR 
VAE B VA 
PA«VAE 

pam* io 
P A 2 " P A t mui.i 
CAU FUNC lNlPA2/lVMiFUNX) 
VRE»VRA-.i 
PRAsVRE/,005 . 


VPE*VPA»«30625 
PPAsVPE/ 0QA3 
PP1 b ,Q26^PPA 


TRT4> F JN2 y 




A FORTRAN 


bUBRQUTlN 


05/11/73 


PAGE 0005 


RVG*2.738/PVS 
QV0=(PVS-PR1 )/RVO 

C N3 = CN3 ’M < (P C^JL?. )*CN7»i7. )*CMg-CN3)*«l 

a - r , r*wr*i , - 

RySsAVE#(l./CN3)*VlM#( UNLH1 . ) #ANZ + 1 . ) 
PGSbPA-PVS 

RSN»RAR*ARM*ANU#AUf'*PAM*\/IM+Rys#l ,79 


RSMsANU*vIm*PAM*AUM*AMM*RAM 

BFM s PGS/R$M 

qao*bfn+bfm+rbp+{ pa-prahfis 


QRO*QRN*C ( l.-QRF \ #AUH*RVM$BSR*HMD#HPR+QRF*QLO/QLN ) 
QrO s QLO+ ( 3P0-QL0 ) /U 


QV 05 QRo+(QV 3 -QRU)/X 


* ut 

DPA-QRD«QPd 
D AS = QLO-QAO 
OLA«QPa-QtO 



03 
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Program 2 . SU BRO UT I N E AUT O . 




A FORTRAN IV ( VfcR 138 ) S n L)RCF LISTING* AJJC SUBROUTINE Q 5 /U /73 PAGE; 0006 


1 

2 

3 

T 

6 

7 


SUBROUTINE AI-TO 

* 

- — 

REAL 12 

AuTQMD'MIC 


AUKjAJL*AUM,AUNiAUo,AUp,AUQ* 

AUR/ AUS, Ayv,AUX* AUZi AUAj AU6* AU8, A1 DAU, E*C > EXE , 
PA #PA1,PB3.P0T*P20 jSTA.VVR^VV9*Y tl$ 


f %1 J 1 2 i PA ,PAI,PDQ,PO 
""y e iZl " ~ 


CONTROL block 
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Program 3. SUBROUTINE HORMON. 



A FORTRAN 


S£ URCE LISTING* H3RMON 


JBRJUriNi 


PAGE 0007 


1 SUDRDUTjNe HCRMOftltAM » AMC# AMP, AMR/ AMT, AMii ANM*CKB . PA* l» FW, 

2 # AGKjANC*ANP*ANRjANT,ANVjANWjAN1#CNA*CNE # 0FN/ 

3 * I >REK) 

zr 'DTHENSltTN' FOFTtT^T' * " * ' ' ' ' ’ " 

6 C AtAL 1 

7 C ******* ft *### ft* #»><,!< ft * ft ftftft ft $ ## ft ftftftftft ftftftftftftftftftft./ 

IT ALOOSTER3N6 C< NTRPL BLOCK 

1 1 £ftft*ftft«ftftftft ^ftftftftftftfcftftftftftftftftftftftftftftftftftsfcftftftftftAftftftftftftftftftftftftftftftftftftftftft*****#^****#**# 

V2 — r 6 8 AT'TR'-'C'K E 1 c N A7TCT03 5 2 » 9 i * ' ' ~ ‘ 

! 3 1F( AMR,LT.O. ^ AMRsOj. 

A CALL FUNC^N ( PA> AM^/ FONT ) 

5 AMl®ArU + ( ANK*AMP*Af*R-A»U}^Z 

•5 AMC «AMt +rA , Ml-WTCTS ( i « XP ( 7 1 / A FIT) ~ 

7 Ak=20»039-19,8*EXP <-,039l#AMC ) . _ 

Vr\ ^UsThte'^sth' fnirritn t ?V7C*K -- * ” 

W C #*#$**#*#*****#*#* ****#*****##*#** *#***#***#*****##fc****#***#**’ !<,i, ' ! '***’ i ' ,!e * , ’ < ** *#*# 

23 CISE-152.-CNA 

If 1 nr«?? i^7sir.?NS?Jii )*AGK+X, J#REK 

26 ANW®ANW+UANR-1. )* 10 .~ANW>*AMV*I 

27 IF(ANW.LT,.0)ANH»0. 

re isw^atvirva^ ■ , ' ' ' 

29 1F( ANP.GT.100. )AWP=100t 

30 |F(ANP.LT.,QX)ANP S .01 

31 ANl®Ar>l»(ANP»ANl)/Z _____ _ 


-MlNt ■*■ \ Mrti-wm, l * \ J. • -c 'r \ 7 

ANM = 4,0-3.3#FXP(".0967*ANC > 
I F ( ANM. LT » .7) ANMs.7 
RE TURN 
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Program 4. SUBROUTINE BLOOD. 



A l-ORTRAM IV ( VER L38l SOURCE LISTING* BLOOD SUBROUTINE Qlllllll P&C£ ° Q ° 6 

1 SUBROUTINE BLOOD , POT, PDYj PQ 1 > P02 , RC 1 t RC 2, 0# RK£/ 

2 * VB ^VIB/VXE,VIH#VP ,VRC> 

3 _RfcAL__I _ _ _ _ . 

5 C KfcU CELLS AND VISCOSITY BLOCK 


7 1 iLUo5"vIscD5lTY 
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Program 5. SUBROUTINE MUSCLE 






3 PAG 


subroutine mlscle(ALO>amm>agm,aup»a4<,bfm,exc.,hm ,i j* 2 mm -» o§a> 

OVA*OVS*92A#PDn*PKl*PK2#PK3iPMO*PMl#PM3jPM*#PM5/ 
PDE> PQM ^PVQjPgnj OOM jRMQj VPFitS jlb l 

-R^T,~mW 

MUSCLE BLOOD FLOW CONTROL AND PQ2 BLOCK 


mu 



0VA=0SA#HM#5, 

DVS«*0VS*nBFM*0VA-RM0)/HM/5 t /BPM-DVS)/Z6 
PV0*57.XA? n ' /,! 
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Program 6. SUBROUTINE AUTORG. 



A FORTRAN IV (VER L38) SOURCE LISTING! AJTQRG SUBROUTINE qS/U/ 73 PAGE OOlO 


SUBROUTINE AcTaRG(AOM> ARM jARI, AR 2/AR3i AIK j A 2KiA3K*BFN, DD8 ,Hm tit 

* gS8'85y^y^Q2H#PP^PDB,POC,l»DO,PnK,PON,PDR,PO^ 

* T . fl , T PaV>POZjPlD,Q024ROa^Z ,Z4 ,27) 

KtrAt — 


1 
2 

3 

4 - 

5 C 

6 0 NON-MUSCLE OXYGEN DELIVERY BLOCK 

7 C AND NON = MUSCLE LOCAL blood flow control block 

9 c"" W "aUTorIgULATio“7rAPID ' " " " - 

0 C- • — — — , ... . — - m wm ...... 

D$V=QSV+( (BFM#OVA^DaB)/HM/5,/9FN-QSV}/Z7 

■fQV=0|V*57 ,"t4" 


RDO»P0TiH*3, 

DqJS? R ablPof /RDO 
Wg 2-A0 M »02M » t 1.-<B>00Q1 " P10) » »3./512, ) 


Ll 32«QQ2+(D03-MG2)*(1."EXP(-I/Z4) ) 
PQT=QQ2*, 00333 
PlQsPQT 

-I r T P 0 T -rUf-r&r^P-tg-trS-.- 


PgDsPQV-POR 

BrPQB + lPQK^PQD+l » -PQB ) fl 


*r{POB,LT 1 .2)PDBn.2 
ARlsARx* ( POB^AR l ) »( 1 . «E X P' ( '«» 1/A IK ) } 

A ¥ S 1 * A % tl. A M I.L A M n 


ARM»AR1*AR2*AR3 


5 

7 autore gulation, intermediate" ” 7777~777~’ ~ 

9 ^ pgAsp5A+7P0N*p55+i7-paA77z' T 

0 IF (PQA, LT. . 5 }P0A»|5 

AR2«*AR2+7Poi-AR2)H<ht-EXP{«I/A2K 


wm ■ 


,)K 


\l X m AUTDREGULATIDN, long-term 


>5 

-36 


37 


ir9t- 


I F(P002 l 94^j92^l92 


QL 


196 


It in fpfe??:^sh, 3 lf 

*0 A frW - ftfrr < - P O C-AR3)»I/A3K 


RETURN 
END 
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Program 7. SUBROUTINE ADH. 


Quality 


0F Poor 



SUBROUTINE: ADH (AH * AHC> AHK, AHMi AHY, AHZ, AH7* AH0, AUP/ CNA, CN8> CNR/ 
* CNZ * I >PRA,Z) 

__ REAL I 

r ” ** " m — — " ” — . •» “* • • - • J - » • 

£ ANT I DIURETIC HORMONE 
C 


CNBsCNA-CNR 
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Program 8. SUBROUTINE MISC1 . 



A FORTRAN IV tV£R 138) S i.UR C E LISTING! M I SCI 


subroutin 


05/11/73 


PAGE 0012 


1 SUBROUTINE M I SC 1 (AHM,AJ4>AU8, I ,5R , 5RK, $TH, TVD^m* VEC, VlC, VTW> 

2 * VVE/ VV6^VV7 # l) 

3 ^ k fcA_U_ J __ __ __ _ 

5 C fc## #***##*$*#* fc*#*#**###**###*#*##******###***#*#*###*#;*##*####*#* 

f C VASCULAR STRtSS RELAXATION BLOCK 


9 C « * tf**#******##*# *#####** #*#*#####>>:**;{< 


'y i 6sVy6t(SR*(VVE*.3011'*VV77VV6)/Z 
'V7*VV7+VV6 #U»-BXpM/SRK> ) 


THIRST AND DRINKING BLOCK 

# # * * # * # * # * * ❖ * # * # # * * * # >* $ * * # * # # * # * # * * * * * $ * * $ # * # j{< # * * # $ 4c # * # * $ ife * * $ # >j( # # # # * # * * * * * # A * s)t * 


TVZ=» ( .Ol^AHN-,009 ) *STH 
TyDsTVDMTV2-tVDi/Z 
IRTVD.LT, 0. )TVD«0. 


|| £********** ***********^************ ************************************** ******* 
is c autonomjc control block_ 

f|“ C AD APT ' AT I &N '~O f ~B' A ' R P R E CEPTORS 

17 AU4*AU4+AU8*I 

t~8~ » RifT*tfRN “ ' * 

[9 END 
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Program 9. SUBROUTINE HEART. 





A PQRTRA 


TT 



SOUR 



LISTING) H 



UBRQUTINP 05/11/73 PA 


SUBROUTINE HEART ( AMR# DHMj HMD# HR ,PA jPMC*PMI?»PMS,»P0T.,PrA#3A3/ 
* QLO/RTp#SVO,VAEiVL£/VPE^VRE,VVE) 

REAL I „ 

HEART HYPERTROPHY qr DETERIORATION BLOCK 


PAGE 0013 
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Program 10. SUBROUTINE CAPMBD. 






SUBROUTINE C APMBM 5FN# CPC* CP 1 , CPP# DF P* 1 '1 Fp iPC ^PCC^PlF^^LC^PPC^ 

# PRP,Prc> p T$,PTT,P\/3, p VS,RVS,TVD,VG * V ID, Vi F, VP* 

* VPDjVTC* VTO,VTL,VT$,VU0> l ,ZliFUN6) 

- TriMENStStN F tf-nStW * * 

REAL I, IFF 


CAPILLARY ^gMRRANlE DYNAMICS BLOCK 


13C 


PTT»(, VT5/12, >**2‘. 
VIF=VTS~VG 

" M^yNClN <VIF;pTS,FUN6> 


TF s PT T ^ P T 5 
Cr| 5 lPP/Vl F 
PTCs.Z5*CPl 
CPPsPRP/VP 



p5GBR^sS? a 79*BFN 
PC=PVG+PV§ 
PCD»PC+PTC-PPC-PIP 
■ Yfxr^ tr * (CFC * PC0»V 

n i n . 1 . n T r H T T 


vtc) rt 


rL0=7.8+ p IF^PTT 
VTLaVTL+(.O04*PL0-VTL)/Z 
IFTVTt.LT,0.AVlL»0. 
VTpsVTC-V i L~ ^ "P 

u Y ► -1/ *P ^ 1 I < T rx %U T 


vTs-VTS+VTD*! 
vpd=vpd+( tvd-vtc+vtl-yup^dfp-vpo) /Z1 

RETURN 

tm~- 
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Program 1 1 . 


SUBROUTINE PULMON. 





-40^ vr^vjJrjH- 
noil 


3 


5 


SUBROUTINE Pi'LMON(CPF/CPP/CPM,DFPi ! » PCP, PFl , PLaIPlF, PQS, PpA/ PPC> 

* PPD/PPi,#PPNjpp r j,»ppR,iVP ,vpo,vpf ,i *i 3) 

R_cAL i . _ _ 

P^UMUNARY OYNApiICS AND FLUIDS BLOCK 
VP=VP+( VPU*I )/Z3 









no 


Program 12. SUBROUTINE MISC2. 





PAGE 0Q16, 


A FORTRAN IV (V£R L38) SCURCF LISTING! MI5C2 SUBROUTINE 


05/11/73 


SUBROUTINE 
REAL I 


hi sc 2 <hpl*hpr,hsl,hsRj. I , PA> PPAj P0T7STHJ Z 10 / Zll> Z13) 


# 5 $ c1 f r5 f 1 ' 3}C^t n" 

c 

C HEART HYPERTRHPHY OR DETER I QRATIOM BLOCK 
C 

fr frTfr y yyg *r sy * rng i reny >y % # ft ft ft ft ft ft - ft ft ft ft ft ft » ft ft * ft ft ft ft * 

Hp|_*HPL + ( ( ( PA / 100. /HSL>#*Z 13) -HP U *1/57600. 

„ HPr*HPR+( C (PPA/15./HSR)ftftZ13)«HPR)ftI/57600i 

C ft* ft ftftftftftftft ft ft ft ft ft ft ft ft ft ftftftftftftft ft ft ftft; ftftftftftft ft ftftftft ft ftftft ft ftftft ft ft * ft ft ft ftftftft ft ftft ft ft ftft ft ftft ft ftft ft: ft ft ft ft ft ft 


C TISSUE EFFECT ON THIRST AND SALT INTAKE 
C 

C ftftft ftft ft ftft* ftftft ftftft ftft**ftft$ftft ftftftftftft ftftftft ft* ft ftftft ftftft ftftftftftftftftftftftftftftftftftftftftftftft.ftftftftftftftftftftft'ft*' 
ST H-(Z1Q-*P 0 T )* Z11 


IF(STH.L,T.I,)STH»I, 
I Fl.STH p GT. 8 * ) STH?8 . 


RETURN 
-EitD 
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Program 13. 


SUBROUTINE PROTEN. 



SUBROUT I Nt P f *OTe\'(CHY,CPG^CPI,CPK^CPPj»CPF 

* DPO/DPY/GPD,GPR/t # IFF 

* VTl,Z /PHD) 

R’Ffll ‘ HT F"P / LFK ' 

TISSUE FLUIDS/ PRESSURES AND GEL BLOCK 


PLASMA AND T I S$UE . FLUI D PRQTg I \] 

rpPL^DPLMVT^CPUoSLT/r" 

“ I F (TC, LT«0"« ) PC *0 

DPC=DPC+(CPK*(CPP-CPI )#PC**PCg-DPC) /Z 
CPl=DPC-DPL 

DLZ=LPK«(CPR-CPP) _ 


iMUHF.U' .O'!* JI-*LZsA,^UL4 

DLP=DLP+(DLZ-0LP5/Z 

PRPsPRP + ^LP-DPa+DPL-DPC-PPQ)*! 


PGX = CHY 5 > Cs ! t 2**0l332 , ! s CPG+CPG 
GPDsGPDM ,0005*<CPI-PGX)*VG»GPD)/Z 

uPRoGPK + gP I 

IFP»IHP+(1»PI-GPD)*I 

RETURN 

END 
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Program 14. SUBROUTINE KIDNEY. 



A P-DRTRAN IV (VER L38> SjURC? LISTING* KIDNEY SUBROUTINE 


05/11/73 


PAGE 0018 


1 

2 

3 

4 

5 C 

6 C 

7 C 

« — m- 


SJBR0UT1NE < JDNEY(AAR/AHM/AM , APD, A RF, AUM, CNE j C NX ICNY/ G&{,/ GFN/ GF R/ 
* GF 2 iGF 3 #GF 4 ,rAP^ I_ vNAE^ ^OiNIO,NlnD>NpijPA /PAR/ 


REAtr l/MAfc,h f -DVMTDiSb^/-Nai RBF ' ,REK, - FM,R ~'-“- 
MONEY DYNAMICS AND excretion block 


r • n a v m 'Li*v/ 1 ; w -j ./ _ 

»trp*vi m /VuD/Z) 


-GfTgt l - GPN T- , 1? 3 - t 

IF(GF3 t GT,l5, )GF?sl5. 
IMGF3.LT,.4)GF3 k.4 

AAR=31 ♦ 67*V 1 'i* ( AUM*ARF + 1 t ■*ARF ) #GF3 



RFN=PAR/RR 

RBF*REK#rI-N 

~AP trg , AA ' R ~ *R PN 


T 


9 

0 

1 
t 

3 

4 
_5 
Wr 

2 

9 


■w 


GLPsPAR-APD 
^FL b GLP-PPC-18. 


p 

G 


IT 


?Fl=GFN 

\b‘ 

r rri 

VtlD 1 

IF(VUD.LT,,OOO2)VU0«.OOO2 


Ns-frEttr fpf tr^rOD 7 81r= GfH^->^S-F-27 ,_ Z: 

<M^r 3ni - CT ” ooilGD T3 142 


J6 
27 

if c 

30 

31 
-§ 


KIDNEY SALT Cu 
-A t S3- Ertfr£ 


^R QL VT E5 ANp T CEU AR ATbft BLOCK) 


NDZ«ipOO,*VUD/AM/ CNE/CNX+CNY 
NQD*NUD+ (FjOZ-NOD J / Z 

NE 'D -Ni ' D^S * H** 1, OD 

NAE=NAE + in£D*I 

RETURN 

EhD 


3 

34 

35 
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Program 15. 


SUBROUTINE IONS. 


tiftiGlNAt 

poo n 


Wfe m 

quality 



RCP 


SUBROUTINE K-NS (AM t CCD/ CKE,C.K I » CNtt/ I , 
* KIR/KPO/NAE/REK/VEC/ VIC/ 

REAL. I/KCU/Kr./KEr, KI/i<IPjKIEjKIR, KOD/ NAE 

ELECTROLYTES ANO cell water block 

VECsVTS+VP+yrF 


t.KhSK.t/ Vfcw 

KCJD- C ,OOOA2#CKE+.OQOl4*AM#CKE)#FEK 
KIR*2§50« + lAC« s f‘CKE 
K I E *K 1 R-K I 


U+USJC'i'.Un-K.tUJ/i 

•KCD#I 
d-kCd«koo 

•KED*I 


7 

C 

8 

C 

9 

V 












118 


Program 16. SUBROUTI NE GELFLD. 





A FORTRAN IV (VfcR L38) SOURCE L 1 S T 1 NG » GELFLD SUBROUTINE Q5/U/73 PAGE 0020 

SUBROUTINE QPLFirKCHV^CPG/CPl/CPRiHYUlFPjPGCiPGH^PGPjPGR^ PQXfPlF* 

* t KtT „ „ „ , PRM,PTC#PTS,PIT/VG iVGD,VIF,VRS,VTSxV2D,FUN6) 

DIMbhSIuN FJ*’6<14) 

Stf rFi 3 — - - - 


HC 


GEL FLUID DYNAMICS 


WHY* 2».2 


PGR*H*CHY 
CPGsGPR/VG 
PGP s « 25 ^PGX 

-Rnjt-=rP&P>pGR 

VlFsVTS-VG 

CALI FUNCTN ( VI F j PTSj» FUN6 ) 


P I F*>P [T-PTs 
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Appendix A - Glossary of Terms 

The following lisf includes all variables used in the model together with the normal 
values of these variables. Independent variables (never calculated by the program) are 
indicated by *. Units used are: volume in liters, mass in grams, time in minutes, chemical 
units in milliequivalenfs, pressure in millimeters of mercury, and control Factors as ratio 
to normal. 

AAR- afferent arteriolar resistance (31.7) 

AGK*- constant concerned with effect of renin on angiotensin formation (0.20) 

AH- antidiuretic hormone secretion rate (3.0) 

ARC- antidiuretic hormone concentration (1.0) 

AHK*- constant used in calculating antidiuretic hormone concentration (7.0) 

AHM- antidiuretic hormone multiplier (1.0) 

AHY- adapted effect of right atrial pressure on antidiuretic hormone secretion rate (0.0) 

A HZ- basic effect of right atrial pressure on antidiuretic hormone secretion rate (0.0) 

AH8- effect of autonomic stimulation on antidiuretic hormone secretion rate (0.0) 

ALO*- maximum aortic arterial oxygen saturation (1.0) 

AM- aldosterone multiplier (1.0) 

AMC- aldosterone concentration (1.0) 

AMM- muscle vascular constriction caused by local tissue control, ratio to resting state (1.0) 
AMP- effect of arterial pressure on rate of aldosterone secretion (1 .0) 

AMR- effect of sodium to potassium ratio on rate of aldosterone secretion (1.0) 

AMT*- time constant of aldosterone accumulation and destruction (60) 
rate of aldosterone secretion (1.0) 


AMI- 
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ANC- angiotensin concentration (1.0) 

ANM- angiotensin multiplier effect on vascular resistance, ratio to normal (1.0) 

< 

ANP- effectLof-renal blood flow on angiotensin formation (1.0) 

ANR- effect of glomerular filtration and sodium concentration on renin formation 
with consequent effect on angiotensin formation (1.0) 

ANT*- time constant of angiotensin accumulation and destruction (15.0) 

ANU- non-renal effect of angiotensin (1.0) 

ANV*- constant used in calculating effect of renin formation on angiotensin formation 

(0.0003) 

ANW- partial effect of renin on angiotensin formation (0.0) 

ANY*- constant used to calculate angiotensin effect on venous volume (-0.2) 

ANZ*- constant used to calculate angiotensin effect on venous resistance (0.4) 

AN! - rate of angiotensin formation (1.0) 

AOM- autonomic effect on tissue oxygen utilization (1 .0) 

APD- afferent arteriolar pressure drop (38.0) 

ARF*- intensity of sympathetic effects on renal function (1.5) 

ARM- vasoconstrictor effect of all types of autoregulation (1.0) 

AR1- vasoconstrictor effect of rapid autoregulation (1.0) 

AR2- vasoconstrictor effect of intermediate autoregulation (1.0) 


AR3- vasoconstrictor effect of long-term autoregulation (1 .0) 
AU- overall activity of autonomic system (1 .0) 

AUB- dffect of baroreceptors on autoregulation (1 .0) 

AUC- effect of chemoreceptors on autonomic stimulation (0.0) 
autonomic stimulation of heart (1.0) 


mt4 


AUH- 
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AUJ- 

AUK*- 

AUL*- 

AUM- 

AUN- 

AUO- 

AUP- 

AUQ*- 

AUR- 

AUS*- 

AUV*- 

AUX*- 

AUY*- 

AUZ*- 

AU4- 

AU6- 

AU8- 

AVE- 

A1B- 

A1K*- 

A2K*- 

A3K*- 

A4K*- 


basic overall autonomic stimulation (1.0) 

time constant of baroreceptor adaptation (0.0005) 

sensitivity of sympathetic control of vascular capacitance (0.21) 

sympathetic vasoconstrictor effect on arteries (1.0) 

effect of CNS ischemic reflex on autoregulation (0.0) 

fractional departure of overall activity of autonomic system from normal (0.0) 

autonomic stimulation of peripheral circulatory sensitivity (1.0) 

sensitivity of sympathetic control of peripheral circulation (1.0) 

autonomic stimulation for heart rate (1.0) 

sensitivity of sympathetic control of heart rate (1.0) 

sensitivity of sympathetic control on heart function (0.3) 

sensitivity of baroreceptors (3.0) 

sensitivity of sympathetic control of veins (0.25) 

overall sensitivity of autonomic control (1.0) 

degree of adjustment of baroreceptor response (0.0) 

adapted baroreceptor response (1.0) 

rate of adaptation of baroreceptors (0.0) 

effect of autonomic stimulation on venous resistance (1 .0) 

sensitivity parameter for baroreceptor drive (1 .0) 

time constant of rapid autoregulation 

time constant of intermediate autoregulation (20.0) 

time constant of long-term autoregulation (11520.0) 

time constant for muscle local vascular response to metabolic activity (1.0) 
muscle blood flow (1.0) 


BFM- 
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BFN- 

CCD- 

CFC*- 

CHY- 

CKE- 

CKi- 

CNA- 

CNB- 

CNR*- 

CNE- 

CNX* 

CNY*- 

CNZ*- 


CN2‘ 

CN3- 

CN7*- 

CPF*- 

CPG- 

CPI- 

CPK*- 


blood flow in non-muscle, non-renal tissues (3.0) 
concentration gradient across cell membrane (0.0) 
capillary filtration coefficient (0.007) 
concentration of hyaluronic acid in tissue fluids (5.0) 
extracellular potassium concentration (5.0) 
intracellular potassium concentration (142.0) 
extracellular sodium concentration (142.0) 

difference between extracellular sodium concentration and set point used to 
calculate antidiuretic hormone secretion rate (3.0). 

reference sodium concentration used in determining effect of sodium on anti- 
diuretic hormone secretion rate (139.0). 

sodium concentration abnormality causing third factor effect (10.0) 

constant used in calculation of renal excretion rate of sodium (2.5) 

constant used in calculation of renal excretion rate of sodium (6.0) 

sensitivity of antidiuretic hormone production rate to extracellular sodium 
concentration (1 . 0) 

constant used in calculation of venous resistance (0.0212) 

dummy variable used in calculation of the effect of capillary pressure on 
venous resistance (0.366) 

constant used in calculation of venous resistance (0.2). 

sensitivity of rate of transfer of fluid across pulmonary capillaries to pressure 
gradient (0.0003) 

concentration of protein in tissue gel (12.5) 

concentration of protein in free interstitial fluid (16.5) 

rate constant used in determining loss of plasma protein through systemic 
capillaries ( 1.6x10"^) 



CPN- 
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concentration of protein in pulmonary fluids (30.0) 

CPP- plasma protein concentration (70.0) 

-CPR*- reference plasma protein concentration governing protein production by liver 
(85.0) 

CV*- venous capacitance (0.0925) 

DAS- rate of volume increase of systemic arteries (0.0) 

DAU- autonomic stimulation drive (1.0) 

DFP- rate of increase in pulmonary free fluid (0.0) 

DHM- rate of cardiac deterioration caused by hypoxia (0.0) 

DLA- rate of volume increase in pulmonary veins and left atrium (0.0) 

DLP- rate of formation of plasma protein by liver (0.007) 

DLZ- undamped plasma protein concentration differential causing protein production 
by liver (0.007) 

DOB- rate of oxygen delivery to non-muscle cells (180.0) 

DPA- rate of increase in pulmonary volume (0.0) 

DPC- rate of loss of plasma proteins through systemic capillaries (0.05) 

DPI- rate of change of protein in free interstitial fluid (0.0) 

DPL- rate of systemic lymphatic return of protein (0.05) 

DPO*- rate of loss of plasma protein (0.007) 

DRA- rate of increase In right atrial Volume (0.0) 

DVS- rate of increase in venous vascular volume (0.0) 

EXC* exercise activity, ratio to normal at rest (1.0) 

EXE- exercise effect on autonomic stimulation (0.0) 

constant concerned with effect of muscle cell PO 2 on aiJ tonomic stimulation 
during exercise (3.0) 


EX1*- 
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FIS*- fistula parameter (0.0) 

GBL*- Goidblatt hypertension parameter (0.0) 

GFN- glomerular filtration rate of undamaged kidney (0.125) 

GFR- glomerular filtration rate (0.125) 

GF1 - value of GFN on previous iteration (0. 125) 

GF2*- constant used in calculation of glomerular filtration rate (0.05) 

GF3- degree of autoregulatory feedback at macular densa (1 .0) 

GF4*- constant controlling the feedback loop for GF3 (5.0) 

GLP- glomerular pressure (62.0) 

GPD- rate of increase of protein in gel (0.0) 

GPR- total protein in gel (143.0) 

HKM*- constant used in calculation of portion of blood viscosity caused by red blood 
cells (0.53) 

HM- hematocrit (41 . 0) 

HMD- cardiac depressant effect of hypoxia (1 .0) 

HMK*- constant used in calculation of portion of blood viscosity caused by red blood 
cells (90.0) 

HPL- hypertrophy effect on left ventricle (1.0) 

HPR- hypertrophy effect on right ventricle (1.0) 

HR- heart rate (72.0) 

HSL*- basic left ventricular strenght (1.0) 

HSR*- basic right ventricular strength (1.0) 

HYL*- quantity of hyaluronic acid in tissues (57.0) 

I- integration step size (0.73) 

IFP- interstitial fluid protein (9.1) 

II- variable integration step size utilized on stable-asymptote 
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12*- normal increment on time (0.003) 

13*- maximum time increment For stable asymptote (20.0) 

KCD- rate of change of intracellular potassium concentration (0.0) 

KE- total extracellular fluid potassium (75.0) 

KED- rate of change of extracellular potassium concentration (0.0) 

Kl- total intracellular potassium concentration (3550.0) 

KID* rate of potassium intake (0.0028) 

KIE- excess potassium concentration causing change in intracellular potassium 
level (0.0) 

KIR- total expected level of potassium in the intracellular fluid under equilibrium 
conditions (3550.0) 

KOD- rate of renal loss of potassium (0.0028) 

LPK* - rate constant for plasma protein production by liver (0.00047) 

LVM- effect of aortic pressure on left ventricular output (1.0) 

MMO- rate of oxygen utilization by muscle cells (60.0) 

M02- rate of oxygen utilization by non-muscle cells (180) 

NAE- total extracel lular sodium (2136.0) 

NED- rate of change of sodium in extracellular fluids (0,0) 

NID*- rate of sodium intake (0. 1) 

NOD- rate of renal excretion of sodium (0.1) 

NOZ- effect of urinary output, aldosterone, and sodium level on renal excretion rate 
for sodi um (0.1) 


OMM*- muscle oxygen utilization at rest (60.0) 
OSA- aortic oxygen saturation (1.0) 



OSV- 


non-muscle venous oxygen saturation (0.7) 
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OVA- oxygen volume In aortic blood (203.0) 

OVS- muscle venous oxygen saturation (0.7) 

02A*- sensitivity of the effect of autonomic stimulation on metabolism (1 .5) 

02M*- basic oxygen utilization in non-muscle body tissues (180.0) 

PA- aortic pressure (100.0) 

PAM- effect of arterial pressure in distending arteries, ratio to normal (1.0) 

PAR- renal arterial pressure (100.0) 

PA1- effective pressure drive on autonomic system (100.0) 

PA2- effective arterial pressure on left ventricle (100.0) 

PC- capillary pressure (18.4) 

PCD- net pressure gradient across capillary membrane (0.45) 

PCE*- capillary pressure exponent (3.0) 

PCP- pulmonary capillary pressure (7.0) 

PDO- difference between muscle venous oxygen POo an< ^ normal venous oxygen Pqo 

( 0 . 0 ) 

PFI- rate of transfer of fluid across pulmonary capillaries (0.0) 

PFL- renal filtration pressure (16.0) 

PGC- colloid osmotic pressure of tissue gel (6. 1 ) 

PGH- absorbency effect of gel caused by recoil of gel reticulum (-4.0) 

PGL- pressure gradient in lungs (15.2) 

PGP- colloid osmotic pressure of tissue gel caused by entrapped protein (4. 13) 

PGR- colloid osmotic pressure of interstitial gel caused by Donnan equilibrium (2.0) 

PGS- pressure difference between arteries and veins (96.0) 
venous pressure gradient (3.7) 


PGV- 
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PGX- activity factor for protein in the interstitial fluid (16.5) 

PIF- interstitial fluid pressure (-6.0) 

PK1*- constant used in calculating muscle cell Pqt ^ rom total volume of oxygen in 

muscle cells (2500.0) 

PK2*- constant used in calculating muscle cell Pq« from total volume of oxygen in 
muscle cells (800.0) 

PK3*- constant used in calculating rate of oxygen transport to muscle cells (2.0) 

PLA- left atrial pressure (0.0) 

PLD- pressure gradient to cause lymphatic flow (0.8) 

PLF- pulmonary lymphatic flow (0.0003) 

PMC- mean circulatory pressure (6.9) 

PMO- muscle cell PO 2 (8 • 0) 

PMP- mean pulmonary pressure (4.6) 

PMS- mean systemic pressure (7.25) 

PM1- effective muscle ce II P 0 , (8-0) 

PM3*- minimum value allowed for PM1 (0.001) 

PM4*- constant used in calculating rate of oxygen transport to muscle cells (-1.0) 

PM5*- constant used in calculating rate of oxygen transport to muscle cells (122.0) 

POA- rate of change of intermediate autoregulation vasoconstrictor effect (1.0) 

POB- rate of change of rapid autoregulation vasoconstrictor effect (1 .0) 

POC- rate of change of long-term autoregulation vasoconstrictor effect (1.0) 

POD- non-muscle venous P 02 minus normal value (0.0) 

POE- sensitivity control for oxygen feedback control loop (1.0) 

POK - sensitivity of rapid system of autoregulation (0.06) 



POM*- 
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sensitivity of oxygen feedback control loop (0.08) 

PON*- sensitivity of intermediate autoregulation (0.3) 

POQ- effective non-muscle cell PO 2 (8.0) 

POR*- reference value of capillary Pq 2 ' n non-muscle tissue (40.0) 

POS- pulmonary interstitial fluid colloid osmotic pressure (12.0) 

POT- non-muscle cell P 02 (8.2) 

POV- non-muscle venous PO 2 (40.0) 

POY*- sensitivity of red cell production (0.0000464) 

POZ*- sensitivity of long-term autoregulation (0.3) 

POl*- constant used in determining oxygen deficit factor causing red cell production (8.25) 
P02- oxygen deficit factor causing red cell production (0.25) 

PPA- pulmonary arterial pressure (15.4) 

PPC- plasma colloid osmotic pressure (28.0) 

PPD- rate of change of protein in pulmonary fluids (0.0) 

PPI- pulmonary interstitial fluid pressure (-10.0) 

PPN- rate of pulmonary capillary protein loss (0.0) 

PPO- pulmonary lymph protein flow (0.009) 

PPR- total protein in pulmonary fluids (0.38) 

PPI- variable used to empirically relate pulmonary arterial pressure and pulmonary 
arterial resistance (0.4) 

PP2- effective pulmonary arterial pressure (15.5) 

PRA- right atrial pressure (0.0) 


PRM- 


pressure caused by compression of interstitial fluid gel reticulum (-5.0) 
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PRP- total plasma protein (208.0) 

PR1- effective right atrial pressure (0.0) 

PTC- interstitial fluid colloid osmotic pressure (4.1) 

PTS- solid tissue pressure (7.0) 

PTT- total tissue pressure (1.0) 

PVG- venous pressure gradient (14.6) 

PVO- muscle venous P 02 (40.0) 

PVS- average venous pressure (3.8) 

PIO- tissue P 02 effective in oxygen utilization (8.0) 

P20- muscle cell PO 2 effective in depressing rate of metabolism (8.0) 

QAO- blood flow in the systemic arterial system (5.0) 

QLN- basic left ventricular output (5.0) 

QLO- output of left ventricle (cardiac output) (5.0) 

QOM- total volume of oxygen in muscle cells (2400.0) 

Q02- non-muscle total cellular oxygen (2400.0) 

QPO- rate of blood flow into pulmonary veins and left atrium (5.0) 

QRF*- feedback effect of left ventricular function on right ventricular function (0.6) 

QRN- basic right ventricular output (5.0) 

QRO- actual right ventricular output (5.0) 

QVO - rate of blood flow from veins into right atrium (5.0) 

RAM*- basic vascular resistance of muscles (96.3) 

RAR*- basic resistance of non-muscular and non-renal arteries (30.52) 

RBF- renal blood flow (1.2) 
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RCD*- rate of change of red cell mass (0.0) 

RC1- red cell production rate (0.000011) 

RC2- red cell destruction rate (0.000011) 

RDO- resistance of diffusion of oxygen from capillaries to cells (555.0) 

REK*- fraction of normal renal function (1.0) 

RFN- renal blood flow if kidney is not damaged (1 .2) 

RKC*- rate constant for red cell destruction (5.8 x 10'^) 

RMO- rate of oxygen utilization by tissues (60.0) 

RPA- pulmonary arterial resistance (1.6) 

RPT- pulmonary vascular resistance (3,0) 

RPV- pulmonary venous resistance (1.4) 

RR- renal resistance (84.0) 

RSM- vascular resistance in muscle (96.5) 

RSN- vascular resistance in non-muscle/ non-renal tissues (32.5) 

RTP- total peripheral resistance (19.4) 

RVG- resistance from veins to right atrium (0.72) 

RVM- depressing effect of pulmonary arterial pressure on right ventricle (1 .0) 
RVS- venous resistance (2.8) 

SR*- intensity factor for stress relaxation (0.5) 

SRK*- time constant for stress relaxation (33.0) 

STA*- overriding value of overall activity of autonomic system AU (0.0) 

STH- effect of tissue hypoxia on salt and water intake (1 .0) 

SVO- stroke volume output (0.07) 
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T- total time elapsed 

TRR- tubular reabsorption rate (0. 124) 

TVD- rate of drinking (0.007) 

TVZ- combined effect of tissue ischemia and central nervous stimulation on thirst 
and drinking (0.001) 

Tl- total time elapsed on previous step 

U*- damping factor for QPO (4.0) 

VAE- excess volume in systemic arteries that causes stretch of arterial walls (0.354) 

VAS- volume in systemic arteries (0.85) 

VB- blood volume (5.0) 

VBD- volume correction factor added to systemic circulation to allow for updating 

blood volume (0.0) 

VEC- extracellular fluid volume (15.0) 

VG- volume of interstitial fluid gel (11.5) 

VGD- rate of change of tissue gel volume (0.0) 

VIB- blood viscosity, ratio to that of water (3.0) 

VIC- cell volume (25.0) 

VID- rate of fluid transfer between interstitial fluid and cells (0.0) 

VIE- portion of blood viscosity caused by red blood cells (1.5) 

VIF- volume of free interstitial fluid (0.55) 

VIM- blood viscosity, ratio to normal (1.0) 

VLA- volume in left atrium (0.40) 

VLE- excess volume in left atrium causing stretch of left atrium and pulmonary veins (0.0) 
VP- plasma volume (3.0) 
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VPA- volume in pulmonary arteries (0.38) 

VPD- rate of change of plasma volume (0.0) 

VPE- excess volume in right atrium causing stretching of the right atrium (0.07) 

VPF- pulmonary free fluid volume (0.012) 

VRA right atrial volume (0.1) 

VRC- volume of red blood cells (2.0) 

VRE- excess volume in right atrium causing stretching of the right atrium (0.0) 

VTC- rate of fluid transfer across systemic capillary membrances (0.0) 

VTD- rate of volume change in total interstitial fluid (0.0) 

VTL“ rate of systemic lymph flow (0.003) 

VTS- total interstitial fluid volume (12.0) 

VTW- total body water (40.0) 

VUD- rate of urinary output (0.001) 

VVE- excess venous vascular volume before stress relaxation correction (0.33) 

VVR- volume of blood in veins at zero venous pressure (2.95) 

VVS- venous vascular volume (3.0) 

VV6- rate of change of vascular stress relaxation effect (0.0) 

VV7- increased vascular volume caused by stress relaxation (0.0) 

VV8- excess volume of blood In the systemic veins after stress relaxation correction (0.31) 
VV9*- reference venous vascular volume (3. 159) 

V2D*- resistance factor which converts pressure drop to rate of change of tissue gel 
volume (0.02) 

X*- damping factor for QVO (10.0) 
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Y*- damping factor for DAU (1.0) 

Z*~ damping factor for AH, DAU, DFP, DLP, DPC, DPL, GFN, GPD, KCD, NOD, 

POA, POB, PPD, TVD, VID, VTC, VTL, VUD, VV6 (1.0) 

Zl*~ damping factor for VPD (1 .0) 

Z3*- damping factor for VP (4.0) 

Z4*- fime constant used to calculate non-muscle cell total cellular oxygen (10.0) 

Z5*- time constant used to calculate volume of oxygen in muscle cells (10.0) 

Z6*- damping factor for OVS (5.0) 

Z 7*- damping factor for OSV (5.0) 

Z8*- time constant of autonomic response (1.0) 

Z10*- constant used to calculate effect of tissue hypoxia on salt and water intake (8.25) 

Zll*- constant used to calculate effect of tissue hypoxia on salt and water intake (4.0) 

Z12*- constant that converts exercise activity to autonomic stimulation (1.24) 

Z13*- constant used in calculating heart hypertrophy (0.625) 



PROGRAM DESCRIPTION GUIDE 


A. IDENTIFICATION 
Program Name 
Programmer’s Names 
Programmer Contact 
Date of Issue 


- Guyton 

- Guyton, White, and Marks 

- V . J. Marks, GE/AGS, Houston 

- April 16, 1973 


B. GENERAL DESCRIPTION 

This model presents a systems analysis of human circulatory regulation 
based almost entirely on experimental data and cumulative present know- 
ledge of the many facets of the circulatory system. The model itself 
consists of 18 different major systems that enter into circulatory 
control. These systems are grouped into l6 distinct sub-programs that 
are melded together to form the total model. 

In spite of the fact that the total model contains almost 100 independ- 
ent variables and over 350 mathematical relations of various types, each 
major system is modeled in a relatively crude way only, with emphasis 
placed on gross correctness, not fine details. It has been found that 
the systems analysis thus developed is successful in predicting the out- 
come of many varied stress experiments. This is only possible because 
of the extreme stability and many built-in compensations of the actual 
circulatory system. Without this inherent stability, each system would 
have to have been modeled in a much more detailed fashion to produce the 
requisite correlation with experiment. 

The model develops circulatory regulation and fluid regulation in a 
simultaneous manner. Thus, the effects of hormonal and autonomic control, 
electrolyte regulation and excretory dynamics are all important and are 
all included in the model. The model does not treat respiration ot 
thermal regulation. 


C. USAGE AND RESTRICTIONS 

Machine and Compiler Required 

Peripheral Equipment Required 

Approximate Amount of Memory 
Required 


XEROX Sigma 3, ANSI Fortran 

Card Reader, Printer, Teletype 

Guyton (Model A) - + 32 AA 
Guyton (Model B) - +393C 
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D. PARTICULAR DESCRIPTION 


Equations used - See the following reference. 

Guyton, A.C., Coleman, T.G., and Granger, H.I., "Circulation: 
Overall Regulation," Annual Review of Physiology , V. 3^: 13-^6, 
1972. 

Definitions of Terms - Appendix A 
Values of Variables - Appendix B 


E. DESCRIPTION OF INPUT 


1. Machine Control Cards 
1J$B 

1 FORTRAN 

Source Cards (See Appendix C for listing of Guyton Model A) 
(See Appendix D for listing of Guyton Model B) 

!E0D 

I0L0AD 

I$R00T 512,, G0 
:$mp 

I$END 

1XEQ 

Data Cards (See Appendix B for printout of input variables) 


2 . Data Cards 

The GUXTON MODEL A reads data variables as illustrated in Appendix B. 


Column 

Format 

Description 

1-13 

E13.6 

Variable Value 

lh-15 

2X 

Blank 

16-20 

15 

Array location (stop reading of input 
data if less than l) 

21-22 

2X 

Blank 

23-26 

Ah 

Variable Name 
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The GUYTON MODEL A utilizes the teletype to input initial data, 
modify specific data, and output requested data. See Subroutine 
INPUT in Appendix C for complete explanation. 

The GUYTON MODEL B reads initial data using the same format as 
does the GUYTON MODEL A, but does not require as many input 
variables because of internal initialization of some variables. 
See Appendix E for those variables required. The GUYTON MODEL B 
does not interact with the teletype and thus requires additional 
data cards. 


Column 

Format 

Description 

Card A+l 



1-80 

20A4 

Variable names of required output 
variables* If columns 1-3 contain 
ALL, then all variables will be 
printed as shown in Appendix E. 

Card A +2 



1-5 

15 

Model time for next output printing 
change of time units, or change of 
specified variables. 

6 

IX 

Blank 

7-10- 

a4 

Model units of time (SECS, MIN, 
HOUR, or DAYS) 

Card A+3 to A+N 



1-5 

15 

Same as card A+2 

6 

IX 

Blank 

7-10 

a4 

Same as card A+2 if units of time 
require changing. Blank to change 
variables values. Same as last 
card to continue run. 

11-14 

A4 

Variable name for which value 
requires changing. 

15-27 

E13.6 

New value of variable 

If columns 7-14 

are blank the 

program will stop. 
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F. DESCRIPTION OF OUTPUT 

See Appendix B for example of GUYTON MODEL A output. See Appendix E 
for example of GUYTON MODEL B output. 


G. INTERNAL CHECKS AND EXITS 

Curve limits are checked with a diagnostic message being printed if 
they are exceeded. 

The GUYTON MODEL B checks input data for invalid requests and exits 
when it finds one. 


H. INDEPENDENT SUBROUTINES 

See Appendix C for listing of all subroutines required by the GUYTON 
MODEL A. 

See Appendix D for listing of all subroutines required by the GUYTON 
MODEL B. 


I. SYSTEM SUBROUTINES 

No special system subroutines required. 


J. COMPLETION OR FINAL CHECKOUT DATE 


3/10/73 
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APPENDIX A 
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FiCL*k 1 Syrtem-analy-Kdiagramiur n s uhitiu (loft he circulation l nil* ire the following volume 
In liters, mas* In gram*, time in minute* cliemical umu m milli«iulv alerts, pr entire in niilhmetr r*. of 
mercury , control factor* in arbitrary unit* but m most Instance* etpres.-*d a* the ratio to normal— for 
icstawr a \ aliic of l rrjxrr^rnl* nerni it Normal v alue* are given on the line* that represent the respec- 
tive variable*. 

Tbe following U a lut of the imjioriunt dependent and indttwr.dent v*mb1e* in ll* analysis (addi 
iwtud variable* are iirwent for tmr)o*e* of calculation but generally hav e no phy mologuxd significance) 

BFB — blood floor In non mu*de non-renal Ume* 


A A R— afferent arteriolar resistance 
ABM — antidiuretic hormone multiplier ratio of 
normal effect 

AW— aldosterone multiplier, ratio of normal 
effect 

AMC — aldosterone concentration 
AMU — mtnele vascular constriction ciu<*vl by 
local tissue control, ratio to renting state 
A W?- effect of arterial jumwure on rate of aldo- 
sterone secretion 

AMR — effect of sodium to potassium ratio on 
aldosterone secretion rate 

AMT — lime constant of aldosterone accumulation 

aiv! destruction 

ABC — angiotensin concentration 
AVlf — angiotensin niultljdler effect on vascular 
resistance ratio to normal 
A B A' — effect ef sodium concentration on rate of 
angiotensin formation 

AtfP— effect of renal blood flow on aoglotensln 
formation 

ANT — time constant of aa*toten*]n«ccu mutation 
and destruction 

ANV— nomenal effect of angiotensin 
AOU — autonomic effect on tissue oxygen utillan- 
Uoa 

APB — afferent arteriolar pressure drop 
ART— Intently of sympathetic effects on renal 
function 

ARM — vasoconstrictor effect of all types of 
au in regulation 

ARI — vasoconstrictor effect of rapid uutoreguU 
tJQO 

A R2-—\ a.'ocorutncioi’ effect* of mtermedUte 

uuioregulaUun 

ARI -vasocon>inctor effect of long term aoto- 
reyruUuon 

A(J — overall activity of autonomic *y stem ratio 
to normal 

AOU -effect of barorecctitors <m aumreRUlatlon 
AV C -effect of cheiimivcetitor* on autonomic 
stimulation 

AUU — autonomic sti mill (turn of heart ratio to 
normal 

■tdff~llmc«jii<aiii of buroreccplor adapt nthm 
AUf -sendllvlty of sy mitailteiir control of 
vtwcul ir catwt itaoce 

AU M — .ympatbeiic vaW)n»lrKlor effect on 
ancrics 

A t\V- -effect of t Khriuic reflex on atuo- 
regglulicm 

AUV — etiuiiuiy conirid «f luiunomka on Iieart 
time turn 

A V >*— sensitivity of mpathrlic control of veins 
4 1*7 inmU*rnuiiiity of autonomic control 
11/ •symtotheiii va-ocxim-triilur effect on 
veins 

AIK- -tune constant «f npul autoregulation 
\2K — time ccn-t mt of intermediate autoregu 
latioft 

AiJC — time cun-4 ml of long-term miorcgulatlon 
A4K — time constant for nm-clc local va-cular 
re*i*»n-e to metabolic nruvity 
JtFM — ium-< 1** 1>W1 lt*m 


CA— capacitance of systemic arteries 
CCD — concentration gradient across cdl mem- 
brane 

CU Y — concentration of hyaluronic add in tissue 
fluid* 

CKE — extracellular potassium concentration 
CK/ - IntnjctlluUr ixst&^um concentration 

CtfA— extracellular sodium concentration 
CNE — sodium concentration abnormality cau-uni 
third factor effect 

CPC — concentration of protein in tissue gd 
CP/— concentration of protein in free Interstitial 
fluid 

CPB — concentration of protein In pulmonary 
fluids 

CPP — plasma protein concentration 
CV — venous capacitance 
DA $*— rate of volume Increase of **y*temJc arteries 
DPP^rate of incretu* hi indmonary free fluid 
DBM — rate of cardiac deterioration rained by 
hypoxia 

DLA — rate of volume Increase in pulmonary vein* 
and left atrium 

DLP — rate of formation of plasma protein by 
Uvtr 

DOB — rate of oxygen delivery la non muscle cell* 
DP A — rate of Increase in pulmonary volume 
DPC — rate of loss of plasma proteins through 
systemic capillaries 

DPI — rate of change of protein la free Interstitial 
fluid 

DPL-~ rate of systemic lymphauc return of 
protein 

DPO— rate of loss of plasma protein 
BRA — rate of increase m right atrial volume 
DVS — rate of increase In venous vascular volume 
EVR — postidonierubr re*i*tance 
£ YC*-e\erCl<e activity ratio to activity at rent 
fflffi^evcrci* effect on autonomic i-tumil'Uion 
GFN — klomerulur filtration rate of undamaged 
kidney 

GFR'" ktomcru|-\r filtntion rate 
Of P-— glaiiurular prewire- 
OPD— rate oi increase of protein m 
f? PR tot il protein m cel 
HM 'hematocrit 

HMD -l trdtoc <lcj>rc*<int effect of hypoxia 
H PI -hypertrophy effect on Irft ventricle 
II PR— h>|*rtruph> effect on heart ratio to 
noTintl 

NR -he-m rite 

-ba-oeWt vcntmtilar -trentjih 
H'oR—U i-'tc -tren^th of i-ujit ventricle 
If l i ~*f|U I nifty ol hyaluronic acid iu 
IFP— mter-uti il lhn<l imitein 
K( D—tj r uf ell llUf ol |kiIv«iutq concentmtlon 

KD ~ -tot il i-urin Hid ir fluid itoti—WM 
KFD -rue of clu«u k e of extrareliul ir thud eon 
<t mr itKitt 

KI -lut d imr io liul ir iwt en nor mm ion 

Kin run oi |h»i i mm mt ike I 
KOI) -rati ul renal hi*- m |sit i^-iuni 
/I U < fleet ot lonu tirv**urrouIWi vuitruuUir 

mpput 

1/1/0— rite ul M\v jell litthr itimi by murfh «vtU 


M02 — rate cf oxy gen utilisation by non mu-cle 
cell* 

BAE — total extracellular sod ttm 
BED — rate of change of sodium in intracellular 
fluids 

BID — rate of sodium Intake 
BOD— rate of renal excretion of sodium 
Oil M— mu*de oxygen utiiuation at rest 
0?v4 — aorde oxy gen e^lu ration 
OSK — non muvde venous oxytrsn -aturaiion 
OVA — oxygen volume in aortic blood 
OV-S — muscle venou* oxygen ^turation 
OZM — ba*ic oxygen utiliiation m non mU'cic 
body tissues 
PA — aortic pressure 

PAW— effect of arterial pressure in distending 
arteries ratio to normal 
PC — cutullary pressure 

PCD — net i /Tenure gradient across capillary 
membrane 

PCP — pulmonary capillary pressure 
PDO — difference between mu<rle venou> oxygen 
Pch and normal venous oxy «tn Poi 
PH — rate of transfer of fluid across pulmonary 
capillaries 

PFL— renal filtration pressure 
^PCC— coUoul osmotic pre-Mire of tt*<ur 
PGII — absorbency effect of gelvau-c*! bv rccoU of 
Ktl reticulum 

PCI. — procure gradient in lun--* 

P6P— colloid osmotic pressure of n—uCMilcuU-^d 
by eturapiwd jiroleln 

PGR — colloid osmotic itfexsurc of intcr-cuihil gel 
caused by Donnan euuihbnum 
PIF — inter* tit I'd fluid ine^sure 
PLA — left atrial jittvure 

PLD — pre»iure gradient to cause ly mphatic flow 
PZJ? — pulmonary ly mphatic flow 
PMO — muscle cell Pot 

POD— no tv-muscle venous Poi ramus normal v atue 
POK — sensitivity of rapid system of autoregula- 
tion 

PO.V — *etu4tl\ity of intermediate autoregulntion 

POS— iiulinonary interstitial fluid colloid osmotic 
pressure 

POT — non muscle cell Po» 

POK— nor. niu*de venou* Poi 
POY— ^en*itiv I ty of red cell production 
POZ— -wn-ativ ity of long term autoregulation 
P02 — oxygen deficit factor causing red cell pro 
duct Ion 

PPA— pulmonary arterial pre-imre 
PPC— pla*m i colloid o*mntie jire—urr 

PPD — HU of change of protein in pulmonary 
fluids 

PPI — pulmonary inter*Utml fluid pressure 
PPB — rate of pulmonary c tpilhiry protein lo*v 
PPO — pulmonary iy mph protein flow 
PPR— total protein m pnlmonary fluids 
/»RA — right atrnl pre**ure 
PRM — prc*xure causal by com predion of inter 
stitlal fluid uel reikniiuni 
PUP — toul plasma protein 
PI C — intcr<nual fluid colloid o-mouc pre-'iure 
PTS— oolid u-^ie pre*- ‘lire 
PTT — lord ti-**ue prc**ure 
l>Cy-|irrHire worn vein* to right atrium 
prfr-vriiiiii. pre* urc gndient 
PVfJ— mu^lr vvnou* Piv* 

?1 ^ --tv^TUKe venou* ptr^-ttre 

OAO —blood flow in tin <>*tenuc arterial system 
Ql V — lgi*ic Mt ventncul ir nutput 
(Jfflwiiitimiiif Irttvrnirsile ^ 


0011 — total voh me of oxygen m mu«cle cell* 
Q02- non muscle total celluhr oxygen 
QPO—r ate of blood flow into pulmonary vein* 
and left atrium 

QRF — feedback effect of left ventricular function 
on right ventricular function 
ORB — ba*!C right ventricular output 
ORO— actual right ventricular output 
QVO — rate of blood flow from veins into sight 
a mum 

RAM — ba«tc \ a*cular resistance of muscle* 
RAR— basic resistance of xton-muscular and non- 
renal irtcrici 
RBF — renal blood flow 
RCI— red cell production rate 
RC2 — red cell destruction rate 
RCD — rate o! change of red cell mas« 

REK — i*rcent of normal renal function 
RF V — renal blood flow if kidney i* not damaged 
RKC— rate factor for red cell destruction 
R WO — rate of oxygen transport to muscle cell* 
RPA — pulmonary Rtenal resistance 
RP7 — pulmonary vascular resistance 
RPV — pulmonary venous resistance 
RR — renal resistance 
RBM — vascular re<i-*tance in mu*cles 
R5Y — vascular resistance in non muscle noo- 
renal tissues 

RVG — re*l»Wnce from veins to right atrium 
RVXf — depressing effect on rig] it ventricle of 
pulmonary artenal pressure 
RV^— venous resistance 
SR — intenxity factor for sirexs relaxation 
^RIC— time constant for stress relaxation 
STB — effect of tissue hypoxia on salt and water 

intake 

SVO — stroke volume output 

TRR— tubular reabsorptlon rate 
TVD — rate of drinking 
VA 9— volume in systemic arteries 
VB — blootl volume 
VEC— extracellular fluid volume 
VG — volume of Interstitial fluid gel 
_VGD — rate of change of tissue gel volume* 

VtB — blood viscosity ratio to that of water 
VIC — celt volume 

VlD—tixz of fluid transfer between interstitial 
fluid and celt* 

VIE — portion of blood vivco-m caused by red 
blood cells 

VIF — volume of free inter>titul fluid 
V7Jf — blood vixco«it> (ratio io normal btood) 

VI A -volume m left atrium 

VP -plasma volume 

VP t — volume in pulmon-uv urteitcx 

t^PD-rate of change ol plx*in i volume 

VPF — pulmonary tree fluid volume 

VRA -rijjiit a' rial volume 

VRC- volume ol red blood cell* 

VI ( -rate of fluid ir tn*fer acro*« -ty stemic expd* 
lary membrane* 

VTD— rate of volume change m total interstitial 
fluid 

VTL —rate of *y itemk ly mph tow 
Vr^— total tntefititial fluid volume 
Vrit — total body tracer 
VL D— rate of tirin ixj output 
VV7— inertrused va*cular volume cansrtl by Mret* 
rtl ixatian 

I' I R — ih nu ui died v i>*cular volume nii<ed by 
*v miMtltetiC ■ttiuuilaiton 
MS -v cnou * van ill ir volume 
7t—imvevon-tant*ji autonomic re«pou*e 


Jfofewtl 

OP POOR 
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***** INPUT PARAMETERS **** 


A< 108 ) 

0*31 7222E+02 


AAR 

A ( 163 ) 

0 • 200000E+00 


AGK 

A ( 203 ) 

0.302806E+01 


AH 

A ( 204 ) 

0.100911E+01 


AHC 

A(247 ) 

0 * 700000E+01 


AHK 

A ( 208 ) 

0.100039E+01 


AHM 

A (288) 

0« 194279E-01 


AHY 

A ( 287 ) 

0#183110E-01 


AHZ 

A ( 205 ) 

0.792152E-01 


AH1 

A ( 206 ) 

0. 120009E+01 


AH2 

A ( 207 > 

0.499961E+01 


AH4 

A ( 201 ) 

1U689E-02 


AH7 

A (202) 

0 * OOOOOOE+OO 


AH8 

A ( 318 ) 

0.100000E+01 


AL0 

A( 161 ) 

0.993475E+00 


AM 

A{ 157) 

0 • 992485E+00 


AMC 

A ( 343 ) 

0 1 994666E + 00 


AMM 

A( 155 ) 

0. 101307E+01 


AMP 

A( 154) 

0* 980250E+00 


AMR 

A ( 238 ) 

0 ♦ 600000E+02 


AMT 

A ( 156 ) 

0.994222E+00 


AMI 

A ( 158 ) 

0. 168723E-01 


AM2 

A( 159) 

0, 103961E+01 


AM3 

A (160 ) 

0. 190455E+02 


AM5 

A( 166) 

0 • 995995E+00 


ANC 

A( 170) 

0. 100303E+01 


ANM 

A ( 164) 

0 • 995531 E+00 


ANP 

A ( 374 ) 

0 ♦ 995531 E+00 


ANR 

A ( 239 ) 

0» 150000E+02 


ANT 

A ( 269 ) 

0. 100303E+01 


ANU 

A ( 3 72 ) 

0.300000E-03 


ANV 

A ( 373 > 

0 1 000000E + 00 


ANW 

A { 37 1 ) 

-.446948E-01 


ANX 

A ( 369 ) 

* ♦ 200000E+00 


ANY 

A ( 370 ) 

0 • 400000E+00 


ANZ 

A ( 165) 

0 • 995826E+00 


A N 1 

A( 167) 

0.418338E-01 


AN2 

A( 168) 

0.110111E+01 


AN3 

A( 169) 

0 o 299697E + 01 


AN5 

A( 188) 

0.998304E+00 


A8M 

A( i 11 ) 

0.376883E+02 


APD 

A (329) 

0. 150000E+01 


ARF 

A (198) 

0.896017E+00 


ARM 

A( 194) 

0.985931E+00 


AR1 

A( 195) 

0.929921E+00 


AR2 

A( 197) 

0 • 977277E+00 


AR3 

A ( 61) 

0 » 988693E + 00 


AU 

A ( 54) 

0. 100375E + 01 


AU6 

A ( 53) 

0. OOOOOOE+OO 


AUC 

At 63) 

0.988693E+00 


AUH 

A ( 60) 

0« 988693E+00 


AUJ 

A ( 227 ) 

0.500000E-03 


AUK 

A ( 324 ) 

0 » 21 OOOOE+OO 


AUL 

A ( 66) 

0 • 990389E+00 


AUM 

A ( 55) 

O.OOOOOOEVOO 


AUN 

A ( 62) 

"• 1 13072E-01 


AU0 



A ( 281 ) 

0 • 988693E+00 


AUP 

A ( 375 ) 

0.100000E+01 


AUQ 

A( 376 ) 

0 « 988693E+00 


AUR 

A ( 377 > 

0* 100000E+01 


AUS 

A<282) 

0 ♦ 300000E+00 


AUV 

A { 226 ) 

0 • 300000E+01 


AUX 

A { 284 ) 

0 « 250000E+00 


AUY 

A ( 228 ) 

0 • 100000E+01 


AU2 

A< 57) 

-t969967E-02 


AU2 

A < 67) 

0.209558E-01 


AU4 

A ( 56 ] 

0 • 990300E+00 


AU6 

A { 58) 

-•484984E-05 


AU8 

A ( 65) 

0.840389E+00 


AU9 

A ( 38) 

0 • 997600E+00 


AVE 

A ( 26) 

0 • 101 126E+01 


A1B 

A ( 242 ) 

0. 100000E+01 


AIK 

A ( 243 ) 

0 • 200000E+02 


A2K 

A ( 244 ) 

0t 115200E+05 


A3K 

A ( 344 ) 

0. 100000E+01 


A4K 

A1332) 

0.994639E+00 


BFM 

A ( 1 86 ) 

0.295463E+01 


8FN 

A 1 1 74 ) 

0« 166661E+04 


Bi 

A 1 131 ) 

••400162E-02 


CCD 

A (230 ) 

0»700000E-02 


CFC 

At 92) 

0. 495104E + 01 


CHY 

A 1 122 ) 

0.498957E+01 


CKE 

A ( 129 ) 

0 • 1 42Q25E+03 


CKI 

A ( 130 ) 

0. 142029E+03 


CNA 

A ( 199 ) 

0.302961E+01 


CNB 

A< 162) 

0.997039E+01 


CNE 

A (245) 

0 » 1 39000E + 03 


CNR 

A ( 210 ) 

0.250000E+01 


CNX 

A (209) 

0 • 600000E + 01 


CNY 

A (246 ) 

0. 100000E + 01 


CNZ 

A ( 39) 

0.2) 2000E-01 


CN2 

A ( 40) 

0.366284E+00 


CN3 

A (225 ) 

0. 200000E+00 


CN7 

A ( 19) 

0 » 632025E+00 


CPA 

A (305 ) 

0.3000005-03 


CPF 

A ( 95 ) 

0. 124609E + 02 


CFG 

A ( 73) 

0. 165321L + 02 


CPI 

A (231 ) 

0. 160000E-06 


CPK 

A (291 ) 

0 « 30 1 352E + 02 


CPN 

A ( 75) 

0 * 70 1 082E + 02 


CPP 

A ( 233 ) 

0.850000E+02 


CPR 

A ( 86) 

0.535762E+02 


CPI 

A ( 224 > 

0 « 825000E"01 


cv 

A ( 49 ) 

-.125924E-01 


DAS 

A ( 59) 

0 * 988056E+00 


DAU 

A (307 ) 

0. 100000E+01 


DAI 

A (298 ) 

-.888210E-08 


DFP 

A ( 1 43 ) 

0.363798E-09 


DFZ 

A (273 ) 

0.554421E-02 


OHM 

A ( 50) 

0.128508E-03 


DLA 

A ( 90) 

0.700026E-02 

m 

DLP 

A (308 ) 

0.699915E-02 

m 

DLZ 

A ( 48) 

0.146389E-02 

m 

DPA 

A ( 87) 

0.532462E-01 

m 

OPC 

A ( 88 ) 

0.661090E-04 


DPI 
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A ( 85) 

0.531801E-01 


DPL 

A (835) 

0.700000E-02 


DP0 

A < 91) 

-.674377E-04 


DPP 

A ( 309 ) 

0.531878E-01 


DPY 

A ( 310 ) 

0.529841E-01 


DPZ 

A ( 1 87 ) 

0» 179681E+03 


D0B 

A ( 51 ) 

0.716209E-03 


DRA 

A ( 286 ) 

0.300000E+01 


DSP 

A ( 47) 

0.102838E-01 


DVS 

A ( 1 04 ) 

0.326489E+00 


EPH 

A ( 3 19 ) 

0* 100000E+01 


EXC 

A ( 328 ) 

O.OOOOOOE+OO 


EXE 

A (350 ) 

0.300000E+01 


EX1 

A (365) 

O.OOOOOOE+OO 


FIS 

A (368 ) 

O.OOOOOOE+OO 


GBL 

A (200) 

0. 125130E+00 


GFN 

A ( 1 14 ) 

0. 125130E+00 


GFR 

At 97) 

0.125131E+00 


GF1 

A ( 27 1 ) 

0.500000E-01 


GE2 

A ( 279 ) 

0* 100525E+01 


GF3 

A ( 280 ) 

0.500000E+01 


GF4 

A ( 1 12 ) 

0.620504E+02 


GUP 

A ( 107 ) 

-. 154896E-04 


GPD 

A ( 151 ) 

0.143459E+03 


GPR 

A ( 3 1 1 ) 

0. 187144E-04 


GPZ 

A ( 1 05 ) 

0.325108E-02 


GP1 

A ( 106 ) 

0.374288E-01 


GP2 

A ( 260 ) 

0.533.330E+00 


HKM 

A ( 173 ) 

0.408241E+02 


HM 

A (272 ) 

0. 100000E+01 


HMD 

A ( 259 ) 

0 . 900000E+02 


HMK 

A < 1 72 ) 

0.408841E+00 


HM1 

A (316 ) 

0. 100243E + 01 


HPL 

A (315 ) 

0. 100509E + 01 


HPR 

A (304 ) 

0*71 7308E+02 


HR 

A ( 2 19 ) 

0. 100000E+01 


HSL 

A (218 ) 

0. 100000E + 01 


HSR 

A ( 236 ) 

0.570000E+02 


HYL 

A ( 2) 

0 • 726600E+00 


I 

A ( 136 ) 

0.670675E+03 


11 

A ( 1 45 ) 

0.300000E-02 


12 

A ( 275 > 

0.200000E+02 


13 

A (150 ) 

0.914784E+01 


IFP 

A ( 1 27 ) 

0. 115573E-03 


KCD 

A (313 ) 

0. 119019E-03 


KCZ 

A ( 1 33 ) 

0 » 750493E+02 


KE 

A (128) 

-. 105172E-03 


KED 

A ( 124 ) 

0.698539E+03 


KE1 

A ( 1 34 ) 

0 o 354853E+04 


KI 

A ( 237 ) 

0.280000E-02 


KID 

A (126) 

0.915527E-02 


K IE 

A { 125 ) 

0.354854E+04 


KIR 

A (153 ) 

0.998025E+01 


KN1 

A ( 152 ) 

0.351305E-01 


KN3 

A ( 123 ) 

0.278960E-02 


K0D 

A ( 89) 

0. 148918E+02 


LPD 

A ( 234 ) 

0.470000E-03 


LPK 

A ( 12) 

0 » 99041 1 E + 00 


LVM 

A (340) 

0.598982E+02 1 


MM0 
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A<266 ) 

0. 179695E+03 


M62 

A ( 120 ) 

0.213630E+04 


NAE 

4(119) 

• • 3623772*02 


NED 

A (220 ) 

0. 100000E+00 


NID 

A ( 1 1 8 ) 

0.103624E+00 


N0D 

A ( 312 ) 

0. 103915E+00 


N0Z 

A ( 346 ) 

0 » 600000E+02 


0MM 

A ( 289 ) 

0.993753E+00 


0SA 

A( 190 ) 

0.695877E+00 


0SV 

A ( 285 ) 

0.300000E+01 


OUT 

A( 185) 

0 • 202845E+03 


OVA 

A (334) 

0.698845E+00 


0 VS 

A( 326 ) 

0. 150000E+00 


02A 

A( 320 ) 

0 * 1 80000E + 03 


02M 

A ( 10) 

0* 997387E+02 


PA 

At 11 ) 

Of 100262E+01 


PAM 

A ( 367 ) 

0.997387E+02 


PAR 

At 52) 

0.997387E+02 


PA1 

A ( 321 ) 

0* 100878E+03 


PA2 

A ( 78) 

0 • 1 837 55E+02 


PC 

A ( 79) 

0* 4541 61 E + 00 


PCD 

A { 232 ) 

0.300000E+01 


PCE 

A ( 306 ) 

0* 697848E+01 


PCP 

A ( 341 ) 

«t679865£«01 


PD0 

A ( 297 ) 

0f298392E-03 


PFI 

A C 1 13 ) 

Of 160071E+02 


PFL 

A ( 93 > 

0f61 1260E+01 


PGC 

A ( 100 ) 

• # 40001 OE+01 


PGH 

A ( 29) 

Of 152458E+02 


PGL 

A ( 96) 

Of 413219E+01 


PGP 

At 94 ) 

Of 198042E + 01 


PGR 

At 42) 

Of 959435E+02 


PGS 

A ( 35 ) 

0. 370362E+01 


PGV 

A ( 98) 

0.165288E+02 


PGX 

At 101 ) 

-.922036E-02 


PG2 

At 72) 

-.598890E+01 


PIF 

A ( 354 ) 

0 • 250000E+04 


PK1 

A ( 363 ) 

Of 800000E+03 


PK2 

A (364) 

0.200000E+01 


PK3 

At 23) 

Of 117865E+00 


PLA 

A ( 293 ) 

Of 298391E-03 


PLF 

At 81) 

0.800060E+00 


PL0 

At 25) 

0 f 20 1 1 79E+02 


PL1 

A t 30 1 ) 

0 • 686006E+01 


PMC 

A ( 338 ) 

0.800240E+01 


PM0 

A ( 303 ) 

0.461073E+01 


PMP 

A { 302 ) 

0 f 724995E + 01 


PMS 

A ( 347 ) 

0.800240E+01 


PM1 

A ( 348 > 

0. 100000E-02 


PMS 

A ( 349 ) 

100000E+01 


PM4 

A (353) 

0» 122000E+03 


PMS 

A (267) 

Of 929492E+00 


P0A 

A (193) 

Of 985898E+00 


P0B 

A (196) 

0.976478E+00 


P0C 

A ( 192 ) 

• • 237595E+00 


POD 

A (342) 

0 • 99456lE J -00 


POE 

A ( 240 ) 

0.600000E-01 


P0K 

A ( 345 ) 

0.800000E-01 


POM 

A ( 24 1 ) 

0 « 300000E+00 


PON 



A(274) 

0 • 800000E+01 


P90 

A (270) 

0 • 400000E+02 


P0R 

A (292) 

0.120541E+02 


P0S 

A ( 191 ) 

0.821768E+01 


POT 

A (261 ) 

0.397624E+02 


P0V 

A ( 268 > 

0 « 464000E"04 


P0Y 

A (262) 

0.300000E+00 


P0Z 

A ( 69) 

0.825000E+01 


P01 

A( 179) 

0.237500E+00 


P02 

A ( 17) 

0 * 153637E+02 


PPA 

A ( 76) 

0.280433E+02 


PPC 

A ( 296 ) 

0.159315E-05 


PPO 

A ( 290 ) 

100054E+02 


PPI 

A (295) 

0.899393E-02 


PPN 

A ( 294 ) 

Of 899208E-02 


PPO 

A (300) 

0.376523E+00 


PPR 

A ( 1 42 ) 

0.185170E-05 


PPZ 

A ( 18) 

0 • 399456E+00 


PPI 

A (322) 

0. 155392E+02 


PP2 

A ( 14) 

0.915550E-01 


PRA 

A ( 93) 

"•501114E+01 


PRM 

A( 149) 

0« 20771 lE+03 


PRP 

A ( 146 ) 

Of 915550E-01 


PR1 

A ( 74) 

0 • 41 3303E + 01 


PTC 

A ( 378 ) 

0 • 600000E + 02 


PTM 

A { 71) 

0.699994E+01 


PTS 

A ( 70) 

0.101104E+01 


PTT 

A ( 77) 

0. 145804E+02 


PVQ 

A ( 335) 

0 • 399320E+02 


PV0 

A ( 34) 

0 • 37951 7E+01 


PVS 

A( 189) 

0 • 800000E+01 


Pie 

A (339) 

0 • 800000E+01 


P20 

A ( 44) 

0.513746E+01 


OA0 

A ( 24) 

0 • 522090E+01 


QLN 

A ( 46) 

0.512487E+01 


QL0 

A ( 337 ) 

0.240003E+04 


O0h 

A ( 264 ) 

0.246777E+04 


Q02 

A ( 30) 

0.512500E+01 


QP6 

A ( 330 ) 

0.600000E+00 


ORF 

A ( 15) 

0« 521973E + 01 


QRN 

A { 45) 

0.512646E+01 


QR0 

A ( 37) 

0.512718E+01 


QV0 

A (327) 

0.394053E+01 


Q 1 

A ( 351 ) 

0. 180000E+02 


02 

A (352 ) 

0. 119632E+04 


Q3 

A ( 31 ) 

0.292328E+02 


Q5 

A (333 ) 

0.963000E+02 


RAM 

A (223) 

0 • 305200E+02 


RAR 

A (265 ) 

0.H8807E+01 


RBP 

A ( 182 ) 

-.834661E-06 


RCO 

A (181 ) 

0.110200E-04 


RC1 

A( 178) 

0. 118547E-04 


RC2 

A (263 ) 

Ot 554944E+03 


RD0 

A ( 1 17 ) 

0. 100000E+01 


REK 

A ( 1 10 ) 

0. 118807E+01 


RFN 

A ( 1 80 ) 

0# 580000E-05 


RKC 

A (336) 

Of 598937E+02 


RM0 

A ( 20) 

Of 158222E+01 


RPA 

A ( 28) 

Of 297457E+01 

m 

RPT 



A t 27) 

0 

139235E+01 

m 

RPV 

A( 109) 

0 

839500E+02 

m 

RR 

A ( 331 ) 

0 

964606E+02 

m 

RSM 

At 184 ) 

0 

324723E+02 

m 

RSN 

At A3) 

0 

193962E+02 

m 

RTP 

At 36 ) 

0 

721443E+00 

m 

RVG 

At 21) 

0 

989118E+00 

m 

RVM 

At 41 ) 

0 

275684E+01 

m 

RVS 

A ( 283 ) 

0 

273012E+01 

m 

RVl 

A (221 ) 

0 

500000E+00 

m 

SR 

A ( 248 ) 

0 

330000E+02 

m 

SRK 

A ( 366 ) 

0 

000000E+00 

m 

STA 

A ( 317 ) 

0 

100000E+01 

m 

STH 

A ( 323 ) 

0 

714458E-01 

m 

SV0 

A < 

1 ) 

0 

OOOOOOE+OO 

m 

T 

A ( 249 ) 

0 

191600E+05 

m 

TM 

At 115) 

0 

124098E+00 

m 

TRR 

A4 216 ) 

0 

100536E-02 

m 

TVO 

At 141 ) 

0 

100391E-02 

m 

TVZ 

A ( 278 ) 

0 

OOOOOOE+OO 

m 

T 1 

A (276) 

0 

400000E+01 

m 

U 

At 

9) 

0 

354072E+00 

m 

VAE 

At 

6 ) 

0 

849072E+00 

m 

VAS 

A ( 171 ) 

0 

500662E+01 

m 

VB 

At 

3) 

a* 

392832E-05 

9 

VBD 

At 121 ) 

0 

150413E+02 

m 

VEC 

At 103) 

0 

115125E+02 

m 

VG 

At 176) 

0 

303328E+01 

m 

VIB 

At 135 ) 

0 

249852E+02 

m 

VIC 

At 132 ) 

m 

120985E-04 

m 

VID 

A ( 1 75 ) 

0 

153328E+01 

m 

VIE 

At 68 ) 

0 

553337E+00 

m 

VIF 

At 177) 

0 

101099E+01 

m 

VIM 

A (314 ) 

m 

400162E-04 

m 

VIZ 

A 1 102 ) 

D 

184407E-03 

m 

VG0 

At 

7) 

0 

401 1 79E+00 

m 

VLA 

At 22) 

0 

117865E-02 

m 

VLE 

A ( 148 ) 

0 

296271E+01 

m 

VP 

A (277 ) 

0 

100000E-01 

m 

VP1 

At 

5) 

0 

379996E+00 

m 

VPA 

At 84) 

m 

205207E-04 

m 

VPD 

At 16) 

0 

737455E-01 

m 

VPE 

A ( 299 ) 

0 

124944E-01 

m 

VPF 

At 

8) 

0 

100458E+00 

m 

VRA 

A ( 183 ) 

0 

204391 E+01 

m 

VRC 

At 13) 

0 

457775E-03 

m 

VRE 

At 80) 

0 

319994E-02 

m 

VTC 

At 83) 

0 

124043E-06 

m 

VTD 

At 82) 

0 

320494E-02 

m 

VTL 

At 147) 

0 

120661E+02 

m 

VTS 

A t 217 ) 

0 

400264E+02 

m 

VTW 

At 137 ) 

0 

317912E-02 

m 

VTY 

At 139 ) 

0 

320024E-02 

m 

VTZ 

At 116) 

0 

103115E-02 

m 

VUD 

At 140) 

0 

103220E-02 

m 

vuz 

At 32) 

0 

325152E+00 

m 

WE 

At 

4) 

0 

327592E+01 

m 

VVS 

A (222) 

0 

2951 37E+01 

m 

VVR 

A ( 2 1 1 ) 

0 

162576E+00 

m 

VV1 
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A { 212 ) 

0. 150500E+0Q 


W2 

A ( 64) 

0 • 207625E+00 


VV4 

A ( 213 ) 

0.253618E-04 


VV5 

A ( 214 ) 

0.827931E-05 


VV6 

A (215 ) 

0. 120510E-01 


VV7 

A ( 33) 

0 • 3131 02E+00 


VV8 

A ( 325 ) 

0.315900E+01 


VV9 

A (250 > 

0.200000E-01 


V2D 

A ( 144 ) 

0 • 100000E+02 


X 

A ( 229 ) 

0. 100000E+01 


Y 

A ( 138 ) 

0* 100000E+01 


Z 

A ( 251 ) 

0.100000E+01 


Zl 

A ( 252 ) 

0. 100000E+01 


Z2 

A ( 253 ) 

0 * 400000E+Q 1 


Z3 

A ( 254 ) 

0.100000E+02 


Z4 

A (255) 

0.100000E+02 


Z5 

A ( 256 ) 

0.500000E+01 


Z6 

A ( 257 ) 

0 » 500000E+01 


Z7 

A (258) 

0. 100000E+01 


Z8 

A( 355 ) 

0 « 120000E+00 


Z9 

A (356) 

0.825000E+01 


Z10 

A ( 357 ) 

0 • 400000E+01 


Zll 

A (358) 

0.124000E+01 


Z12 

A (359) 

0 ♦ 625000E+00 


Z13 

A (360) 

0 • 000000E+00 


Z14 

A ( 361 ) 

O.OOOOOOE+OO 


Z 15 

A (362) 

0 • OOOOOOE+OO 


Z16 

A ( 379 ) 

O.OOOOOOE+OO 



A ( 380 ) 

O.OOOOOOE+OO 



A ( 381 ) 

O.OOOOOOE+OO 



A ( 382 ) 

O.OOOOOOE+OO 



A ( 383 ) 

O.OOOOOOE+OO 



A ( 384 ) 

O.OOOOOOE+OO 



A (3851 

O.OOOOOOE+OO 



A ( 386 ) 

O.OOOOOOE+OO 



A ( 387 ) 

O.OOOOOOE+OO 



A( 388) 

O.OOOOOOE+OO 



A ( 389 ) 

O.OOOOOOE+OO 



A (330) 

O.OOOOOOE+OO 



A ( 391 ) 

O.OOOOOOE+OO 



A ( 392 ) 

O.OOOOOOE+OO 



A ( 393 ) 

O.OOOOOOE+OO 



A (394) 

O.OOOOOOE+OO 



A ( 395 ) 

O.OOOOOOE+OO 



A (396) 

O.OOOOOOE+OO 



A (397) 

O.OOOOOOE+OO 



A (398) 

O.OOOOOOE+OO 



A (399) 

O.OOOOOOE+OO 



A ( 400 ) 

O.OOOOOOE+OO 
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APPENDIX C 



17 


1 C 

2 C 

3 C 

4 C 

5 C 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 C 

58 


R67093# 8 LARRY NEAL — HATH; NASA FILE 


PROGRAM GUYTON 

CIRCULATORY DYNAMICS - CIRCE 
CIRCE 1 

REAL LVM#I# IFP#LPD# KE# KE1 #KOD#K IR# K IE# K I# KCD# KED# KN1# KN3 
REAL NA£#NEO# NID-NOD* I 1# LPK#K ID#M02# N0Z# KCZ#HPL#HPR# 12# I3#MM0 
DIMENSION f-UNl ( 14 ) # FUN2 ( 1 4 ) # FUN3 { 1'4 ) # FUN4 ( 1 4 ) # FUN6 < 14 )#FUN7< 14) 
COMM0N/ARRAY/T# 1/ V8D# VVS; VP A/ VAS# VLA# VRA# VAE# PA# P AM# LVM# 

* VRE# PRA# QRN# VPE# PPA# PP1# CPA# RP A# RVM# VLE# PLA# QLN# PL 1# 

¥ A1B#RPV#RPT#PGL#QP0#Q5 # WE# VV8# PVS# PGV# RVG# QVO# AVE 

COMMON/ ARRAY/CN2# CN3# RVS# PGS# RTP# QA0# GR0# QLO# DVS# DP A# DAS# DLA# ORA# 
¥ PA1#AUC# AU8# AUN# AU6# AU2# AU8#DAU# AUJ# AU # AU0# AUH# VV4# 

¥ AU9# AUM# AU4# VIF#P01#PTT#PTS#PIF#CPI#PTC#CPP#PPC#PVG 

C0MMON/ARRAY/PC # PCD# VTC# PLD# VTL# VTD# VPD# DPL# CPI# DPC# DPI#LPD# DLP# 

* DPP,CHY#PRM/PGR#CPG#PGP#GFl# p GX#PGC#PGH#PG2#VGD#VG # 

* EPH# GP1# GP2# GPD# AAR# RR # RFN# APD# GL p # PFL# GFR# TRR# VUD 
C0MM0N/ARRAY/REK#N0D#NED#NAE#VEC# CKE#K0D#KE1#KIR#KIE#KCD#KED# CKI# 

* CNA# CCD# VI D# KE #KI #VIC#Il #VTY#Z # VTZ# VU2# TVZ# PPZ# 

¥ DFZ# X #12 # PR1# VTS# VP # PRP# IFP# GPR# KN3# KN1# AMR# AMP 

COMMON/ ARRAY/ AMI# AMC# AM2# AM3# AM5# AM # CNE# AGK# ANP# AN1# ANC# AN2# AN3# 

* AN5# ANM# VB # HM 1 # HM #B1 # VI E# VI B# VIM# RC2# P02# RKC# RC1# 

* RCD# VRC# RSN# OVA# BFN# DOB# AOM# P10# OSV# POT# POD# P08# AR1 
COMMON/ ARRAY/ AR2# POC# AR3# ARM# CNB# GFN# AH7# AH8# AH # AHC# AHl# AH2# AH4# 

* AHM# CNY# CNX# VV 1# VV2# VV5# VV6# VV7# TVD# VTW# HSR# HSL# NID# 

¥ SR # VVR# RAR# CV # CN7# AUX# AUK# AUZ# Y # CFC# CPK# PCE# CPR 

COMMON/ ARRAY/LPK# DPO# HYL# KI D# AMT# ANT# POK# PON# AIK# A2K# ASK# CNR# CNZ# 
¥ AHK# SRK# TM # V2D# Z 1 # Z2 # Z3 #Z4 # Z5 # Z6 #Z7 #Z8 # HMK# 

¥ HKM# P0V# POZ# RDO# Q02# RBF# M02# POA# POY# ANU# POR# GF2# HMD 

COMMON/ ARRAY/DHM#POQ# 13 #U #VP1#T1 # GF3# GF4# AUP# AUV# RV1# AUY# OUT# 
¥ DSP# AHZ# ANY# OS A# PP I# CPN# P0S# PLF# PP0# PPN# PPD# PF I # OFP# 

¥ VPF# PPR# PMC# PMS# PMP# HR # CPF# PCP# DA 1# OLZ# DP Y# DPZ# GPZ 

COMM0N/ARRA Y/NOZ# KCZ# VI Z#HPR# HPL# STH# ALO#EXC# 02M# PA2# PP2# SVO# AUL# 
¥ VV9# 02 A# Ql #EXE# ARF#QRF#RSM#BFM#RAM#SVS#PVO#RMO#Q0M# 

¥ PMO# P20# MM0# PD0# POE# AMM# A4K# P0M# 0MM# PM1# PM3# PM 4# EX1 

COMM0N/ARRAY/Q2 #Q3 #PM5#PK1#Z9 # Zl 0# Z1 1# Z12# Z13# Zl 4# Z1 5# Z1 6# PK2# 
¥ PK3# FIS# STA# PAR# SBL# ANY# ANZ# ANX# ANV# ANW# ANR# AUQ# AUR# 


¥ AUS# PTM #DUMMY<22)#TITLE(400) 

DATA FUNK1 ) # FUNI ( 2 ) # FUNI { 3 ) # FUNI ( 4 > # FUNI ( 5 ) # FUNK 6 ) # FUNi ( 7 ># 
¥FUN1 ( 8 ) # FUN1 ( 9 ) # FUN1 ( 10 ) # FUN1 f 11 >#FUN1 ( 12)# FUNI I13)#FUN1(14)/ 
¥0»# 1» 04# 60. #1*025# 125.# .97# 160.# .88# 200 •# *59# 240. # 0 •# 240. # 0 •/ 
DATA FUN2 ( 1 ) # FUN2 ( 2 ) # FUN? < 3 > # FUN2 < 4 ) # FUN2 ( 5 ) # FUN2 ( 6 ) # FUN2 ( 7 > # 
¥FUN2(8)#FUN2(9)#FUN2< 10 ) # FUN2 ( 1 1 ) #FUN2 l 12 > # FUN2 ( 13 ) #FUN2 ( 14 ) / 
¥-100.#0»0#«6.#0.0#-3.#»75#-l.#2i6#2»#9.8#8.#13.5# 1000*# 13.5/ 
DATA FUN3 ( 1 ) # FUN3 ( 2 > # FUN3 i 3 > # FUN3 t 4 ) # FUN3 { 5 ) # FUNS ( 6 ) # FUNS ( 7 > # 
¥FUN3(8)#FUN3( 9)#FUN3( 10 ) # FUN3 ( 1 1 ) #FUN9 < 12 ) # FUN3 113 ) #FUN3{ 14 )/ 
¥0*0# 1*06# 20* # *97# 24.# .93#30»# «8#38 . # . 46#45. # 0 •# 45 •# 0. / 

DATA FUN4( 1 ) # FUN4 < 2 )# FUN4 < 3 ) # FUN4 < 4 )# FUN4 < 5 )# FUN4 < 6 )#FUN4(7 )# 
¥FUN4( 8)#FUN4( 9)#FUN4( 10 >#FUN4( 11) #FUN4< 12)#FUN4( 13)#FUN4< 14 >/ 
¥-100.#0.#-4.#0.#-l.#3«6#3»# 9.4#6«#ll«6#10»#13»5#1000»#13.5/ 
DATA FUN6 { 1 ) # FUN6 ( 2 > # FUN6 < 3 > # FUN6 < 4 ) # FUN6 < 5 )# FUN6 < 6 > # FUN6 ( 7 )# 
¥FUN6 < 8 ) # FUN6 ( 9 ) # FUN6 ( 10 ) # FUN6 1 1 1 ) # FUN6 ( 12 ) # FUN6 ( 13 ) #FUN6 ( 14 >/ 
¥-100.#10000.#0.#70.# »4#9.3#.8#3»3#1»2#1»3#1»6# .43# 100»#0./ 
DATA FUN7 ( 1 ) # FUN7 < 2 ) # FUN7 ( 3 ) # FUN7 ( 4 > # FUN7 ( 5 ) # FUN7 ( 6 > # FUN7 ( 7 > # 
¥FUN7(8)#FUN7( 9) #FUN7( 10)#FUN7( tl)#FUN7( 12 > # FUN7 < 13 ) #FUN7 < 14)/ 
¥0.#7.#30t#6.25#60»#3.#100.#l»#160.# *15# 400.# »05#400»# *05/ 


WRITE ( 102# 5 ) 



i8 


59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 
117 


WRITE. ( 6/ 5 ) 

5 FORMAT </♦ GUYTON MODEL FROM WHITE'/ 

* ' REFER TO GE-AGS USER GUIDE TIR 741-MED-3017 ' // ) 
WRITE ( 102/60) 

60 FORMAT) 10X/ 'KEY-IN CODES' /10X/ '—— — — —• // 
*5X/'001 « INITIALIZE FROM CARDS'/ 

*5X/'002 * CHANGE VARIABLES'/ 

*5X/'0O3 ■ PRINT OUT VARIABLES'/ 

*5X/'004 » PRINT OUT COMPLETE ARRAY'//) 


90 CALL INPUT 
C* 


PRTIME=PTM 


IF C I .GT. 0.5) 1 = 0.5 
IOUT-OOT 

100 IFIT .LT. PRTIME) GO TO 300 
CALL INPUT 
PRTIME- PR TIME+PTM 


C 

C100 IF < OUT .EQ. 3* ) 
C IF(DSP .EG. 3.) 

C 

300 T = T + 1 2 
C 

CALL MEMO 

♦ 

•* 

* 

* 

* 

* 

* 

C 

120 CALL AUTO 
♦ 

* 

* 

C 


CALL OUTPUT 
CALL DSPLAY 


( AMM/ ANM/ ANU/ ANY/ ANZ/ ARM/ AUH/ AUM/ AUY/ AVE/ BFM/ BFN/ 
CN2/ CN3/ CN7/ CV / DAS/ OLA/ DP A/ DRA/DVS/ FIS/ HMD/ HPL/ 
HPR/HSL/HSR/I2 /LVM/PA /PAM/PA2/PC / PGL/ PGS/ 3 LA/ 
PPA/ PP1/ PP2/ PRA/ PR1/ PVS/ OAO/ QLN/ QLO/ QPO/ QRF/ QRN/ 
QRO/ QVO/ RAM/ RAR/ RBF/ RPA/ RPT/ R PV/ RSM/ RSN/ RVG/ RVM/ 
RVS/U /VAE/ VAS/VBD/VIM/ VLA/VLE/VP / VP A/ VPE/ VRA/ 
VRC/ VRE / WE/ VVR/ VVS/ VV7/ VV8/X / FUN1/ FUN2/ FUN3/ 
FUN4 ) 

( AU / AUB/ AUC/ AUH/ AU J/ AUK/ AUL/ AUM/ AUN/ AU0/ AUP/ AUQ/ 
AUR/ AUS/ AU V/ AUX/ AUZ/ AU4/ AU6/ AU8/ AlB/ DAU/ EXC/ EXE/ 
EX1/I2 /PA /PA1/P0Q/POT/P2O/STA/VVR/VV9/Y / Z/ 

Z8 / Z12 ) 


C 


C 


c 


c 


I F I I3«LE»I2)G0 TO 168 
IF ( ABS < DAU-AU J ) »GT . DAI ) GO TO 100 
110 IF ( ABS ( QAO-QLO ) »GT * * 2 ) GO TO 100 
IF ( ABS ( QAO-QPO ) .GT . . 2 ) GO TO 100 
IF (ABS(QAO-QRO) .GT».4)G0 TO 100 


168 CALL HORM0N 
* 

« 


(AM /AMC/AMP/AMR/AMT/AM1/ANM/CKE/PA/Z/FUN7/ 
AGK/ANC/ANP/ANR/ANT/ANV/ANW/AN1/CNA/CNE/GFN/ 
I /REK) 


170 CALL BLOOD 
♦ 


(HKM/HM /HMK/I /P0T/POY/PO1/P02/RC1/RC2/RCD/RKC/ 
VB / VIB/ VIE/VIM/VP /VRC) 


180 CALL MUSCLE 
* 

* 


(AL0/AMM/A0M/AUP/A4K/BFM/EXC/HM /I /MMO/OMM/OSA/ 
0VA/0VS/O2A/PD0/PK1/PK2/PK3/PM0/PM1/PM3/PM4/PM5/ 
POE/P OM/PVO/P20/QOM/RM0/VPF/Z5 /Z6) 


CALL 


AUTORG 


( AOM/ ARM/AR1/AR2/ AR3/A1K/ A2K/A3K/BFN/D0B/HM >1/ 
MO2/OSV/OVA/O2M/POA/P0B/P0C/POO/P0K/P0N/P0R/P0T/ 
P0V/POZ/P10/Q02/RD0/Z / Z4 zZ7> 


C 


* 

* 
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118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 


call 

4 

ADH 

(AH /AHC/AHK/AHM/AHY/AHZ* AH7/ AH8/AUP,*CNA/CNB/CNR/ 
CNZ/I /PRA/Z) 

CALL 

4 

MISC1 

(AHM/ AU4/AU8/ I /SR / SRK/ STH/ TVD/ T VZ# VEC/ V ICz V TW/ 
WE/ VV6/VV7/Z ) 

CALL 

4 

HEART 

(AUR/DHM/HMDjHR /I iPA /PMC/ PMP/PMS/ P0T/PRA/ QA0/ 
QL0/RTP/SV0/VAE/VLE/VPE/VRE/VVE) 

130 CALL 

CAPMBD 

(BFN/CFC/CPI/CPP/DFP/I /IFP/PC /PCD/PIF/PLO/PPC/ 


PRP/PTC/PTS/PTT/PVG/PVS/RVS/TVD/VG /VID/VIF/VP/ 
VPD/VTC/VTD>VTL/VTS/VUD/Z *Z1 / FUN6 ) 


I«I*1»2+T-T1 
I1=ABS< VP1/VPD/I ) 

IF(I1.LT. I) I = H 

IF ( 13 + T“T1 » LT * I ) I = I3 + T-T1 

T=I+T1 

T1=T 

IF (T.LT.PRT1ME) G0 TO 200 
CALL INPUT 
PRTIME-PRTIME+PTM 


C 

C 

C 

C 

C 

C 


C 

C 


I F ( 8UT . EQ . 4 • ) 

I F ( OUT . LQ . 4 . ) 
IF(DSP.EQ.4. ) 

200 CALL PULM0N 

♦ 

CALL MISC 2 
135 CALL PRflTtN 

4 

4 

142 CALL MONEY 

4 

* 


160 CALL IONS 
♦ 


CALL PUT0UT 

CALL OUTPUT 
CALL DSPLAY 

(CPF/ CPP/CPN/DFPz I /PCP/PFI/PLA/PLF/POS/PPA/PPC/ 
PPD/PPI/PPN/PP0/PPR/VP /VPD/VPF/Z >Z3) 

(HPL/HPR/HSL# HSR/ I/PA/PPA/P0T/STH/Z1O/Z11/Z13) 

(CHY/CPG/CPI/CPK/CPP/CPR/CPl/DLP/DLZ/DPC/DPIz dpl/ 
DPP/DPY/GPD*GPR/ I / IFP/LPK/PC / PCE* PGX/ PRP/ VG / 
VTL/Z / DP0/ PPD ) 

(AAR/ AHM/ AM / APD/ ARF/ AUMz CNE/CNX/ CNY/ GBL/ GFNz GFR/ 
GF2/GF3/GF4/GLP/ I / NAE/ NEO/N ID/ N0D/ N0Z/ P A /PAR/ 
PFL/PPC/RBF/REK/RFN/RR /STH/TRR/VIH/ VUO/Z ) 

(AM /CCD/CKE/CKI/CNA/I /KCD/KE /KED/KI /KlD/KIE/ 
KIR/K0D/NAE/REK/VEC/VIC/VIO/VP / VPF/VTS/Z ) 


C 


C 


140 CALL GELFLD 

* 


(CHY/CPG/CPI/GPR/HYL/ IFP/PGC/PGH/PGP/PGR/PGX/PIF/ 
VTS/PRM/PTC/PTS/PTT/VG /VG0>VIF/VRS/V20/FUN6) 


G0 T0 100 
END 

SUBROUTINE HEM0 

4 
4 
4 
4 
4 
4 
4 

real 


( AMM/ ANM/ ANU/ ANY/ ANZ/ ARM/ AUH/ AUM/ AU Y/ AVE/ BFM/ BFN/ 
CN2/ CN3/ CN7/ CV / DAS/ DL A/ DPA/ DR A/ DVS/ FIS/ HMD/ HPL/ 
HPR/HSL/HSR/I2 /LVM/PA /PAM/PA2/PC / PGL/ PGS/ »LA/ 
PP A/ PP1/ PP2/ PR A/ PR 1/ PVS/ QA0/ QLN/ OL0/ QP0/ QRF/ QRN/ 
QR0/ QV0/ RAM/ RAR/ RBF/ RPA/ RPT/ RPV/ RSM/ RSN/ RVG/ RVM/ 
RVS/U /VAE/VAS/VBD/ VIM/ VLA/VLE/ VP / VPA/ VPE/ VRA/ 
VRC/VRE/VVE/VVR/VVS/ VV7/ VV8/X /FUN1/FUN2/FUN3/ 
FUN4) 


I2/LVM 
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177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 


C CIRCULATORY DYNAMICS BLOCK 
C HEMODYNAMICS 
C 

VBD“VP+VRC-WS-VAS-VLA-VPA-VRA 

VVS«VVS+DVS*I2+VBD*«3986 

VPA*VPA+DPA*I2+VBD*»155 

VAS=VAS+DAS*I2+VBD*»261 

VLA=VLA*DLA*I2+VBD*.128 

VRA«=VRA + DRA*I2+VBD*«0574 

VAE“V AS- ♦ 495 

PA-VAE/. 00355 

P AM= 1 00 * /PA 

PA2=PA/AUH 

CALL FUNCTN(PA2*LVM>FUN1) 

VRE“VRA-.l 

PRA-VRE/.005 

CALL FUNCTN IPRA>QRN*FUN2 ) 

VPE«=VPA». 30625 
PPA*VPE/.0048 
PP1=»026*PPA 
IF (PP1.LT.0. )PPl-0. 

RPA=PPHMM-.5 > 

PP2*PPA/AUH 

CALL FUNCTN (PP2jRVM/FUN3) 

VLE=VLA-.4 

PLA=VLE/.01 

CALL FUNCTN(PLA,QLN/FJN4) 

RPV=1 ./(PLA + 20. ) / • 035 7 

rpt*rpv+rpa 

PGL»PPA-PLA 

OP0»PGL/RPT 

anu-anm 

IF ( ANU*LT« *8 )ANU*.8 
VVE«VVS-VVR-( ANU-1 • )*ANY 
VV8»VVE-V V7 

I F C VV8.LT. .0001 )VV8=. 0001 

PVS«V V8/CV 

PR1*PRA 

IF (PRA.LT.O. )PR1*0. 

RVG«2»738/PVS 
QV0«(PVS-PR1 J/RVG 

CN3«CN3+( ( (PC-17. )*CN7+17. ) *CN2-CN3 )*. 1 
AVE- ( AUM-1 . ) ♦ AUY+ 1 • 

RVS*AVE*( 1* /CN3 )*VIM* ( ( ANU-1 • )*ANZ+1 • > 

PGS-PA-PVS 

RSN=RAR*ARM*ANU*AUM*PAM*VIM+RVS*1.79 

bfn«pgs/rsn 

RSMs>ANU*VIM*PAM*AUM*AMM*RAM 

BFM=PGS/RSM 

QA0=BFN+BFM+RBF+( pa-pra )*FIS 
QLe*LVM*QLN+AUH*HSL*HMD*HPL 

QRQsQRN* ( ( 1 .-QRF ) *AUH*RVM*HSR*HMO*HPR+GRF»QL0/GLN ) 

OP0=QLO+ ( QP0-QL0 ) /U 

QV0»QR0+( QV0-QR0J/X 

DVS*QA0-QVO 

DPA-OR0-QP0 

DAS-QL0-QA0 

DLA*QP0«QL0 

DRA»Q V0-QR0 
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RETURN 

END 

SUBROUTINE AUTO ( AU / AUB/ AUC/ AUH/ AU J/ AUK/ AUL/ AUM/ AUN/ AU0/ AUP/ AUQ/ 

* AUR/ AUS, AUV/ AUX/ AUZ/ AU4/ AU6/ AU8/ AlB/ DAU/ EXC/ EXE/ 

* EX1/I2 /PA /PA1/P0O/P0T/P2Q/STA/VVR/VV9/Y / Z/ 

* Z8 / Z12 ) 

REAL 12 

C 

C AUT0N0MIC C0NTR0L BL0CK 
C 


120 EXE*( 8.-P20 )*EX1+(EXC-1 • )*Z12 
POQ — P 0 T 

IF (P0Q.GT.8» )P0Qos. 

IF (P0Q.LT. 4. )P0Q»4. 

PAl»PA*P0Q/8.-EXE 

AUOO. 

IF ( Pa 1 • LT ♦ 80 . )AUC*.03*< 80.-PA1 ) 

I F ( PA1 • LT » 40. ) AUC = 1 * 2 
AUB-0 * 

IF(PA1.LT.170«) AUB* « 014286* ( 170 «*PA1 ) 

IFIPA1.LT.40. ) AU8“1 *83 

123 A1B*< AU8-1 . I *AUX+1 . 

124 AUN=0 

IF(Pa 1.LT.50. >AUN=.2*(50.-PAl ) 

IF ( PA 1 . LT * 20 . ) AUN = 6.0 
AU6-A1B-AU4 
AU8=AUK*< AU6-1. ) 

DAU=0AU+( AUC+AU6+AUN-DAU ) /Z/Y 
AUJ*AUJ+ ( OAU-AU J ) * 1 2*6 « /Z8 
IF ( AUJ.LT.O. ) AUJ-0. 

IF ( AUJ-1 . ) 126/127/127 

126 AU»AUJ**AUZ 
G0 TO 128 

127 AU* ( AUJ-1 . ) *AUZ+1 * 

128 IF(STA.GT..0OOOl)AU=STA 
AU0*AU-1 . 

AUP«AU0*AUO+1 » 

AUH«*AU0*AUV + 1» 

AUR*AU0*AUS+1 • 

VVR*VV9-AUL*AUP 

AUMs. 15+. 854AUP 

RETURN 

END 

SUBROUTINE H0RM0N { AM / AMC/ AMp, AMR/ AMT/ AMI / ANM/ CKE/ P A/ Z/ FUN7/ 

¥ AGK/ ANC/ANP/ ANR/ ANT/ ANV/ ANW/AN1/CNA/CNE/GFN/ 

* I /REK) 

REAL I 

c 


aldosterone control block 


C************ ********>Hf.3f*^»t-****iHf****‘Hf*>f.>HHHf***if***^***if******4-4.*](.*^^4.^4.^ 
168 AMR=CKE/CNA/. 00352-9. 

IF ( AMR.LT.O. ) AMR = 0 » 

CALL FUNCTN < PA/ AMP/FUN7J 
AM1saH1+( ANM* AMP* AHR-AM1 ) /Z 
AMC=AMC+(AM1-AMC)*( 1 ♦ -EXP | - I/AMT ) ) 

AM-20. 039-1 9. 8*EXP(-. 0391 *AMC> 
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C************* ¥************** **¥*******¥**¥**********¥****¥**¥* ********** 

C 

C ANGIOTENSIN CONTROL BLOCK 
C 

C************** ************************************************ ********** 
CNE-152.-CNA 
IF ( CNE»LT • 1 * )CNE«1» 

ANRal < 17»75-GFN*CNA >*AGK+1 « >*REK 
ANW«ANW+( ( anr-i. )*10.-ANW )*ANV*I 
IF( ANw.lt . .0) ANW*0. 
anp*anr+anw 

I F ( ANP * GT • 100 » ) ANPo 100 • 

IF(ANP«LT * .01 ) ANP- « 01 
ANl=ANl+{ ANP-AN1 ) /Z 
ANOANC-M AN1-ANC>*( l.-EXP (-I/ANT) J 
ANM=*4 »0-3«3*EXP(». 0967*ANC ) 

IF( ANM.LT. .7 ) ANM*> • 7 

RETURN 

END 

SUBROUTINE BLOOD (HKMjHM ,HMK*I * POT/ P6Y. P01 / P02/RC1 » RC2* RCD * RKC* 

* VB /VIB, VIEjVIM/VP / VRC ) 

REAL I 

C 

C RED CELLS AND VISCOSITY BLOCK 

C — — — --- • ...... 

C BLOOD VISCOSITY 

C — --- — — «»—••—•,« — 

170 VB-VP+VRC 

HM=100.*VRC/VB 
VIE=HM/(HMK»HM)/HKM 
VIB“VIE4-1 .5 
V IM* . 3333*V IB 

O---— .......... ........................................... 

C RED BLOOD CELLS 

C — • . 

RC2=RKC*VRC 

P02-P01-P0T 

IF < P82.LT. .2375) P02=. 2375 

RC1=P0Y*P02 

RCD=RC1-RC2 

VRC®VRC+RCD*I 

RETURN 

END 

SUBROUTINE MUSCLEt ALS^ AMM,A0 M*AUPj A4Kj BFM, EXCjHM *1 >MM0>0MM.0SA* 

* 0VA>0VSiO2A#PDO>PKl < iPK2jPK3>PM0/PMl>PM3/PM4^ PM5> 

* P0Ej P0M.PV0* P20jQ0MjRM6j VPF*Z5 /Z6» - 
REAL I t MMQ 

C 

C MUSCLE BLOOD FLOW C0NTROL AND P02 BLOCK 
C 

180 0SA»AL0«VPF*.5 
0VA*0SA*HM*5. 

0VS=0VS+( 1BFM*0VA-RM0 )/HM/5»/BFM*0VS) /Z6 
P V0 = 5 7 . 14*0 VS 

RM0n( PV0-PM0) *PM5/( PM1**PK3»PM4 } 

Q0M«=Q0M+<RM0-MM0)*< 1 . -EXP < -I/Z5 > ) 

PM0«PK2/( PK1-Q0M) 

PM1=PM0 

IF ( PM1«LT *PM3)PM1*PM3 
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373 
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379 

380 

381 
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383 
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389 
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391 
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393 
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P20=PM0 

1F(P20.6T*8.)P20«8. 

A0M*( AUP-1. U02A+1. 

MM0*A0M+0MM+EXC*( 1 » - ( 8 . 0001 -P20 ) **3 » /51 2* > 

PD0=PV0-4O. 

P0E=P0M*PD0+1 . 

IF (P0E.LT. «OO5)P0E*»OO5 
AMM«AMM+{ P0E-AMM ) *( 1 . -EXP ( -I/A4K ) ) 

RETURN 

END 

SUBROUTINE AUT0RG ( A0M* ARM, AR1 1 AR2, AR3, A IKj A2K., A3IO BFN* D0B> HM »lt 

* M02-»0SVj0VA^02M#P0A^P0B/P0CjP0D/P0K#P0N#P0R#P0T/ 

* P0VjP0Z#P10j»O 02^ RD0>2 iZ4 >Z7) 

REAL I * M02 

C 

C NON-MUSCLE 0XYGEN DELIVERY BLOCK 

C AND NON-MUSCLE LOCAL BL00D FLOW CONTROL BLOCK 

C AUTOREGULATIONiRAPID 

c— ........................ ....................................... 

0SV*0SV+( ( BFN+0VA-D0B )/HM/5»/BFN-0SV)/Z7 

P0V=0SV*57. 14 

RD0*=P0T**3» 

IF ( RDO«LT *50« ) RD0»5O • 

D0B»<P0V-P0T)*3161./RD0 
M02=A0M*02M*( l.-(8.OOOl-P10)**3»/512. ) 

002*002+ (D0B*M02}*| 1 . -EXP ( -I/Z4 ) ) 

P0T*Q02*. 00333 
P10=P0T 

IF ( POT »GT»8 • )P10*8. 

P0D=POV-P0R 

P0B=P0B+( P0K+P0D+1* -P0B ) /Z 

IFtP0B.LT. .2)P0B=. 2 

ARl«ARl+( P0B-AR1 )*(1.-EXP(-I/A1K)) 

ARM*AR1+AK2*AR3 

C — ............................ ...... 

C AUT0REGULATI0N# INTERMEDIATE 

C ............................................................ 

P0A*P0A+( P0N+P0D+1 .-POA )/Z 
IF ( P0A.LT. .5)P0A**5 
AR2»AR2+(P0A-AR2) + ( 1 . -EXP ( - 1/ A2K ) ) 

C-- — ..... ............... 

C AUT0RLGULATI0N* LONG-TERM 

IF <P0DI194,192,192 
192 Pec=P0Z*P0D+l» 

GO T 0 196 

194 P0C -P0Z+P0D* • 33 + 1 * 

196 IF (P0C*LT. .3)P0C*.3 

AR3 = AR3 + < P0C-AR3) +I/A3K 

RETURN 

END 

SUBROUTINE ADH (AH > AHO AHK* AHM,» AHY/ AHZ> AH7> AN8> AUPjCNA# CNB> CNR, 

* CNZ,I » PRA-* Z ) 

REAL 1 
C 

C ANTIDIURETIC H0RM0NE 

C 


CNB*CNA-CNR 
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AHZ* » 2*PR A 

AHY«AHY+< AHZ-AHY)*. 0007*1 
AH8-AUP-1 * 

I F { AH8 * LT • 0 • > AH8*0. 

IF ( CNB • LT * 0 * )CNB*0. 

AH«AH + < CNZ*CNB+AH8-AHZ+AHY«4H >/Z 
lF(AH.LTt0. >AH«0. 

AHC-AHC+( .3333*AH-AHC)*U«-EXP<*I/AHK) ) 

AHM-6. * (1. -EXP (-0*1 8084 AHC) > 

IF( AHM.LT* .3) AHM-.3 

RETURN 

END 

SUBROUTINE MISCi ( AHM* AU4* AU8* I *SR *SRK*STH* TVD* TVZ*VEC*VIC*VTW* 

* VVE*VV6*VV7*Z ) 

REAL I 

C 

C************** ****************** ***************************************** 

c 

C VASCULAR STRESS RELAXATION BLOCK 
C 

C ****************************************************** ******************‘ 

VV6»VV6+(SR*< VVE-.301 > -W7-VV6 ) /Z 
VV7aW7+VV64{ l.-EXP(-I/SRK> ) 

Q****** ****************** ****************************************** ****** i 

c 

C THIRST AND DRINKING BLOCK 

C** ********************************************************** ****** ****** 1 
C ' 

TVZ«< • 01 *AHM“ • 009 ) ♦STH 
TVD*TVO+(TVZ-TVD>/Z 
IF ( TVD*LT »0« ) TVD*0» 

VTW-VIC+VEC 

C**** ****** ****** x***** ****** ************** **** ****************** ********> 

C 

C AUTONOMIC CONTROL BLOCK 

C ADAPTATION OF 8AR0RECEPT ORS 

C 

c********** ******* ********* ****** ********* ************** ****** ** ********* t 

AU4=AU4+AU8+I 

RETURN 

END 

SUBROUTINE HEART < AUR* OHM* HMD * HR ,1 ,PA * PMC* PMP., PMS* POT* PRA * Q AO* 

* QL0*RTP*SV0*VAE*VLE*VPE*VRE*VVE) 

REAL I 
C 

C HEART HYPERTROPHY OR DETERIORATION BLOCK 

C 

C— 

C HEART VICIOUS CYCLE 

C— — ........ — -- — 

DHM=< POT-6. >4.0025 

hmd«*hmd+dhm*i 

IF (hMD.GT.l. >HMD«1» 

O-- -- 

C MEAN CIRCULATORY PRESSURES 

C-— — ' 

PMC»( VAE + VVE+VRE + VPE+VLE)/.U 
PMS=( VAE+VVE+VR£>/. 09375 
PMP«( VPE+VLE)/* 01625 
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********4**#** ******** **¥.*¥¥*****♦** ******************** * ***** ****** I 

C 

C HEART RATE AND STR0KE VOLUME BLOCK AND TOTAL PERIPHERAL RESISTANCE 
C 

******************************* ************************************* 
HR*<32*+40.*AUR+PRA*2»)*< < HMD-1 •) 4.5+1 » > 

RTPMPA-PRA >/QA0 

SV0=*QLO/HR 

RETURN 

END 

SUBROUTINE CAPMBD ( BFN* CPC/ CPI * CPP> DFP , l >IFP*PC * PCD* PIF/ PLD> PPCj» 

* PRP*PTC*PTS.*PTT>PVGiPVS.*RVS.*TVD*V6 /VID/VIF* VP, 

* VPD* VTC*VTD*VTL*VTS*VUD*Z ,Z1 ,FUN6> 

REAL I* IFP 

C 

C CAPILLARY MEMBRANE DYNAMICS BLOCK 
C 

130 PTT«( VTS/12. )**2» 

VIF«VTS-VG 

CALL FUNCTN < VIF, PTS, FUN6 ) 

PIF=PTT-PTS 

CPI»IFP/VIF 

PTC».25*CPI 

CPP* PRP/VP 

PPC*.4*CPP 

PVG=RVS*1«79*BFN 

PC= PVG+PVS 

PCD*PC+PTC-PPC-PIF 

VTC«VTC+{ CFC*PCD-VTC) /Z 

PLD=»7*8+PIF-PTT 

VTL«=VTL + { »004*PLD-VTL >/Z 

IFIVTL.LT.O. ) VTL«0. 

VTD*VTC-»VTL-VID 

VTS»VTS+VTD*I 

VPD=VPD+(TVD-VTC+VTL-VUD-DFP-VPD)/Zi 

RETURN 

END 

SUBROUTINE PULM0N ( CPF* CPP* CPN*DFP* I *PCP*PFI*PLA*PLF*POS*PPA*PPC* 

* PPD*PPI*PPN*PPO*PPR*VP *VPD*VPF*Z *Z3) 

REAL I 

C 

C PULMONARY DYNAMICS AND FLUIDS BLOCK 
C 

VP »VP+< VPD*I )/Z3 
C 

200 PCP-.45*PPA+.55*PLA 

PPI-2.-.150/VPF 
CPN«sPPR/VPF 
P0ScCPN4.4 
PLF*{PP1+11. )*.0003 
PPO»PLF*CPN . 

PPN=( CPP-CPN)*. 000225 
PPD*PPD+t PPN-PP0-PPD)/Z 

IF(PPR+PPD*I-»025»LT»0. )PPO»< • 025*»PPR)/I 

PFI«(PCP-PPI+POS-PPC)*CPF 

DFP«DFP+< PFI-PLF-DFP ) /Z 

I F ( VPF+DFP*I-«001*LT.0. IDFP*« .001-VPF>/I 
VPF«VPF+DFP*I 

PPR»PPR+PPD*I 
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RETURN 

END 

SUBROUTINE MISC2 < HPL/HPR/ HSL* HSR/ I / PA/ PP A/ POT/ STH/ ZlQ* Z1 1/ Z13 ) 

REAL I 
C 

C******************************** 1 

c 

C HEART HYPERTROPHY OR DETERIORATION BLOCK 
C 

C****** **************** **************** ****** ***************************** 
HPL-HPL+< { <PA/100./HSL)**Z13)-HPL)*I/57600. 
hPR«HPR+( ( ( PPA/ 15 • /HSK )**Z1 3 )»HPR)* 1/57600. 

C ******************************************************************** ***** 

c - • 

C TISSUE EFFECT ON THIRST AND SALT INTAKE 
C 

C ******************************************** ************* **************** 
STH*<Z1O-P0T)*Z11 
IF ( STH.LT • 1 * ) STH=1 « 

I F ( STH » GT • 8 » ) STH=8 • 

RETURN 

END 

SUBROUTINE PROTEN ( CHY/ CPG/ CPI / CPK/ CPP/ CPR/ CPi /DLP/DLZ/DPC/ DPI / DPL/ 

* DPP/DPY/GPD/GPR/ I /IFP/LPK/PC /PCE/ PGX/ PRP/ VG / 

* VTL'Z /DP0/PPD) 

REAL I/IFP/LPK 

C 

C TISSUE FLUIDS/ PRESSURES AND GEL BLOCK 
C 

C — * ...................... 

C PLASMA AND TISSUE FLUID PROTEIN 

C .... . 

135 DPL®DPL+( VTL*CPI-DPL)/Z 
IF ( PC* LT«0 » ) PC®0 • 

DPC»DPC+( CPK* (CPP-CPI )*PC**PCE-DPC)/Z 

DPI«DPC-DPL 

DLZ«LPK*(CPR«CPP) 

IF ( CPP.GT.CPR )DLZ*4.*DLZ 
DLP»DLP+ ( DLZ-DLP > /Z 
PRPcPRP+{DLP-DP0+DPL-DPC-PPD>*I 

C — — — ... .................... „ 

C GEL PROTEIN DYNAMICS 

C — mm mm 

141 PGX»CHY**2* .01332*CPG+CPG 

GPD»GPD+( . 0005* { CPI -PGX)*VG-GPD)/Z 

GPR*GPR+GPD*I 

lFP-IFP+< DPI-GPD)*I 

RETURN 

END 

SUBROUTINE KIDNEY ( AAR/AHM/ AM / APD/ ARF/ AUM/CNE#CNX/CNY/GBL*GFN/GFR/ 

* GF2/GF3/GF4/GLP^ I /NAE/NED/NID/N0D/N0Z/PA /PAR/ 

* PFL*PPC/RBF/REK/RFN/RR /STH/TRR/ VIM/VUD/Z ) 

REAL I/NAE/NED/NID/NOD/NOZ 

C 

C KIDNEY DYNAMICS AND EXCRETION BLOCK 
C 

142 GF3=< (GFN/. 125-1. >*GF4)+1 . 

IF ( GF3.GT *15. >GF3*15. 

IF(GF3»LT • •4)GF3***4 
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AAR*31.67*VlfHM AUM* ARF+1 • -ARF >*GF3 
RR«AAR+5l • 6 6* VIM 
PAR*P A**GB L 
RFN=P AR/RR 
RBF=REK*RFN 
150 APD=AAR*RFN 
GLP«PAR-APD 
PFL-GLP-PPC-18. 

gfi=gfn 

GFN»GFN+(PFL*»00781“GFN)*GF 2/Z 
IF (ABS(GFN-GF1 > .GTt .OO2>G0 TB 142 
GFR»GFN*REK 

TRR«.8*GFR+.025*REK-.001*REK/AM/AHM 
VUD-VUD+ t GFR»TRR-VUD ) /Z 
I F ( VUD»LT. . 0002 JVUDc 0002 

C — ............ ............. ...................... 

C KIDNEY SALT BUTPUT AND SALT INTAKE 
C (SEE ALSB ELECTROLYTES AND CELL WATER BLOCK ) 

N0Z*1OOO.*VUD/AM/(CNE/CNX+CNY ) 

N0D*»N0D + ( NOZ-N0D » /Z 

NED»N1D*STH-N0D 

NAE=NAE+NED*I 

RETURN 

END 

SUBROUTINE IONS (AM #CCD*CKEjCKI*CNAjI >KCD>KE jKED^KI >KID#KlE-» 

* KXR.>K0D>NAE#REK>VECjVICjVID*VP /VPFaVTS/Z) 

REAL I/KCD^KE/KEDjKI/KID/KIEiKIRjKOD^NAE 

C 

C ELECTROLYTES AND CELL WATER BLOCK 
C 

160 VEC*VTS+VP+VPF 
CKE=KE/VEC 

KOD»( •00042*CKE+.00014*AM*CKE)»REK 

KIR=2850#+140»*CKE 

KIE*KIR-KI 

KCD*KCD+(KIE*»013-KCD >/Z 

Kl»KI+KCD*I 

KED=KID-KCD»KOD 

KE“KE+KED*I 

CKI*=KI/VIC 

CNA*NAE/VEC 

CCD-CKI-CNA 

VlD«VID+( .01*CCD-VID)/Z 

V IC*V IC+VID*I 

RETURN 

END 

SUBROUTINE GELFLD (CHY/CPG/CPI>GPRjHYL# IFP#PGCiP6h.»PGP.»PGR,PGXjPIF.> 

♦ VTS/PRM/PTC*PTS>PTTjVG /VGDjVIF#VRS/V2D>FUN6) 

REAL IFP 

C 

C GEL FLUID DYNAMICS 
140 CHY«HYL/VG 

PRN«*«5» 9*CHY+24*H 

PGR«.4*CHY 

CPG»GPR/VG 

PGP»*25*PGX 

PGC»PGP+PGR 

VIF»VTS«VG 
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649 


CALL FUNCTN ( VIF# PTS# FUN6 ) 

650 


PlF-PTT-PTS 

651 


CPI-IFP/VIF 

652 


PTC-.25+CPI 

653 


PGH-PIF+PTS+PRM 

654 


VGD*V2D*( PIF+PGC-PTC-PGH) 

655 


VG-VG+VGD 

656 


I F < VG »LT* 0 » ) VG-0 « 

657 


IF{ •012*LT«ABS(VGD) > G0 TO 140 

658 


RETURN 

659 


END 

660 


SUBROUTINE FUNCTN < TH# POL# TAB ) 

661 


DIMENSION TABU4) 

662 


N-14 

663 


DO 110 I=1#N#2 

664 


IF(TAB<I>-TH) 110# 120# 1 10 

665 

110 

CONTINUE 

666 


GO TO 140 

667 

120 

POL-TAB ( 1+1 ) 

668 

130 

RETURN 

669 

140 

NN-N-2 

670 


DO 150 I » 1 # NN# 2 

671 

150 

IF(TAB(I) • LT • TH .AND. TAB(I+2) *GTt TH) GO TO 160 

672 


WRITE ( 6# 100 ) TH 

673 

100 

FORMAT ( 5X# 1 ***** CURVE LIMITS EXCEEDED ***** »#G12»6//) 

674 


IF(TH .LT. TAB ( 1 ) ) P0L»Ta8(2> 

675 


I F ( TH .GT. TAB ( N-l > ) POL-TAB(N) 

676 


GO TO 130 

677 

160 

P0L“T AB ( I +1 >+(TAB( I +3) -TAB < 1+1 ) )*(< TH-TAB ( I > )/<TAB< 1 + 2) -TAB ( I ) ) ) 

678 


GO TO 130 

679 


END 

$80 


SUBROUTINE INPUT 

681 


COMMON/ARRAY/ A ( 400 ) 

682 


DIMENSION N V < 50 ) 

683 


DATA NYES/'YE'/ 

684 

1 

WRITE { 102# 100 ) 

685 

100 

F 0RMAT ( //2X# ' *** INITIALIZE# CHANGE# OR OUTPUT DATA (YES OR NO)*/) 

686 


READ( 101# 200 ) NOP 

687 

200 

FORMAT ( A2 ) 

688 


IF ( NOP »NE» NYES) GO TO 999 

689 


WRI TE ( 102# 250 ) 

690 

250 

F0RMAT(2X# •*** INPUT CODE*) 

691 


READ ( 101# 300 ) KEY 

692 

300 

FORMAT ( 13 > 

693 


GO TO ( 10#20#30#40)#K£Y 

694 

10 

WRITE ( 102#400 ) 

695 

400 

FORMAT ( //2X# '*** INPUT INITIAL DATA ***'//) 

696 


CALL PUTIN 

697 


WRITE < 102# 450 ) 

698 

450 

FORMAT ( 2X# * *** YOU MAY CHANGE INITIALIZED DATA IF DESIRED') 

699 


GO TO l 

700 

20 

WRITE ( 102#500 ) 

701 

500 

F0RMAT(2X# '*** NUMBER OF VARIABLES TO BE INPUT ( 13 )'// ) 

702 


READ ( 101# 300 ) NVAR 

703 


WRITE (6#600) 

704 

600 

F0RMAT(//2X# 'VARIABLES TO BE CHANGED') 

705 


WRITE ( 102# 650 ) 

706 

650 

F0RMAT(6X # '*** ARRAY NUMBER# VALUE ( 13# E13 . 6 )• ) 

707 


DO 25 I»1#NVAR 
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708 


READ t 101# 700 ) N/V 



709 

700 

F6RMAK I3/£13»6 > 



710 


WRITE ( 6/ 800 ) Ni V 



711 


WRITE ( 102*800) N/V 



718 

800 

FORMAT <5X/ '*** A(*/I3/')«'/El3»6) 



713 

25 

A(M = V 



714 


G0 T0 1 



715 

30 

WR I TE 1 1 02/ 900 ) 

- IF LESS THAN 

1/RET 

716 

900 

F0RMATI2X/ U*» VARIABLE T0 BE PRINTED 0UT(I3) 

717 

*UWN«//) 



718 

35 

RE AO ( 10 Is 300 ) N VAR 



719 


I F < NV AR «LT • 1) G0 T0 1 



720 


WRITE 16/902 ) 



721 

902 

FORMAT ( //2Xz * VARIABLES TO BE PRINTED 0UT'> 



722 


WRITE < 102*800) NVAR/ A ( NVAR ) 



723 


WRITE ( 6/800 ) NVAR/A(NVAR> 



724 


60 T0 35 



725 

40 

WRITE ( 102/ 903 > 

LINE PRINTER'//) 


726 

903 

F 0RMAT ( //2X/ ' *** COMPLETE ARRAY IS OUTPUT 0N 


727 


CALL PUT8UT 



728 


GO TO l 



729 

999 

RETURN 



730 


END 



731 


SUBROUTINE PUTIN 



732 


C8MM0N/ARRAY/A( 400) /TITLE (400) 



733 


WRITE ( 6/300 ) 


//) 

734 

300 

F0RMAT( 1H1/39X/ ' ***** INPUT PARAM 

E T E R S ****** 

735 

1 

READ (5/ 100 ( X/N/T 



736 

100 

FORMAT ( El 3. 6/ 2X, I5/2X/A4) 



737 


IF ( N tLT • 1) RETURN 



738 


A ( N )»X 



739 


T ITLE ( N ) »T 



740 


WRITE (6/400) N, A (N )/ TITLE( N) 



741 

400 

FORMAT (50X/ ' A( '/ 13 / * > */Ei2»6/» - */A4) 



742 


GO TO 1 



743 


END 



744 


subroutine putout 



745 


C0MM0N/ARRAY/AI 400) /TITLE (400 ) 



746 


WRITE ( 6/ 100 ) A ( 1 ) 



747 

100 

FQRMAT l //2X# ****** autPUT AT 'j F l0*4i» MINUTES ***** /7 ' 


748 


WRITE (6/200) (I/A (I )/ TITLE! I)/ 1-1/378) 



749 

200 

F0RMAT( 3( 15X/ 1 A ( • /I3j • ) »/E12»6/< - »/A4)) 



750 


RETURN 



751 


END 
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APPENDIX D 



31 


1 C 

2 C 

3 C 

4 C 

5 C 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 
3b 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 C 


R67093;8 LARRY NEAL — MATH; NASA FILE 


PROGRAM GUYTON 

CIRCULATORY DYNAMICS - CIRCE 
CIRCE1 

REAL LVM; I;IFP;LPD;KE;KEl.*K0D.»KIR;KrE;KI;KCD;K£D;KNl;KN3 
REAL NAE; NED; NID;NOD; I 1 ;LPK; K ID;M02; N8Z; KCZ; HPL> HPR; 12; 13; MMO 
DIMENSION FUNK 14 ); FUN2( 14) ;FUN3< 14 );FUN4 ( 14 ) ;FUN6 ( 14 I $ FUN7 ( 1 4 ) 
COMMON/ ARRAY/T; It VBD; VVS ; VP A/ VAS; VLA; VRA t VA Et PA;PAM;LVM; 

* VRE; PR A; URN* VP£, PP At PP It CP At RPA; RVM; VLE; PLA; QLN; PL1; 

* A1B;RPV; RPT;PGL;QP0;Q5 ; WE; VV8; PVS; PG V; RVG; GV0; AVE 
COMMON/ ARRAY/ CN2; CN3; RVS; PGS; RTP; QA0; OR0; QLO; DVS; DP A; DAS; DLA; ORA; 

* PA1; AUC;AUB; AUN;AU6;AU2;AU8;DAU;AUJ;AU ; AUO; AUH; VV 4; 

* AU9; AUM; AU4; VIF; P01; PTT; PTS; PIF; CPI* PTC; CPP; PPC; PVG 
COMMON/ARRAY/PC ; PCD; VTC; PLD; VTL; VTD; VPD; DPL; CPI; DPC; DPI; LPD; DLP; 

* DPP;CHY;PRM;PGR;CPG;PGP;GF1;PGX;PGC;PGH;PG2; VGD;VG ; 

* EPH; GP1; GP2; GPD; AAR; RR ; RFN; APD; GLP-* PFL; GFR; TRR.* VUD 
C0MM0N/ARRAY/REK;NOD;NtD;NAE; VEC; CKE; KOD; KE1 ; KI R; KI E; KCD; KED ; ZKlt 

* CNA;CCD; VID;KE ;KI ;VIC;I1 ;VTY;2 ; VTZ; VUZ; TV2; PPZ; 

* OFZ;X ;I2 ;PR1;VTS;VP ; PRP; IFP; GPR; KN3; KN1; AMR; AMP 
C0MM0N/ARRAY/AM1; AMC; AM2; AM3; AM5; AM ; CNE; AGK; ANP; AN1; ANC; AN2; AN3; 

* AN5; ANM; VB ; HM 1 ; HM ;B1 ; V I E; VI B; V I M; RC2; P02; RKC; RC 1 ; 

4 RCD; VRC; RSN; OVA; BFN; D0B; A0M; PI 0; OSV; POT; P0D; P0B; AR1 

C0MMON/ARRAY/AR2;P0C;AR3; ARM;CN8;GFN; AH7; AH8; AH ; AHC; AHl; AH2; AH4; 

* AHM;CNY;CNX; VV1; VV2; VV5; W6; VV7;TV0; VTW;HSR;HSL;NID; 

* SR ; VVR; RAR; CV ; CN7; AUX; AUK; AUZ; Y ; CFC; CPK; PCE; CPR 
COMM0N/ARRAY/LPK;DP8;HYL;KID; AMT; ANT;Pf!K;P0N> AIK; A2K; A3K;CNR; CNZ; 

* AHK; SRK; TM ;V2D;Z1 ;Z 2 ; Z3 ;Z 4 ; Z5 ; Z6 *Z7 ;Z8 ; HMK; 

4 HKM; P0 V; PO Z; RDO; Q02; RBF; M02; PDA; P0Y;ANU;P0R;GF2; MMD 

COMMON/ ARRAY/ DHM; POQ; 13 ; U ;VPl;Tl ; GF3; GF4; AUP; AU V; RV 1 ; AU Y; OUT; 

* DSP; AHZ; AH Y; OSA; PP I; CPN; POS; PLF; PPO; PPN; PPD; PF Ii OFP; 

* VPF;PPR;PMC;PMS;PMP;HR ; CPF; PCP; DA 1; DLZ; DPY; DPZ; GPZ 
COMMON/ ARKAY/NOZ;KCZ; VIZ; HPR; HPL; STH; ALO/EXC; 02M; PA2;PP2;SV3; AUL; 

* VV9; 02A; U1 ;EXE; ARF;ORF;RSM;BFM;RAM;OVS;PVO;RM0;QOM; 

* PM0;P2O;MM0;PD0;POE; AMM; A4K;PPM;RMM;PM1;PM3;PM4;EX1 
C0MM0N/ARRAY/Q2 ; Q3 ;PM5;PK1;Z9 ; ZlO; Z1 1; Zl 2; Z1 3; Z1 4; Z15; Z16; PK2; 

* PK3;FIS;bTA; PAR;GBL; ANY; ANZ; ANX; AN V; ANW; ANR; AUQ; AUR; 

* AUS; A3 78; DUMMY (22 ) ;TITLE( 400 > ; DUMNY ( 40 > 
C0MM0N/NUMBER/NOUMY (22) ;DMMY 

DATA FUN1 <1 );FUN1 (2);FUN1 (3);FUN1(4)>FUN1 (5 );FUN1 (6 )/FUNl(7 ); 
*FUN1(8);FUN1(9);FUNH10);FUN1(11);FUNU12>;FUN1(13)>FUN1<14>/ 

*0*; 1*G4;60*;1* 025; 1 25 • ; . 9 7; 1 60 . ; » 88; 200 • t • 59; 240 • ; 0 • t 240 . ; 0 • / 

DATA FUN2 ( 1 ) ; FUN2 ( 2 ) ; FUN? < 3 );FUN2(4 ) ; FUN2 ( 5 > ; FUN2 ( 6 ) ; FUN2 ( 7 >; 
*FUN2( 8);FUN2(9);FUN2( 101;FUN2( 11 |;FUN2( 1 2 ) ; FUN2 ( 1 3 ) ; FUN2 < 14)/ 
*-100«;0*0;-6.;0*0.»-3»,**7 t 5;-l..»2»6;2*.*9*8;8.;13.5; 1000 •; 13*5/ 

DATA FUN3 (1 >;FUN3(2 );FUN3 (3 );FUN3(4 );FUN3(5 >;FUN3(6 );FUN3 (7 ); 
*FUN3< 8) ;FUN3(9) ;FUN3( 10); FUN3( 11 ) ;FUN3( 12 ) ; FUN3 ( 13 ) *FUN3 ( 14 ) / 

* 0 • 0; 1 »06;20»; *97; 24 »; »93;30»; • 8; 38 « ; ♦ 46* 45 • ; 0 ♦ ; 45 • ; 0. / 

DATA FUN4 ( 1 );FUN4 (2 );FUN4 (3 >;FUN4( 4 );FUN4 (5 );FUN4(6 );FUN4(7 ); 
*FUN4( 8) ;FUN4( 9);FUN4( 10 );FUN4( 11 );FUN4( 12 );FUN4< 13)/FUN4( 14 )/ 

*— 100*;0*;— 4»;0.;— 1«;3»6;3*; 9.4;6*;il*6; 1 0 • ; 13«5>1000**13«5/ 

DATA FUN6 ( 1 ); FUN6(2 );FUN6(3 );FUN6< 4); FUN6(5 );FUN6(6 );FUN6(7 ); 
*FUN6< 8) ;FUN6( 9 ) ; FUN6 ( 10 J ; FUN6 ( 1 1 );FUN6{ 12 );FUN6( 13);FUN6( 14 )/ 
*-100.; 10000*; Q*;70.; .4; 9, 3; *8>3*3;1*2#1*3;1*6;.43;100*;0*/ 

DATA FUN7 ( 1 ) ; FUN7 < 2 ) ; FUN7 { 3 ) ; FUN7 ( 4 ) ; FUN7 ( 5 > ; FUN7 ( 6 ) / FUN7 ( 7 ); 
*FUN7( 8) ;FUN7( 9) ;FUN7( 10 );FUN7( 11 ) ;FUN7< 12 ) ; FUN7 ( 13 ) *FUN7 ( 14 >/ 

*0 •; 7*;3Q*;6*25; 60 * ; 3»; 100*; l.*160*; *15; 400.; .05*400*/ *05/ 
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59 

WRITE <102/5) 


60 

WRITE (6/S) 


61 

5 FORMAT (/* GUYTON MODEL FROM WHITE’/ 

62 

* ' REFER TO GE-AGS USER GUIDE TIR 741-MED-301 7 • // ) 

63 

90 CALL PUTIN 


64 

C 


65 

C90 CALL INPUT 


66 

C 


67 

I F < I .GT. 0 *b ) 3 

>0.5 

68 

100 I F ( OUT tEG# 30 

CALL PUT9UT 

69 

C 


70 

C100 I F ( OUT *EQ« 30 

CALL OUTPUT 

71 

C I F ( DSP »EQ# 3 • J 

CALL DSPLAY 

72 

C 


73 

T«T+I2 


74 

C 


75 

CALL HEM0 

( AMM/ ANM/ ANU/ ANY/ ANZ/ ARM/ AUH/ AUM/ AUYz AVE/ BFM/ BFN/ 

76 

¥ 

CN2/CN3/CN7/CV / DAS/ DLA/ DPA/ DRA/ DVS/ FI S/ HMD/ HPL/ 

77 

* 

HPR/HSL/HSR/ 12 /LVM/PA /PAM/PA2/PC / PGL/ P6S/ PL A/ 

78 

* 

PPA/PPl/PP2/PRA/PRl/PVS/QAO/QLN/QL0/QPO/QRF/yRN/ 

79 

* 

QRO/ QV0/ RAM/ RAR/ RBF/ RPA/ RPT / RP V/ RSM/ RSN/ RVG/ RVM/ 

80 

* 

RVS/U /VAE/VASi VBO/ VIM/VLA/VLE/VP / VRA/ VPE/ VRA/ 

81 

* 

VRC/ VRE/ VVE/ VVR / VVS / VV7/ VV8/ X / FUNi/ FUN2/ FUN3/ 

82 

* 

FUN4 ) 

83 

C 


84 

120 CALL AUTO 

( AU / AUB/ AUC/ AUH/ AUJ/ AUK/ AUL/ AUM/ AUN/ AU0/ AUP/ AUQ/ 

85 

* 

AUR/ AUS/ AUV/ AUX/ AUZ/ AU4/ AU6/ AU8/ A1B/ DAU/ EXC/ EXE/ 

86 

* 

EX1/I2 /PA /PA1/POQ/P0T/P2Q/ STA/ VVR/ VV9/ Y / Z/ 

87 

* 

Z8 / Z12 ) 

88 

c 


89 

IF< I3»LE. 12 )G0 TO 168 

90 

IF < ABS( DAU-AUJ) . 

GT.DADGO TO 100 

91 

110 IF UBS(QA0-QL0) 

. GT « • 2 ) GO TO 100 

92 

IF (ABS(UAO-QPO) 

• GT • » 2 ) GO TO 100 

93 

IF ( ASS ( QAO-ORO ) 

• GT • • 4 ) GO TO 100 

94 

C 


95 

168 CALL H0RM0N 

( AM / AMC/ AMP/ AMR/ AMT/ AMI / ANM/ CKE/ PA/ Z/ FUN7> 

96 

* 

AGK/ ANC/ANP* ANR/ ANT/ANV/ ANW/ AN1/CNA/CNE/GFN/ 

97 

* 

I /REK) 

98 

C 


99 

170 CALL BLOOD 

(HKM/HM /HMK/I /POT/P0Y/P01/P02/ RC1/ RC2/ RCD/ RKC/ 

100 

* 

VB / V IB/VIE/VIM/VP / VRC ) 

101 

C 


102 

180 CALL MUSCLE 

( AL9/ AMM/ AOM/ AUP/ A4K/BFM/EXC/HM /I / MMR/ 0MM/ 6SA/ 

103 

* 

0VA/OVS/O2A/PDO/PK1/PK2/PK3JPM0/PM1/PM3/PM4/PM5/ 

104 

* 

POE/P0M/PV0/P20^ Q0M/RM0. VPF* Z5 /Z6 ) 

105 

C 


106 

CALL AUT0RG 

( AOM/ ARM/ ARl/ AR2/ AR3/ A1K/A2K/ A3K#BFN/D0B/HM >1/ 

107 

¥ 

MO2/0SV/0VA/02M/POA/P0B/P0C/P3D/POK/P0N/POR/P0T/ 

108 

* 

POV/PQZ/P10/Q02/RO0/Z / Z4 /Z7) 

109 

C 


110 

CALL ADH 

< AH /AHC/AHK/AHM/AHY/AHZ/AH7/AH8/AUP/CNA/CNB/CNR/ 

111 

¥ 

CNZ/I /PRA/Z) 

112 

C 


113 

CALL MISC1 

( AHM/AU4/AU8/I /SR / SRK/ STH/ T VD/ T VZ/ VEC/ VIC » VTW/ 

114 

* 

VVE/VV6/W7/Z) 

115 

C 


116 

CALL HEART 

< AUR/DHM/HMDiHR /I /PA / PMC/ PMP/ PMS/ P0T/ PRA/ QA0/ 

117 

* 

QL0/RTP/SV0/VAE/VLE/VPE/VRE/WE) 
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118 

119 

120 
121 
122 
123 
120 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 


C 


C 


C 

C 

C 

C 

C 


C 

C 


C 


C 


C 


C 


C 

C 

C 

C 


130 CALL CAPMBO ( BFN*CFC*CPI/CPP*DFP* J *IFP*PC *PCD*PIF*PLD*PPC* 

* PRP*PTC*PTS*PTT*PVG*PVS*RVS*TVD*VG *VID*VIF*VP* 

* VPD*VTC* VTD*VTL*VTS*VUD*Z *Z1 * FUN6 ) 

I«I*1 •2+T-T1 
I1»ABS( VP1/VPD/I ) 

I F C 1 1 • L T'« I ) I'«U 

I F ( 1 3+T«T 1 • LT • I > I = I3 + T-T1 

T-I+Tl 

Tl-T' . 

IF ( OUT • EQ »4 • ) CALL PUT0UT 

IF( 0UT.EQ.<f # ) CALL OUTPUT 
I F ( DSP* EQ »4 • ) CALL DSPLAY 


200 CALL PULMBN 
* 

call MISC2 

135 CALL PR0TEN 
* 

* 

142 CALL KIDNEY 
* 

¥ 

160 CALL I 0NS 
* 

140 CALL GELFLD 
* 

GB T 0 100 

end 

SUBROUTINE HE MO ( AMM* ANM* ANU* ANY* ANZ* ARM, AUH* AUM* AUY* AVE* BFM* BFN* 

* CN2*CN3*CN7*CV * DAS* DLA* DPA* DRA* DVS* FIS# HMD# HPL* 

* HPR*HSL*HSR* 12 *LVM*PA *PAM*PA2*PC / PGL* PGS* PLA* 

* PPA*PP1*PP2*PRA*PR1*PVS*QA0*QLN*QL0*QPO*QRF*QRN*' 

* QR0* QV0* RAM* RAR* RBF* RPA* RPT* RP V* RSN* RSN* RVQ* RVM* 

* RVS* U *VAE*VAS* VBD* VI M* VL A* VLE* VP * VPA* VPE* VRA* 

* VRC* VRE* VVE* VVR* VVS* VV7* W8* X * FUN1*FUN2* FUN3* 

* FUN4 ) 

DIMENSION FUN1(14)*FUN2(14)*FUN3( 14)*FUN4(14) 

REAL 12* L VM 

CIRCULATORY DYNAMICS BLOCK 
HEMODYNAMICS 

vbd»vp+vrc-ws-vas-vla-vpa-vra 

VVS«VVS+DVS*I2+VBD*«3986 

VPA«VPA+DPA^I2+VBD**155 

VAS»VAS+0AS*I2+VBD**261 

VLA»VLA+DLA*I2+V8D**128 

VRA«VRA+DRA*I2+VBD»»0574 

VAE»VAS«* 495 


(CPF*CPP*CPN*DFP*I *PCP*PFI*PLA*PLF*POS*PPA*PPC* 
PPD*PPI*PPN*PP8*PPR*VP * VPD* VPF* Z *Z3) * 

(HPL*HPR*HSL#HSR* I* PA* PPA* POT* STH* Z10* Z1 1* Z13) 

(CHY*CPG*CPI*CPK*CPP*CPR*CP1*DLP*DLZ*0PC*DPI*DPL* 
DPO*DPY*GPD*GPR* I *IFP*LPK*PC *PCE*PGX*PRP* VG * 
VTL#Z *PP0) 

( A AR* AHM* AM * APD* ARF* AUM* CNE* CNX* CNY* GBL* GFN* GFR* 
GF2*GF3*GF4*GLP> I *NAE*NED*NID*N0D* N0Z*PA *PAR* 
PFL*PPC*RBF*REK*RFN*RR *STH* TRR* VIM* VUD* Z ) 

(AM *CCD*CKE*CKI*CNA*I *KCD*KE *KED*KI *KID*KIE* 
KIR*K0D*NAE*REK*VEC*VIC*VID*VP * VPF*VTS*Z) 

(CHY*CPG*CPI*GPR*HYL*IFP*PGC*PGH*PGP*PGR*PGX*PIF* 
PRM*PTC*PTS*PTT*VG * VGD* V IF* VRS* VTS* V2D* FUN6 > 



177 

178 

179 

180 
181 
182 

183 

184 

185 
186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 C 

234 C 

235 C 


PA»VAE/ *00355 

PAM*100*/PA 

PA2-PA/AUH 

CALL FUNC?N(PA8#LVM#FUN1 1 
VRE-VRA-. 1 
PRA-VRE/.005 

CALL FUNCTN ( PRA# QRN/FUN2 ) 

VPE«VPA-*30625 
PPA=VPE/«0048 
PP1».026*PPA 
IF (PPl*LT.O. )PPl-0. 

RPA*PPl**(-.5> 

PP2-PPA/AUH 

CALL FUNCTN IPP2#RVM#FUN3) 

VLE-VLA-. 4 
PLA-VLE/.01 

CALL FUNCTN(PLA#QLN/FUN41 

RPV-1 ./(PLA+20. 1/.0357 

RPToRPV+RPA 

PGL-PPA-PL A 

QP0*P6L/RPT 

ANU-ANh 

IF (ANU.LT* •8)ANU*«8 
VVE«VVS-VVR-< ANU-1. >*ANY 
VV8«VVE-W7 

IF ( VV8.LT • .0001 )VV8 = . 0001 

PVS=V V8/C V 

PR1-PRA 

IF (PRA.LT.O. )PR1«0. 

RVG«2»738/PVS 
QVP=< PVS-PR1 ) /RVG 

CN3«CN3+{ ( (PC-17. ) *CN7+17 • ) *CN2-CN3 > * • 1 
AVE«M AUM-1. ) *AUY+1 . 

RVS=AVE*< l./CN3)*VIM*( ( ANU-1* )*ANZ+1« ) 

PGS-P A-PVS 

RSN=RAR*ARM*ANU*AUM*PAM*VIM+RVS*1 .79 
BFN=PQS/RSN 

RSHsANU*VIM*PAM*aUM*AMM*RAM 

BFMbPQS/RSN 

QAO=BFN+BFM+RBF+( PA-PRA>*FIS 
QLP=LVM*QLN*AUH*HSL*HMD*HPL 

QRO=QRN*< ( 1 .-QRF)*AUH*RVM*HSR*HMD*HPR+QRF*QL0/QLN) 

QP0=QL0+( QP0-QL0 >/U 

GV0»QR0+( GV0-QR01/X 

DVS-QA0-QV0 

DPA»GR0-QP0 

DAS-GL0-GA0 

DLA-GP0-GL0 

DRA«QV0-GR0 

RETURN 

END 

SUBROUTINE AUT0 ( AU # AUB# AUC# AUH# AU J# AUK# AUL# AUM# AUN# AUO# AUP# AUG# 
¥ AUR# AUS# AUV# AUX# AUZ/AU4# AU6# AU8# AlB# OAU# EXC# EXE# 

¥ EX1#I2 # PA # PA1#P0Q#P0T#P20#STA# VVR/W9/Y # Z# 

* Z8 # Z 1 2 ) 

REAL 12 

AUTONOMIC CONTROL BLOCK 
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120 EXE-< 8.-P20 |*EX1+ (EXC-1* >*Z12 
P0Q=P9T 

IF (P0Q.QT.8. )P0Q«8# 

IF ( POO »L T • 4 * ) P0Q*4 • 

PAl«PA*P0Q/8#-EXE 
AUOO • 

I F ( PA 1 • LT » 80 * >AUC*.03*< 80.-PA1 ) 

IF(PA1.LT.40. )AUC*1»2 
AUB*0. 

IF(PA1.LT.170* )AUB««014286*(170.<-PA1) 

I F ( PAi • LT • 40 * ) AUB“1 *83 

123 A1B = ( AUB-1 • )*AUX + 1 . 

124 AUN«0 

IF(PAl.LT»b0. )AUN**2*<50#-PA1 ) 

IF ( PA 1 « LT • 20 • ) AUN«6 » 0 
AU6-A1B-AU4 
AU8“AUK*( AU6-1. ) 

0AU*DAU+( AUC+AU6+AUN-DAU) /Z/Y 
aUJ=AUJ+(DAU-AUJ)*I2*6./Z8 
IF ( AU J. LT * 0 * ) AU J®0 • 

IF( AUJ-1. > 126# 127/ 127 

126 AU«AU J**AUZ 
GO T0 128 

127 AU“ ( AUJ-1 •) *AUZ+1 * 

128 IF<STA.GT..00001)AU=STA 
AU0=AU-1. 

AUP=AU0*AUU+1 « 

AUH»AU0*AUV+1 • 

AUR*AU0*AUS+1 * 

VVR=VV9-AUL*AUP 

AUN=. 15+.85*AUP 
RETURN 

END 

SUBROUTINE H0RM0N ( AM t AMC# AMP# AMR # AMT# AMi# ANM# CKE # P A# Z# FUN7# 

* AGK# ANC# ANP# ANR# ANT t ANV# ANW# AN1 # CNA# CNE# GFN# 

* I # REK ) 

DIMENSION FUN7 ( 14 ) 

REAL I 
C 

£*********** * ****** ******** ******** ¥* ****** ** ****** ****¥***************** 


C 

C ALDOSTERONE CONTROL BLOCK 
C 

C * ** * ¥ ** * * * ** * **¥ ** ** * * * # ** * * ¥ * ¥ * * * *#¥¥ ¥ * * * ****** * * ** * * * * * * * * **** ** ** * * ** 


C 

c 

c 

c 

c 


168 AMR=CKE/CNA/. 00352-9. 

I F < AMR * LT »0 » i AMR=0» 

CALL FUNCTN < PA# AMP#FUN7) 

AMleAMl+< ANM*AMP*AMR-AM1 ) /Z 
AMC«AMC+( AM1-AMC ) *< 1. -EXP (- I /AMT) 1 
AM=20» 039-1 9. 8*£XP(«. 0391 *AMC ) 

*******************¥***********¥***¥* + ****¥********* * 4 ** 4 **'* ****** ****** 


ANGIOTENSIN CONTROL BLOCK 




i****** 


CNE*=152 * -CNA 
I F ( CNE.LT . 1 . ) CNEal • 

ANR*( ( 17.75-GFN*CNA UAGK+1. >*REK 
ANW»ANW+( ( ANR-1 • ) *10 • -ANW )*ANV*I 
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IF(ANW.LT».0)ANW**0. 

ANPoANR+ANW 

1 F < ANP • GT • 1 00 • ) ANP-100* 

IF( ANP.LT..01 )ANP».01 

AN1*AN1+ ( ANP- AN 1 >/Z 

ANC-ANC+t AN1-ANC) *( l.-EXP (- I/ANT] ) 

ANM<*4.0-3«3*EXP t-.0967*ANC) 

IF{ ANM.LT* *7) ANM«=*7 

RETURN 

END 

SUBROUTINE 8L00D (HKMaHM jHMK/I j POT* P8Y# P01 >P02 > RCl# RC2 # RCD> RKCj 

* VB • VIB/VIEa VIM* VP >VRC) 

REAL I 

C 

C RED CELLS AND VISCOSITY BL0CK 

C — — — — — — — — — —————— — —— — — — — — — — — mmmmmmm 

C BL00D VISCOSITY 

C— 

170 VB= VP+VRC 

- hm*ioo.*vrc/vb 

VlE»HM/(HMK-HM)/HKM 

VIB=VIE+1«5 

VIM=.3333»VI8 

c— ————————————————————————————— 

C RED BLOOD CELLS 

C — — — — —— — — — — — — — ——— — — — — ——— — — — — — —— 

RC2=RKC*VRC 

P02=PO1-POT 

IF{ P02.LT. .2375) P02»* 2375 

RC1*P0Y*P02 

RCDoRCl-RC2 

VRC=VRC*RCD*I 

RETURN 

END 

SUBROUTINE MUSCLE < AL0> AMMj A0M* AUP.» A4K, BFM# EXC> HM j I * MM0 > 0MM j 0SA> 

* 0VA> 0VS.» 02A.»PD0..PKi,PK2>PK3* PM0* PM1 * PM3 j PM4 # PM5.F 

* P0£/R0MjPV0/P20.*Q0M/RM0>VPF>Z5 >Z6) 

REAL I/MM0 

C 

C MUSCLE BLOOD FLOW CONTROL AND P02 BLOCK 
C 

180 0SA=AL0-VPF*.5 
0VA*0SA*HM*5. 

0VS*0VS+! ( BFM*OVA-RM0 )/HM/5./BFM«.0VS) /Z6 
PV0=57. 14*0VS 

RM0«( PV0-PM6 ) *PM5/ ( PMl4*PK3“PM4 ) 

Q0M*O0M+[RMe-MM0)*tl.»EXP(-I/Z5) ) 

PM0-PK2/( PK1-Q0M) 

PM1®PM0 

IF< PMl.LT.PM3 > Pi'll <*PM3 
P20»PM0 

IF( P20.GT »8« )P20»8* 

A0Mo(AUP-l. )*02A+1 * 

MM0 = AOM*0MM4EXC*< 1 8. 0001-P20 ) **3»/5l2» ) 

PD0aPV0-4O. 

P0E=P0M*PD0+1 • 

IF(P0E.LT* «OO5)P0E«»OO5 
AMM=AMM+( P0E«AMM)*( 1 . -EXPI-I/A4K) ) 

RETURN 
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381 
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383 

384 

385 

386 

387 

388 

389 

390 

391 
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END 

SUBROUTINE AUT0RG ( ABN/ ARM* ARl > AR2> AR3> A 1 K> A2K > A3K> BFN> DOB# HM > I> 

4 M02>0SV>0VA>02M>PeA>P0B>P0C>P0D>P0K>PeN>P0R>P0T> 

* P0V>P0Z>P10>Q02>RD0, 7 . >74 >Z7) 

REAL I>M92 
C 

C NON-MUSCLE OXYGEN DELIVERY BLOCK 
C AND NON-MUSCLE LOCAL BLOOD FLOW CONTROL BLOCK 


C AJT0RE'GULATION>RAP1D 


! PSV=0SV+< (BFN40VA-D0B >/HM/5*/BFN-0SV)/Z7 

P0V-0SV457. 14 
RD0=P0T443* 

IF( RD0.LT. 50. )RD0=5O* 

D0B*(P0V-P0T)*3161«/RD0 
M02<=A0M*02M*( 1>-(8.OOO1-P10)**3./512» ) 

Q02=Q02+(D0B-M02)4< 1 . -EXP { -I/Z4 ) ) 

P0T“Q 02 4*00333 
P1O=P0T 

IF(P0T»GT*8*)P1O=8. 

P0D-P0V-P0R 

P0Bs=P0B+( P0K4P0D+1 *-P0B ) /Z 
IF ( P0B*LT • « 2 ) P0B = *2 
AR1*AR1+(P0B-AR1)4(1*-EXP(-I/A1K) ) 

ARM»AR1*AR2*AR3 

c--— --------- — - ... — .............. , 

C AUT0REGULATI0N> intermediate 

P0A=P0A+< P0N4P0D+1 *-P0A )/Z 
IF (P0A*LT. *5 >P9A®*5 

AR2=AR2+< P0A-AR2 >4( 1 . -EXP ( -I/A2K ) ) 

C- . - — — — ,« — — . « . 

C AUT0REGULATI0N/L0NG-TERM 



IF ( P0D) 194> 192> 192 
192 P0C=P0Z*P0D+1 « 

G0 T0 196 

194 P0C=P0Z*P0D*«33+1* 

196 IF (P0C.LT* *3 )P0C=*3 

AR3-AR3+I P0C-AR314I/A3K 

RETURN 

END 

SUBROUTINE AD« (AH > AHC> AHK> AHM> AHY> AHZ> AH7> AH8> AUP> CNA> CN8> CNR, 
4 CNZ> 1 > PRA > Z ) 

REAL I 
C 

C ANTIDIURETIC H0RM0NE 
C 

CNB=CNA-CNR 

AHZ**24PRa 

AHY=AHY+( AHZ- AH Y ) 4*000741 
AH8=AUP-1 * 

I F ( AH8 ♦ LT ♦ 0 • ) AH8-0. 

IF ( CNB.LT.O* ) CNB®0* 

AH=AH+( CNZ4CNB+AH8-AHZ + AHY-AH )/Z 
I F l AH.LT.O. )AH®0. 

AHC*AHC+< *33334AH-AHC)4( 1,-EXP (®I/AHK) ) 

AHM-6*4 ( 1 « -EXP( -0* 18084AHC > ) 
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IF( AHM.LT . .3) AHM*.3 

RETURN 

ENO 

SUBROUTINE M1SC1 ( AHMs AU4; AU8' I >SR » SRtO STH » T VO > T VZj VEC* V IC> V TW t 

* WE* VV6*VV7*Z ) 

REAL I 

C 

C************** ******** ******** ********************************** *4 ***** 

C 

C VASCULAR STRESS RELAXATION BLOCK 
C 

c** ************************************************************** ******* 
VV6*VV6+( SR*{ VVE-*30l >-W7-W6 )/l 
W7=VV7 + VV6* ( 1 • -EXP l-I/SRK ) 1 

C************************ ****************************************** ***** 
C 

C THIRST AND DRINKING BLOCK 

C ************************************************************ *********** 
C 

TVZ*< *01*AHM-.009)*STH 
TVD=TVD+<TVZ-TVD)/Z 
IF(TV0*LT .0. ) TVDeO. 

V TW = V IC+VEC 

C******************** ******************************************** ******* 
C 

C AUTONOMIC CONTROL BLOCK 

C ADAPTATION OF BAR0RECEPTORS 

C 

C**** ********************************************** ********************* 
AU4*AU4+AUg*I 
RETURN 
END 

SUBROUTINE HEART ( AUR/ DHM/HMD* HR , I >PA jPMCjPMP/PMSjPOTjPRAjQAOj 

* QLO#RTPjSVO# VAEj VLE/VPE^ VRE/VV t ) 

REAL I 

C 

C HEART HYPERTROPHY OR DETERIORATION BLOCK 

C 

C— — — — — — ...... 

C HEART VICIOUS CYCLE 

0 »— — — 

DHM=<P0T»6. ) * *0025 

hmd=hmd+dhm*i 

IF {HMD.GT.l. )HMD»1* 

C----- - <* . - — — ...... 

C MEAN CIRCULATORY PRESSURES 

PMC«( VAE+VVE+VRE+VPE+VLE) /. 11 
PMS=( VAE+VVE+VRE) /. 09375 
PMP*{ VPE+VLE) /« 01625 

C****** ********************************************** ****** ************** 

c 

C HEART RATE AND STROKE VOLUME BLOCK AND TOTAL PERIPHERAL RESISTANCE 
C 

C4MMMMMHMF ************ ****** ********************************************** 

HR® ( 32 ♦ +40 ♦ *AUR+PRA*2 • )*{ (HMD-1 # ) **5+1. ) 

RTP 3 (PA-PRA )/QA0 

SV0®QL0/HR 

RETURN 
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2ND 

SUBROUTINE CAPMBD(6FN/CFC/CPI/CPP/DFP/I /IFP/PC /PCD/ PIF/PLD/PPC/ 

* PRP/,PTC/PTS/PTT/PVG/PVS/RVS/TVD/VG / V ID/ V IF/ VP/ 

* VPO/ VTC/ VTD/VTL/ VTS/VUO/Z /Z1/FUN6) 

DIMENSION FUN6 (14) 

REAL I/IFP 
C 

C CAPILLARY MEMBRANE DYNAMICS BLOCK 
C 

130 PTT*( VTS/12. )**2. 

VlFoVTS-VG 

CALL FUNCTN t VIF/PTS/FUN6 > 

pif=ptt«pts 

CPI “IFP/V1F 

PTC*.25*CPI 

CPP*PRP/VP 

PPC«.4*CPP 

PVG=RVS*1./9*BFn 

PC*PVG+PVS 

PCD=PC+PTC-PPC«PIF 

VTC«VTC+(CFC4PCD-VTC>/Z 

PL0«7*8 + PIf-*PTT 

VTL“VTL+l .004*PLD“VTL >/z 

IF( VTL.LT. 0. )VTL=0. 

vtd*vtc-vtl«vid 

VTS»VTS+VTD*I 

VP0*VPD+(TVD-VTC+VTL-VUD-0FP-VPD)/Z1 

RETURN 

END 

SUBROUTINE PULMON (CPF/CPP/CPN/ DFP/ I /PCP/PFI/ PLA/PLF / P0S/ PPA/PPC/ 

* PPD/PPI/PPN/PP0/PPR/VP /VPD/VPF/Z /Z3) 

REAL I 

C 

C PULMONARY DYNAMICS AND FLUIDS BLOCK 
C 

VP=VP + < VPD*I ) /Z3 
C 

200 PCP«.45*PPA+.55*PLA 

PPI»2.-ol50/VPF 
CPN*PPR/VPF 
P0S=CPN».4 
PLF = < PPI + H . ) * « 0003 
PP0=PLF*CPN 
PPN=( CPP-CPNJ*. 000225 
PPD=PPD+ ( PPN-PP0-PPO ) /Z 

I F < PPK+PPD*I-.025»LT.0. )PPD=< .025«PPR)/I 
PFI=( PCP-PPI+POS-PPC) *CPF 
DFP»DFP+(PFI-PLF-DFP J/Z 

I F ( VPF+DFP*I-»001 .LT.O. )DFP»< .001»VPF)/I 

VPF*VPF+DFP*I 

PPR«*PPR+PPD*I 

RETURN 

END 

SUBROUTINE MISC2 ( HPL / HPR/ HSL > HSR/ I / P A/ PPA/ P0T/ STH/ Z1 0/ Z1 1/ Z13 ) 

REAL I 
C 

C 44 * 4 ****** * * 4* 4 4 ** *♦ *4 44 44 4* * 4 44 4* 4 * **** 4 4 44 44 4* 44*4 ****** *4**4* 44 444* 44 

C 

C HEART HYPERTROPHY OR DETERIORATION BLOCK 
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C ************ ****** 44 ****** *** 4 ** 44 4 * ****** **♦**#*♦*» 4 ** 4 4 * ****** 44 4 ** 4 ** 

HPl=HPL+( ( (PA/100»/HSL)**Z13)-HPL>*I/57600« 

HPR=HPR+< < (PpA/15»/HSR)**Zi3>-HPRI*I/57600. 

C*4** ************ ***********************************,**.****#************* 

c 

C TISSUE EFFECT 0N THIRST AND SALT INTAKE 
C 

STH-( Z1O-P0T) *Z11 
I F ( ST H * LT • 1 • ) STH=1 . 

IF (STH*GT » 8 • ) STH*8 • 

RETURN 

end 

SUBROUT INE PR0TEN ( CHY> CPG#CPl .» CFtOCPP# CPR / CPl j DLPj DLZ.* DPCj DPI *DPL.» 

* DP0>DPY>GPD>GPR*I jIFPjLPKjPC * PCEiPGX/PRP>VG » 

* VTL/Z #PPD> 
real IjIFPjLPK 

c 

c TISSUE FLUIDS* PRESSURES AND GEL BLOCK 
C 

c- — .... 

C PLASMA AND TISSUE FLUID PR0TEIN 



135 DPL=DPL+< VTL*CPI-DPL ) /2 
IF < PC« LT »0 * ) PC = 0 • 

DPC=DPC+ < CPK* ( CPP«CPI ) *PC**PCE*DPC ) /Z 

DPI»DPC-0PL 

DLZ»LPK*(CPR-CPP) 

IF< CPP • GT »CPR )DLZ a 4«*OLZ 

DLP*DLP+(OLZ-DLP)/Z 

PRP=PRP+ ( DLP-DP0+DPL-DPC-PPD ) *1 

C — .............................. mm mmmmmmm .......... mmmmmmmmmmmmmmmmmmm. 

C GEL PR6TEIN DYNAMICS 

C — - — — . 

141 PGX=CHY **?* • o 1 332*CPG + CPG 

6PD*GPD + ( . 0005* < CPI -PGX ) * VG-GPD ) /Z 

GPR=GPR+GPD*1 

IFP*IFP+{ OPI-GPD ) *1 

RETURN 

END 

SUBR0UT I Nt KIDNEY ( AAR* AHM*AM » A PD/ ARF* AUM* CNE*CNX* CNY* GBL* GFN* GFR* 

* GF2*GF3*GF4*GLP* I >NAE*NED.»NIO*N0D.»N6Z*PA #PAR* 

* PFL* PPCjRBFjREKjRFNj RR t STH* TRR* V IM> VUD * Z ) 

REAL I>NAE*NED>NID*N0D*N0Z 

c 

c KIDNEY DYNAMICS AND EXCRETI 0 N BL 0 CK 
C 

142 GF3 = (( GFN/ • 125*1 •) *QF4 ) + 1 • 

IF(GF3«6T»15» ) GF3®1 5 * 

IF ( GF3«LT» «4) GF3=«4 
AAR=31.67*VIM*( AUM*ARF+1.-ARF)*GF3 
RR=AAR+51»66*VIM 
PARsPA-GBL 
RFN-PAR/RR 
RBF *REK *RFN 
150 APD-AAR4RFN 
GLP=PAR-APD 
PFL*GLP-PPC»18. 
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GF 1 *GFN 

GFN*GFN+( PFL*»00781-GFN }*GF2/Z 
IF UBS(GFN-GF1 ) .GT..OO2JG0 T0 142 
GFR»GFN*REK 

TRR*.8*GFR+.025*REK-.001*REK/AM/AHM 
VUD-VUD-M GFR-TRR-VUDJ/Z 
IF< VUD.LT*. 0002) VUD*. 0002 



C KIDNEY SALT OUTPUT AND SALT INTAKE 
C (SEE ALSO ELECTROLYTES AND CELL WATER BLOCK) 

C- — ----- 

N0Z = 1000.*VUD/AH/ (CNE/CNX+CNY > 

N0D»N0C«-< N0Z-N0D)/Z 

NED«NID*STh-N0D 

NA£*NAE+NED*I 

RETURN 

END 

SUBROUTINE IONS (AM /CCDjCKEjCKIjCNA/I t KCP/ K£ /KEDfKI /KID/KIE/ 

* KIR/KOD>NAE>REK.*VEC/VlC/VlD/VP /VPF/VTS/Z) 

REAL I/KCD/KE/KED/KI/KID/KIE/KIR/K0D/ NAE 

C 

C ELECTROLYTES AND CELL WATER BLOCK 
C 

160 VEC«VT$+VP+VPF 
CKE=KE/VEC 

K0D“ ( •00042*CK£+.00014*AM*CKE>*REK 
KIR°2850«+140»*CKE 
K I E*K IR*K I 

KCD®KCD+(KIE¥.013-KCD)/Z 

Kl=KI+KCO*I 

KED*KID-KCU-KOD 

KE*KE +ked*i 

CKI*.KI/V1C 

cna*nae/vec 

CCD=CKI-CNA 

VlD*VlD+( •01*CCD-VID)/Z 

VIC=VIC+VID*I 

RETURN 

END 

SUBROUTINE GELFLD t CHY/CPG/CPI/GPR>HYL/ IFP/PGC/PGH/PGP/PGR/PGX/PIF/ 

* PRMjPTC/PTS^PTT^ VG / VGD> VIF/VRS/VTSj V2D/ FUNS ) 
DIMENSION FUN6 ( 14 ) 

REAL IFP 
C 

C GEL FLUID DYNAMICS 
140 CHY«HYL/VG 

PRM*-5*9*CHY+24»2 

PGR«.4*CHY 

CPG-GPR/VG 

PGP«.25*PGX 

PGCaPGP+PGR 

VIF-VTS-VG 

CALL FUNCTN ( VI V, PTS/ FUN6 > 

PIF-PTT-PTS 

C p I “I FP/V IF 

PTC*.2b¥CPl 

PGH-PIF+PTS+PRM 

VGD»V2P* ( PIF+PGC-PTC-PGH ) 

VG=VG+VGD 
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1F( VG.LT.O. )VG*0« 

IF( .012*LT.ABS( VGD) > GO TO 140 

RETURN 

END 

SUBROUTINE FUNCTN ( TH/ POL* TAB ) 

DIMENSION TAB ( 14 ) 

N*1 4 

DO 110 I = 1 j N/ 2 
I F ( TAB ( I) »TH ) 1 10* 120/ 110 
110 CONTINUE 
GO TO 140 - 
120 POL-TABU + l) 

130 RETURN 
140 NN*N»2 

DO 150 I*l/NN/2 

150 IF < TAB ( 1 ) »LT • TH .AND. TABU+2) *GT t TH) GO TO 160 
WRITE ( 6/ 100 ) TH 

100 FORMAT I 5X/ ' ***** CURVE LIMITS EXCEEDED ***** *r G12«6//) 

IF ( TH .XT* TAB 1 1 ) ) P0L*TAB(2) 

I F( TH .GT'. TAB ( N-l ) ) POL»TAB(N) 

GO TO 130 

160 P0L=T AS ( 1 + 1 )+ (TAB l 1+3 >-TAB< 1 + 1 ) ) *1 < TH-TAB( I ) )/( TA8( 1+2) -TAB < I ) ) ) 
GO TO 130 
END 

SUBROUTINE PUTIN 

COMMON/ARRAY/A ( 400) /TITLE I 400 > / COL ( 20 ) / ALPHA ( 20 ) 

COMMON/ NUMBER/KjNOt 20 ) * NTIMEC/UNITS 
DATA ALL/'ALL 1 // BLANK/ • •/ 

DO 1 J“l/ 400 
A ( J J =0 « 

1 TITLEI J>=BLANK 

2 READ ! 5/ 100 ) VALUE/ NUMBER/ SYMBOL 

100 f-ORMAT (E13.6/2X/I5/2X/A4) 

IF(NUMBER.EQ.O) GO TO 3 
A (NUMBER) “VALUE 
TITLE (NUMBER) “SYMBOL 
GO TO 2 

3 READ ( 5/ 200 ) l ALPHA { J ) / J = 1 / 20 ) 

200 FORMAT (20A4) 

IF( ALPHA( 1 ) .NE. ALU GO T0 4 
READ ( 5/ 300 ) NTIMEC/UNITS 

WRITE (6/ 30) UNITS/ < TI TLE ( J ) / A I J > / J=»l/ 377 ) 

GO TO 31 

4 DO 5 K-l/20 

IF(ALPHA(K).EQ» BLANK) GO TO 6 

5 CONTINUE 

6 IF { K « LT • 20 ) K=K-1 
DO 10 J=»l/K 

L*1 

7 IF( ALPHA( J) .EQ.TITLE(L) ) GO TO 9 
L*=L + 1 

IF( L.LT.401 > GO TO 7 
WRITE ( 6/ 8 > ALPHA ( J > 

8 FORMAT ( 1H1////44X/ 'THE VARIABLE »A4/'WAS ILLEGALLY CALLED FOR.’) 
COLI J)-0. 

NO ( J ) »1 

ALPHA(J)“TITLE(1) 

GO TO 10 
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9 C0L(J)«A(L) 

N0 ( J ) «L 

10 CONTINUE 

READ ( 5/ 300 ) NTIMEC, UNITS 
300 FORMAT ( 15; IX, A4 ) 

IP(K.GT»9) GO T0 20 

GO TO ( 11, 12,1 3, 14, 1 5, 16, 17, 18, 19), K 

11 WRI TE ( 6; 21 ) UNI TS, ALPHA < 1 ), COL ( 1 ) 

G0 T0 31 

12 WRI TE ( 6, 22 ) UNI TS, < ALPH A ( J ) , J*i, 2 ) , < C0L ( L > / L-l, 2 > 

G0 T0 31 

13 WRI TE ( 6, 23 ) UNI TS, < ALPHA ( J ) , J«1 , 3 ) , ( COL ( L ) , L«l, 3 ) 

GO TO 31 

14 WRI TE ( 6; 24 ) UNITS, < ALPHA < J ) , J=i , 4 ) , ( C0L ( L )/ L«l, 4 > 

G0 T0 31 

15 wR I TE ( 6, 25 ) UNI TS, < ALPH A ( J > , U = i , 5 ) # < COL < L > , L«1 / 5 ) 

GO T0 31 

16 WRI TE ( 6/ 26 ) UNITS, < ALPHA < J ) , J«1 ; 6 ) , ( COL ( L > , L'l, 6 ) 

GO TO 31 

17 WRITE ( 6,27 ) UNI TS, ( ALPHA { J ) , J*l, 7 ) , ( C0L < L ) , L*l, 7 ) 

GO TO 31 

18 WRITE ( 6, 28 ) UNITS, ( AL p HA< J ) / J»l, 8 > , (COL ( L > / L«*l; 8 ) 

GO TO 31 

19 WRITE ( 6,29 ) UNI TS, < ALPHA ( J ) , J«1 ,9 ) , ( COL < L ) , L*l, 9 ) 

GO TO 31 

20 WRI TE ( 6; 30 ) UNI TS, ( ALPH A ( J J , COL ( J ) , J-l , K ) 

21 FORMAT ( 1H 1, 56X, A4, 7X, A4//59X, 1H0, 2X,E13»6 ) 

22 FORMA T ( 1H 1 , 49X; A4/3X, 2( 4X/A4, 6X )//52X,2H0 , 2 ( 1 X, El 3 • 6 > ( 

23 F0RMAT<lhl,42X,A4,3X,3l4X,A4,6X)//45X,2HO , 3( 1X,E13*6 ) ) 

24 F0RMAT< 1H1,35X,A4,3X,4(4X,A4,6X)//38X,2H0 , 4 ( 1 X, E 1 3 • 6 ) > 

2b FORMATt 1H1,28X,A4,3X,5<4X/A4,6X)//31X,2H0 , 5 < 1 X, E 1 3 . 6 > ) 

26 F0RMATUH1,21X, A4,3X,6(4X,A4,6X)//24X,2H0 , 6 < IX, E 13 • 6 > ) 

27 FORMAT( 1H1,14X,A4,3X,7(4X,A4,6X>//17X,2H0 , 7 ( IX, E 13 • 6 ) ) 

28 FORMAT ( 1H1, 7X, A4, 3X, 8 ( 4X , A4, 6X ) //10X, 2H0 , 8 ( 1 X, E 1 3 . 6 ) ) 

29 F0RMAT( 1H1, 1X,A4,3X,8(4X, A4,6X),4X, A4//4X/2H0 , 9 < IX, E 13 , 6 ) ) 

30 FORMAT { 1H1/63X/2H0 , A4//b ( 4X, A4,2H» ,E13.6,3X>> 

31 RETURN 
END 

SUBROUTINE PUTOUT 
C 

COMMON/ ARRAY/ A (400), TITLE (400) j COL (20), ALPHA (20) 
COMMON/NUMBER /K,N1,N2;N3,N4;N5;N6;N7;N8;N9;N10; 

* Nli,Nl2-»N13;N14;N15;N16;N17;N18;N19;N20; 

* NTIMEC, UNITS 

DATA SECS/* SECS */,TMlN/ 'MINS' //HOUR/ 'HOUR'/, DAYS/ 'DAYS'/ 
DATA NTlMEP/l/;N/0/; NN/ 1/;ALL/’ ALL V, BLANK/' »/ 

C E0UIVALENCE(A( 1 ),T) 

T = A ( 1 ) 

IF(UNITS.tQ.SECS) GO TO 2 
IF(UNITS.EQ*TMIN) GO T0 3 
IF(UNITS«EQ»H0UR> GO TO 4 
I F ( UN ITStEQtDAYS) GO TO 5 
WRITE (6,1) UNITS 

1 F0RMAT( 1H1, 47X, ' YOU CANNOT ASKED FOR TIME UNITS OF 'A4,'.'/ 

* 45X; 43HUSE EITHER “SECS" , « MI NS « , "HOUR" , OR "DAYS".) 
GO TO 66 

2 NTIME*=T*6Q. 

lF(NTIME.LT.NTIMEP) GO TO 65 
IF(NTIME.LT.(N+1)*60) G‘0 TO 6 
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800 
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807 

808 

809 

810 
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812 

813 

814 
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816 

817 

818 
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820 
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822 
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825 


N-N + l 

GR6SSU*TMIN 
G0 T0 6 

3 NTIME-T 

IFJNTIME.LT. NTIMEP) 60 T0 65 
IFJNTIME.LT* (N+l)*60) 60 T0 6 
N“N+1 

GR0SSU=>H0UR 
G0 T0 6 

4 NTIME«T/60. 

IFJNTIME.LT. NTIHEP) 68 T0 65 
IFJNTIME.LT. <N+1)*24) 60 T0 6 
N*»N + 1 

GR0SSU»DAYS 
G0 T0 6 

5 NTIME-T/1440. 

IFJNTIME.LT. NTIMEP) G0 T0 65 

6 IFJ ALPHAJ 1 ) .NE.ALL) G0 T0 7 

WRITE (6/50) NTI ME/ UNITS/ ( TI TLE ( U ) / A ( J ) / J«*l/ 377 ) 
60' T0 51 

7 G0 T0J3O/29/28/ 27/26/25/24/ 23/22/21/ 

* 20/ 19/18/ 17/ 16/ 15/14/13/ 12/ 11 )/K 

11 C0L ( 20 ) “A ( N20 ) 

12 C0L ( 19 ) =A ( N19 ) 

13 C0L ( 18 ) »A JN18 ) 

14 C0LJ 17)=AJN17> 

15 C0L J 16) -AJN16) 

16 C0L J 15) =A (N15 ) 

17 C0LJ14)«A(N14) 

18 C0LJ 13)»A(N13) 

19 C0L ( 12 ) =A ( N12 ) 

20 C0L ( 1 1 ) =A ( N1 1 ) 

21 C0L ( 10 ) <*A ( N10 ) 

22 C0L(9)=AJN9) 

23 C0L ( 8 ) = A ( N8 ) 

24 C0L ( 7 > «A ( N7 ) 

25 C0LJ6)=A(N6) 

26 C0L J 5 >«A ( N5 ) 

27 C0L ( 4 ) = A ( N4 ) 

28 COL ( 3 ) = A ( N3 ) 


29 

30 


31 


32 


33 


34 


35 


36 


37 


38 


C0L(2)*A(N2) 

C0LJ 1 )“AJ N1 ) 

IFJK.GT.9) 60 T 0 40 

60 TO < 31/32/ 33/34/ 35/36/37/38/39>/K 
WRITE ( 6/41 ) NTIME/C0LJ1) 

G0 T0 51 
WRITEJ6/42) 

G0 T0 51 
WRITEJ6/43) 

G0 T0 51 
WRITE (6/44) 

G0 T0 51 
WRITE ( 6/45 ) 

60 T 0 51 
WRITE ( 6/ 4b ) 

G0 T 0 51 
WRI TE ( 6/ 47 ) 

GO T0 51 
WR I TE (6/48) 


NTIME/ (C0L( J)/J«*l/2) 
NTIME/ (C0LJ J) / J“1 / 3 ) 
NTIME/ (C0L( J)/J»l/4) 
NTIME/ ( C0L ( J ) / J 3 l / 5 ) 
NTIME/ ( COL ( J ) / J*1 / 6 ) 


NTIME/ ( COL ( J ) / J a i / 7 ) 


NTIME/ (C0LJ Jl/J-1/8) 
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) 


G0 T8 51 

39 WRITE ( 6/49 ) NTIME* < COL* J) * J«l*9 ) 

GO T0 51 

40 WRITE ( 6/ 50 > NTIME*UNITS* < ALPHA! J ) * COL t J ) * J-l* K ) 

41 FORMAT* 55X* I5*2X*E13*6) 

42 FORMAT ( 4SX* I5*1X* 2( 1X*E13»6 ) ) 

43 FORMAT^!/* 15* 1 X* 3 ( IX* E13 . 6 ) ) 

44 F0RMAT{34X*I5*1X*4(1X*E13»6) ) 

45 F0RMAT(27X*I5*1X*5*1X*£13.6) ) 

46 FORMAT ( 20Xj I5*1X*6I1X*E13»6 ) ) 

47 F0RMAT(13X*I5*1X*7(1X*E13.6) I 

48 F 0RMA T ( 6X * I5*1X*8(1X*E13*6> ) 

49 FORMAT* I5*1X*9( 1X*E13»6 ) } 

50 F0RMAT<///6OX*I5* 1X*A4//5(4X* A4*2H= *El3.6*3X>) 

51 NTIMEP=NTIME+1 
IF*N.LT»NN) G0 T0 53 
WRITE(6*52) N*GR0SSU 

52 F0RMAT* I4*1X*A4) 

NN=N+1 

53 IF*NTIME.LT»NTIMEC) 60 TO 65 

54 READ ( 5* 400 ) NTI MEC* CUN I TS* SYMBOL* CVALUE 

400 FORMAT * I5*1X* A4 * A4 *El3.6) 

IF ( SYMBOL ♦ EG • CUNI TS ) GO TO 66 
IF* CUNITS.NE. BLANK) GO TO 59 
00 55 MN= 1 * 400 

IF* SYMBOL. EGUTITLE(MN) ) GO TO 57 

55 CONTINUE 

WRITE *6*56) SYMBOL 

56 F6RMATI/26X* 'THE VARIABLE 'A4*'WHICH YOU WANT TO CHANGE » 

¥ 'DOES NOT MATCH ANY EXISTING VARIABLE*') 

GO TO 66 

57 WRITE I 6* 58 ) NTI ME* UNI TS* SYMBOL* A ( MN ) * CVALUE 

58 F0RMAT(/16X* 'AT«I5*1X*A4* • INTO THE SIMULATION* THE VALUE OF ' 

* A4* • WAS CHANGED FROM *E13«6*' TO ' E 1 3 • 6* ' . * / ) , 

A*MN)»CVALUE 
GO TO 54 

59 IF *UNITS.LQ.CUNITS) GO TP 65 
WRITE* 6*60) UNITS*CUNITS 

60 FORMAT* //34X* 'THE TIME INCREMENT FOR OUTPUT HAS BEEN ' 

¥ 'CHANGED FROM 'A4>* TO 'A4*'.»//l 

UNI TS=CUN ITS 

IF*UNITS.EQ»SECS) GO TO 61 
IF* UNITS. EQ. THIN) GO TO 62 
IF* UNITS. EQ. HOUR) GO TO 63 
NTIMLP*T/1440* + 1 . 

GO TO 65 

61 NTIMEP=T¥60*+1. 

N*T 

GO TO 64 

62 NTI MEP*T+ 1 • 

N*T/60* 

GO TO 64 

63 NTIMEP»T/60* + l . 

N*T /1 440 » 

64 NN«N+1 

65 RETURN 

66 STOP 
END 
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0 • 300000E-03 

anw 

tx 

0. OOOOOOE+OO 


Of 

995578E + 00 

AUG 

Of lOOOOOE+Ol 

AUR 

ft 

0 

988731E+00 

AUS 

0,1000005+01 










THE TIME INCREMENT FSR OUTPUT HAS BEEN 'CHANGEO FROM SECS TO HOUR ♦ 


1 HOUR 


T 

ae 

Of 6520885+02 

I 

n 

0*1141 05E+02 

VBD 

s 

*0* 232339E-03 

VVS 

■ 

0*3274355+01 

VPA 

* 

0*3794485+00 

VAS 

k 

0*8489025+00 

VLA 

ft 

Ot 400638E+00 

VRA 

ft 

0. 100293E+00 

VAE 


0 f 353902E+ 00 

PA 


0 * 996907E + 02 

PAM 

* 

Of 100310E+01 

L VM 

c 

0 f 99 1 1 1 6E+00 

VRE 

s 

Of 292838E-03 

PRA 

■ 

0 f 5856755-01 

GRN 

(* 

0*514056E+01 

VPE 

■ 

0* 731975E-01 

PPA 

s : 

0. 1524955+02 


ft 

0 • 3964 87E + 00 


b 

0 f OOOOOOE+OO 

RP A 

ft 

0» 158813E+01 

RVM 

s 

0*991 1 335+00 

VLE 

E 

0*> 'JB068E-03 

PLA 

ft 

0. 638068E-01 

GUN 

* 

0*51 4252E+01 


s 

0. OOOOOOE+OO 

A1B 

ft 

0* 1013615+01 

RPV 

ft 

0*1396115+01 

RPT 

ft 

0* 298424E + 0 1 

pgl 

a 

0.1518575+02 

GP0 

ze 

0.509045E+01 


» 

0 • OOOOOOE+OO 

VVE 

ft 

Of 323430E+00 


ft 

0*3119 1 8E+00 

PVS 

ft 

0* 378082E+01 


ft 

0* OOOOOOE+OO 

RVG 

h 

0*724181 E+0O 

ove 

ft 

0 * 5099 1 3E+0 1 

AVE 

ft 

0 • 99 74 94E + 00 

CN2 

■ 

0 • ?1 2000E-01 

CN3 

ft 

0#366472E+00 



RVS 

0 

277737E+01 


a 

0 

959099E+02 

RTP 


0* 1 9 47 84E + 02 

QA0 

0* 51 1500E + 01 

QRO 

a 

0 

509459E+01 

OfcO 

0 

509106E+Q1 

D VS 

a 

0 

lb8691E-01 

DPA 


0 * 4 1 3704 E-02 

DAS 

"0* 239344E**01 

DLA 

a 

-0 

607491E-03 

DRA 

0 

453568E-02 


9 

0 

996839E+02 

AUC 


0. OOOOOOE+OO 

AUB 

0* 100454E+01 

A UN 

a 

0 

OOOOOOE+OO 

AU 

0 

992931 E+OO 


M 

0 

OQQQOOE+OO 

AU8 


* 0 • 353444E-05 

0 AU 

0 • 992931 E+00 

A U J 

a 

0 

988291 E+00 

0 

98829 1 E+QO 

AU0 

a 

-0 

117092E-01 

AUH 


0.996487E+00 


O.OOOOOOE+OO 


a 

0 

OOOOOOE+OO 

AUM 

0 

990047E+00 

AU4 

m 

0 

206422E-01 

V IF 


0*5541 48E+00 

P01 

0 . 825000E+01 

PTT 

a 

0 

101 161E+01 

RTS 

0 

698778E+01 

PIF 

a 

"0 

5976 1 7E+0 1 

CPI 


0 • 1652 12E+02 


0 » 4 1 3031E + 01 

CPP 

a 

0 

701763E+02 

RPC 

0 

28Q705E+02 

PVO 

m 

0 

146493E+02 

PC 


0 • 18430 1E+02 

PCD 

0 * 464957E+00 

VTC 

a 

0 

325470E-02 

PEG 

0 

812717E+00 

VTL 

a 

0 

325087E-02 

vtd 


-0.349137E-05 

VPD 

-0*694443E-04 

DPL 

a 

0 

537005E-01 


0 

OOOOOOE+OO 

DPC 

a 

0 

537444E-01 

DPI 


0* 438876E-04 


O.OOOOOOE+OO 

DLP 

a 

0 

696714E-02 


0 

OOOOOOE+OO 

CHY 

a 

0 

494994E+01 

PRM 


•■0 • 500462E+01 

PGR 

0. 197997E+01 

CPG 

a 

0 

124577E+02 

P«F 

0 

413084E+01 


a 

0 

OOOOOOE+OO 

PGX 


0* 165234E+02 

PGC 

0.61 1082E+01 

PGH 

■a 

-0 

399301 E+01 


0 

OOOOOOE+OO 

VG0 

B 

-0 

530434E-04 

VG 


0* 1 1 5 1 52E + 02 


0 . OOOOOOE+OO 


a 

0 

OOOOOOE+OO 


0 

QOOOOOE+OO 

GPD 

8 

•0 

262686E-04 



0 • 31 9625E + 02 

RR 

0 ■ 84 61 96E+02 


a 

0 

1 1 78 10E+01 

APO 

0 

376552E+02 

GLP 

a 

0 

620355E+02 

PFL 


0 ♦ 159650E+02 

GFR 

0* 125106E+00 

TRR 

a 

0 

124075E+00 

VUD 

0 

103092E-02 

REK 

9 

0 

100000E+0 1 

N0D 


0* 103743E+00 

NED 

-0*374323E-02 

NAE 

tx 

0 

213598E+04 

VEC 

0 

150419E+02 

CKE 

n 

0 

498939E+01 

K0D 


0.278986E-02 


O.OOOOOOE+OO 

KIR 

a 

0 

354851E+04 

K IE 

-.0 

1342772-01 

KCD 

9 

-0 

174560E-03 

KED 


0* 164702E-03 

CKI 

0 • 1 42007E+03 

CNA 

a 

0 

142002E+03 

CCD 

0 

5081 18E-02 

V10 

m 

0 

5081 18E-04 

KE 


0.750520E+02 

KI 

0 • 354853E+04 

VIC 

m 

0 

249890E+02 


0 

126200E+02 


9 

0 

OOOOOOE+OO 

Z 


0* 100000E+0 1 


0 ♦ OOOOOOE+OO 


a 

0 

OOOOOOE+OO 

TVZ 

0 

964638E-03 


9 

0 

OOOOOOE+OO 



0. OOOOOOE+OO 

X 

0 ■ 100000E+02 

12 

a 

0 

300000E-02 


0 

585675E-01 

VTS 

9 

0 

120694E+Q2 

VP 


0 • 295957E+01 

PRP 

0*207700E+03 

IFP 

a 

0 

915521 E+01 

QPR 

0 

143454E+03 


9 

0 

OOOOOOE+OO 



0 • OOOOOOE+OO 

AMR 

0 . 98 1927E+00 

AMP 

a 

0 

101547E+01 

AMI 

0 

999865E+00 

AMC 

9 

0 

9942 82E+00 



0 * OOOOOOE+ 00 


O.OOOOOOE+OO 


a 

0 

OOOOOOE+OO 

Ab 

0 

993988E+O0 

CNE 

9 

0 

999338E+01 

AGK 


0 • 200000E+00 

ANP 

0.9961 91 E+00 

A N1 

a 

0 

996191 E+00 

ANC 

0 

995638E+00 


9 

0 

OOOOOOE+OO 



0 • OOOOOOE+OO 


O.OOOOOOE+OO 

ANM 

X 

0 

100290E+01 

V« 

0 

500364 E4-0 1 


a 

0 

OOOOOOE+OO 

HM 


0*4084 76E+02 


O.OOOOOOE+OO 

VIE 

9 

0 

155821 E+01 

VIB 

0 

305821 E+01 

VIM 

m 

0 

10 1 930E+0 1 

RC2 


0.U8544E-04 

P02 

0 • 237500E+00 

RKC 

a 

0 

580000E-05 

RC1 

0 

110200E-04 

RCD 

m 

-0 

834421E-06 

VRC 

a 

0 * 20 4386E+01 

RSN 

0.325487E+02 

OVA 

a 

0 

202968E+03 

8FN 

0 

294666E+0L 

DOB 

m 

0 

1 796 15E+03 

A0M 


0 • 998244E+00 

pie 

0 * 800000E+01 

0SV 

a 

0 

6955O0E+0O 

POT 

0 

821676E+Q1 


m 

-0 

2591 40E+00 

P0B 


0*9844522+00 

ARl 

0.9844 52E + 00 

AR2 

a 

0 

923648E+00 


0 

97 4345E+00 

AR3 

9 

0 

977264E+00 

arm 


0 * 889430E + 00 


0 ♦ 300662E+0 1 

gfn 

a 

0 

125106E+00 


0 

OOOOOOE+OO 


9 

0 

OOOOOOE+OO 

AH 


0.301409E+01 

AHC 

0.100451E+01 


a 

0 

OOOOOOE+OO 


0 

OOQOOOE+OO 


9 

0 

OOOOOOE+OO 

ahm 


0. 996464E+O0 

CNY 

0 • 600000E+01 

CNX 

a 

0 

250000E+01 


0 

OOOOOOE+OO 


9 

0 

OOOOOOE+OO 



O.OOOOOOE+OO 

VV6 

-0.297081E-03 

VV7 

a 

0 

114378E-01 

TVO 

0 

964638E-03 

VTW 

9 

0 

400301E+02 

HSR 


0* 100000E+01 

HSL 

0.100000E+01 

NID 

9 

0 

IOOOOOE+OO 

SR 

0 

500000E+00 

V VR 

a 

0 

295146E+01 

RAR 


0.305200E+02 

cv 

0.825000E-01 

CN7 

9 

0 

200000E+00 

Agx 

0 

300000E+0 1 

AUK 

9 

0 

500000E-03 

AUZ 


0* 100000E+01 

Y 

0 • 1 OOOOOE + 0 1 

CFC 

a 

0 

700000E-02 

CPK 

0 

160000E-06 

PCE 

* 

0 

300000E+0 1 

CPR 

X 

0 • 850000E+02 

LPK 

0.470000Ea03 

DP0 

a 

0 

700000E-02 

HYL 

0 

570000E+02 

KID 

a 

0 

2B0000E-02 

AMT 


0* 600000E + 02 

ANT 

0 ♦ 150000E+02 

POK 

a 

0 

600000E-01 

PON 

0 

300000E+00 

AlK 

9 

0 

lOOOOOE+01 

A2K 

« 

0* 200000E + 02 

A 3K 

0 . 1 1 5200E+05 

CNR 

ts 

0 

139000E+03 

CNZ 

0 

100000E+01 

AHK 

9 

0 

700000E+01 

SRK 


0* 330000E + 02 


O.OOOOOOE+OO 

V 2D 

X 

0 

200000E-01 

Z1 

0 

IQOOOOE+Ol 


e 

0 

OOOOOOE+OO 

Z3 


0 * 400000E + Q 1 

Z 4 

0 . 100000E+02 

Z5 

9 

0 

ioooooe +02 

26 

0 

500O00E+01 

Z 7 

a 

0 

500000E+01 

Z 8 


0. 100000E+01 

HMK 

0.900000E+02 

hkm 

9 

0 

533330E+00 


0 

397409E+02 

P0Z 

9 

0 

300000E+ 00 



0*55 47 83E + 03 

002 

0 « 24 6749E + 04 

RBF 

a 

0 

11 78 10E+01 

M02 

0 

179684 E + 03 

P0A 

n 

0 

922258 E+00 

P0Y 


0 * 464000E-04 

ANU 

0. 100276E+01 

P0R 

■ 

0 

400000E+02 

GE2 

0 

500000E-01 

HMD 

a 

0 

lOOOOOE+O 1 

OHM 


0*S54189E-02 


0* 800000E+01 

13 

a 

0 

200000E+02 

U 

0 

400000E+01 

VP1 

a 

0 

100000E-01 

T1 


0 • 65 2088E + 02 


0* 1005 13E+01 

GF4 

a 

0 

500000E+01 

A UP 

0 

98829 i E+00 

AUV 

n 

0 

300000E+00 



O.OOOOOOE+OO 

AUY 

0 . 250000E+00 

OUT 

a 

0 

300000E+01 

DSP 

0 

300000E+01 

A HZ 

n 

0 

117135E-01 

A HY 


0* 191842E-01 

0SA 

0.9937 84E+00 


a 

-0 

100652E+0? 

CPN 

0 

304262E + 02 

POS 

B 

0 

12 1705E+02 

PLF 


0.280428E-03 

PP0 

0.853236E-02 

PPN 

c 

0 

894377E-02 

PPD 

0 

41 U02E-03 

PFI 

a 

0 

3 1 87 73 E*» 03 

dfp 


0 * 38 34 49E-04 

VPF 

0.128699E-01 

RPR 

a 

0 

382965E+00 

PMC 

0 

68 3 1 46E + 0 1 

PMS 

V 

0 

722800E + 0 1 

PMP 


0 . 45 43 73E + 0 1 

HR 

0*71 64 88E + 02 

CPF 

a 

0 

3000002^03 

PCP 

0 

689736E+01 

DAI 

a 

0 

lOOOOOF+O 1 

DLZ 

n 

0 « 6967 14E*"02 


O.OOOOOOE+OO 


a 

0 

OOOOOOE+OO 


0 

000000E+ 00 

N0Z 

a 

0 

103743E+00 


9 

O.OOOOOOE+OO 

B 

O.OOOOOOE+OO 

hpr 

a 

0 

1 00509E + 0 1 

HPL 

0 

100241E+01 

sth 

K 

0 

Iv^OOOE+Ol 

AL0 

B 

0 . ioooooe +01 

EXC «= 

0 • 100000E+01 

02M 

a 

0 

180000E+03 

02A 

0 

100045E+03 


B 

0 

153036E+02 

$VB 

a 

0.7105S8E-01 

AUL « 

0 » 2 1 OOOOE+OO 

VV9 

B 

0 

31 5900E+01 

0 

150000E+00 


* 

0 

OOOOOOE+OO 

EXE 

B 

0.679493F-02 

ARF * 

0. 150000E+01 

ORF 

B 

0 

600000E+00 

RSM 

0 

968703E+02 

BFM 

B 

0 

990085E+00 

RAM 

as 

0 ♦ 963000E+02 

0 VS « 

0 • 698000E + 00 


tt 

0 

398837E+02 


&ftfclNAl PA& 1S$ 
OF POOR QUALITY 



RHO 

■ 

0* 598809E+02 

P&0 

at 

-0* 116302E+00 

OHM 

m 

0*600000E+02 


m 

0*OOOOOOE+O0 

Z10 

u 

0* 825000E + Q1 


m 

0 *OOOOOOE+O0 

STA 

m 

o.ooooonE+oo 


m 

0* OOOOOOE+OO 

AUR 

m 

0* 988291 E + 00 


Q0M 

* 

0. 239996E+04 


at 

0.990696E+00 

PHI 

E 

0* 799709E + 01 


u 

o. ooooooe+oi 

211 

K 

0* 400000E+Q1 


E= 

0 • OOOOOOE+OO 

PAR 

m 

0*996907E+02 

anv 

m 

0*3000002-03 

A US 

m 

0* 1000002+01 


PH0 

m 

0* 799709E+Q 1 

AMM 

M 

0* 990696 £ + 00 

PM3 

m 

0* 1 OOOOCE-02 

FM5 

■m 

0. 122C00E+03 

Z 12 

m 

0. 124000L+01 

PK2 

m 

0.800000E+03 

GBL 

m 

0 • OOOOOOE+OO 

ANVv 

m 

0 • OOOOOOE+OO 


P20 « 0*7997 09 £+01 

A4K « 0* 1 OOOOOE + 0 1 

ph4 « -o* loooooe+oi 

PKl * 0 • 250000E+04 

713 * 0 • 625000E+00 

PK3 ■ 0 • 20Q000E+01 

ANY ■ -0* 200000E+00 
» 0 ♦ 996191 E+00 


MHO 

M 

0 • 598946E+02 

POh 

a 

0 • 800000E- 01 

E X 1 

M 

0* 300000E + 01 


m 

0 • OOOOOOE+OO 


m 

0* OOOOOOE+OO 

FIS 

m 

0* OOOOOOE + OO 

ANZ 

M 

0* 400000E+00 

AUO 

a 

0* i0OO00E+01 




VERIFICATION PLAN AND PROCEDURE FOR 


WHITE'S VERSION, OF GUYTON'S MODEL 


The input-output subroutines are designed to accept 
data from punched cards and display desired output on 132- 
column listing paper. 

The user should first determine exactly what experiment 
he wishes to do. Thus he must decide what parameters he wishes 
to monitor (as many as twenty) and how often the values of 
these (dependent) variables are needed in the output (e.g. 
each simulated second, minute, hour, or day.) The program is 
flexible enough to allow the user to change the frequency of 
output at any predetermined time or times. Then the user 
must decide what independent variables are to be changed, at 
what time or times these changes are to take place in the simulatic 
and what the new values of these variables should be. There 
is no limit to the number of variables which may be changed 
at any given time, nor is there a limit to the number of 
times changes may occur. Finally, he must decide at what 
time the experiment is to be terminated. 

The input cards should be arranged as follows: 

STEP 1 . 

For each of the nearly 400 variables , a card should be 
read giving the initial value, array number, and symbol of 
that particular variable. The initial value should appear 
in the first 13 columns in El 3. 6 notation. Thus the decimal 
point should be in column 3 followed by six digits, the letter 
"E" m column 10, and the exponent right- justified to column 13. 
The array number of the variable being initialized should appear 
as an integer right- justified in columns 18-20, and the symbol 
for that variable should appear left- justified in columns 23-26. 
These cards should be read in one after another, one variable 
per card. A blank card should follow the last variable 
initialized. 


- ‘ ASSSHM. ttAfiR #QT 
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STEP 2 . 

Following this blank card, one card should be read containing 
the symbol (s) for the variable (s) to be monitored. The order 
in which the values of these variable (s) appear as output will 
be the same as the order in which the sumbol(s) for these 
variables are read in at this step. The symbol for the 
first variable desired in output should appear left- justified 
in columns 1-4, the symbol for the second variable left- justified 
in columns 5-8, the symbol for the third variable left-justified 
m columns 9-12, and so on, each symbol left- justified in a 
field of four columns, with a maximum of 80/4=20 symbols. 

For best results, it is suggested that no more than nine 
variables be monitored at a time so that the output appears 
in nice column form. If the values of all variables are 
desired, simply punch "ALL" in columns 1-3 of this card. 

STEP 3 . 

The next card read should contain the time at which the 
first {or next) change of independent variable is to be made, 
in the units of time that the user desires the output to 
appear until that change is made. For example, "8 HOUR" 
would cause the output to apxiear each hour up thru 8 hours, then 
a change of variable (s) (or a change in time units for output) 
would be expected. The time should appear as an integer 
right- justified in columns 1-5 and the units in either "SECS", 
"MINS", "HOUR", or "DAYS" m columns 7-10. 

STEP 4 . . 

Following the time card there should be the card or cards 
which change the values of the desired independent variable (s). 

For each variable to be changed, one card should be read with 
the symbol of that variable left- justified in columns 11-14 
and the new value of that variable in E13.6 format in columns 
; 15-27, with the decimal in column 17 followed by six digits 
and an "E” m column 24 with the exponent right- justified 
in column 27. The user may change the values of as many 
independent variables at this time as he wishes, one change 
of variable per card, according to the instructions above, 
one card after another. 

Steps 3 and 4 may now be repeated as often as desired; 
step 3 to give the time at which the next change is to occur 
and the units of time for output until that change occurs, 
and step 4 to make the desired changes. 

To terminate the experiment at a predetermined time, 
read this time in according to the format given in step 3 
and follow this "termination time” card with a blank. 


03 



EXERCISE STRESS EXPERIMENT 


CARDS 1-377: (Variable initializing cards) 



1 T 


377 AUS 


CARD 378: 

— (BLANK) — 


CARD 379: 

VUD PVO PMO PA AUP QLP BFM MMO 

CARD 380: 

30 SECS 

CARDS 381-386: 


EXC 

6 . 000000E 

01 

A4K 

0 . 025000E 

00 

Z 

5 . OOOOOOE 

00 

Z8 

3. 000000E 

00 

Z5 

1. OOOOOOE 

00 

Z 6 

1. OOOOOOE 

01 


CARD 387: 


120 SECS 



CARD 388 


EXC 1. 000000 OE 00 


CARD 389: 

300 SECS 

CARD 390: 

13 2 . OOOOOOE 01 


CARD 391: 

10 MINS. 

CARD 392: 

— (BLANK) — 


Note: For this experiment, the initial value of 13 is 0. 



NEPHROSIS EXPERIMENT 

CARDS 1-377 (Initializing cards) 



— (BLANK) — 

CARD 379: 

VUD VG VTS VP PRP PIF PA QLO 

CARD 380: 

24 HOUR 

Z 

CARD 381: 

DPO . 050000E 00 

CARD 382: 

128 HOUR 

CARD 383: 

DPO . 021000E 00 



6 . 


CARD 384: 

312 HOUR 

CARD 385: 

— (BLANK) — 



